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SUMMARY 


Test data are given, in duplicate, for approximately 110 dis- 
tinct tension joints that were formed by uniting 24S-T sheets 
with protruding-head rivets. With the aid of a general empirical 
“stress-concentration” factor, the data are well correlated with a 
modified form of the elementary theory of tension joints. The 
modification lies in the assumption that rivets of unequal diam- 
eter in a joint carry equal bearing stresses, rather than equal 
shearing stresses, when the ultimate load is approached. Ex- 
tensive data are given to demonstrate that this is a good approxi- 
mation for the case of 17S-T rivets in 24S-T sheets. An algebraic 
formulation of the elementary theory is presented for the purpose 
of rendering the problem of optimum proportions of a multiple- 
row joint accessible to mathematical analysis. Explicit rules are 
developed for laying out efficient joints. 


1. INTRODUCTION 


. MAIN STEPS in the stress analysis of a riveted 
tension joint are the checking of the rivets and the 
checking of the sheet. 

The first step is usually performed by adding the 
ultimate loads of the respective rivets and ascertaining 
if the sum is greater than the design load for the joint. 
The data given herein show that this is not a safe pro- 
cedure when mixed sizes of rivets occur in a joint, be- 
cause the smallest rivets are apt to shear prematurely. 
Other investigators’ ? who have considered load dis- 
tributions in riveted joints have confined their studies to 
joints with rivets of a single size. Their results have 
shown that the front row of rivets carries a dispropor- 
tionate share of the load when the deformations are 
elastic, but that the load distribution approaches uni- 





that joints be designed with the hypothesis that the 
bearing stresses on all rivets are equal. 

Extensive data on the allowable loads for 24S-T and 
A17S-T rivets at three diameters edge distance have 
been given by G. Holback.* The present investigation 
supplements Holback’s results with data for 17S-T 
protruding-head rivets at two diameters edge distance. 

The checking of the sheet requires experimental in- 
formation on the influences of strain concentrations at 
rivet holes. Schnitt and Koegler‘ have studied the net 
effects of perforations on the static strengths of 24S-T 
sheets. Further studies of these effects have been made 
by Holleman,® Drake,® and many others. However, 
the published data on tests of high-strength aluminum- 
alloy joints are meager. The most complete account of 
such tests has been given by Holleman,® who has pre- 
sented data for an assortment of 20 distinct joint speci- 
mens. The sparsity of joint test data indicates that 
little attention has been directed to the effect of bearing 
stress upon strain concentrations at rivet holes. The 
test data presented herein show that this effect fre- 
quently acts adversely upon a joint. Consequently, 
the interaction between bearing stress and rivet-hole 
strain concentration must be considered if accurate 
stress analyses of joints are to be obtained. 


2. NOMENCLATURE 


In the following exposition of notations and terminology, 
attention is focused on the thinner of the sheets in a lap splice 
that joins two sheets of unequal thickness. Rivet rows, beginning 







formi ‘ . with the rear row, are consecutively numbered 1, 2, 3,..., ” 
a, Mi ae progresses. However, when the (Fig. 1). 
e front row are smaller than those in the back m = margin distance (Fig. 1) 
tows, a rivet failure can occur before the loads have be- e = edge distance (Fig. 1) P 
1 L. | SMe equalized. As a practical method for implicitly d; = initial diameter of rivet holes in the ith row (Fig. 1) 
taking this phenomenon into account, it is proposed ?: = pitch of rivets in the ith row (Fig. 1) 
S = spacing between rows (Fig. 1) 
Received April 29, 1946. Revised and received October 28, W = width of sheet (Fig. 1) 
2pNER & 1946, t = thickness of the sheet 
* Structural Research Engineers. n = number of rows of rivets in the joint 
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Fic. 1. Figure showing nomenclature. 
C; = rivet factor for the ith row (ratio of net area to gross 


area on the row). For a joint of infinite width, this 
becomes 1 — d;/p;. In all cases, Co = 0 
o = uniform tensile stress at some distance from the joint 
o% = average tensile stress in the sheet on the ith row of 
rivets (based on net cross-section area of the sheet 
on the ith row) 


Ii = bearing stress on rivets in the ith row 
ti = ultimate value of f; (see Sect. 6) 
Fi, = ultimate tensile stress of the sheet stock. For ascer- 


taining efficiencies of test specimens, this has been 
determined by test coupons cut from the parent 
stock of the joint test sheets. For more general 
calculations, the value Fi. = 70,000 Ibs. per sq.in. 
has been accepted as typical for unclad 24S-T (the 
ANC-5 value is 66,000 Ibs. per sq.in.) 

Fy, = ultimate bearing stress of the sheet stock (see Sect. 6) 

Fy,, = ultimate shearing stress of rivets in the ith row. This 
can vary from row to row, if, for example, steel and 
aluminum rivets are used in the same joint. In this 
study, the value F,.; = 38,000 Ibs. per sq.in. has 
been accepted as typical for 17S-T protruding-head 
rivets, when rivet cross-section areas are based on 
initial hole diameters (the ANC-5 value is 34,000 
Ibs. per sq.in.) 

A ” fn/ Fru 


E = joint efficiency = (¢/ Fy) maz: 

E; = theoretical efficiency of an idealized joint 

Ro = empirical efficiency correction factor for a perforated 
sheet 

R = empirical efficiency correction factor for a joint 

K =R/R, 

vi = (1 — G)/(1 — CG), ro = 1/11 — Ch), tn = 1 
n 

H =)" 


i=] 
EXPERIMENTAL INVESTIGATION 


3. Pitot TEstTs 


In order to obtain the information required for out- 
lining a suitable test program, a few pilot tests were 
conducted. These are described below. 


3.1 Principle of Similitude 


The ultimate stress for a structure derives its dimen- 
sions from characteristic stresses of the material (e.g., 
Young’s modulus, yield stress, ultimate tensile or shear 
stress, etc.) and, therefore, the linear measurements of 
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the structure enter as dimensionless ratios in the for. 
mulas for ultimate stress. This fact permits a reduction 
in the number of test specimens required to obtain g 
complete set of data, because all ratios of dimensions 
can be.varied while one dimension (e.g., rivet diameter) 
is held constant. However, for the study of riveted 
joints, the validity of this principle may be questioned, 
because the process of rolling metal sheet may cause 
differences in crystal structure among sheets of unequal 
thickness. This, in turn, may influence strain concen. 
trations at rivet holes. Consequently, pilot tests were 
devised to exhibit scale effects. These tests were per- 
formed on tandem-riveted lap joints with two identical 
rows of rivets. : The test results are given in Table J, 
These data show that the principle of similitude holds 
for riveted joints in 24S-T sheet; i.e., joints with the 
same ratios of dimensions (d/t, p/d, and S/p) furnish 
approximately the same efficiencies, irrespective of their 
sizes. 


3.2 Effect of Sheet Thickness 


The data of Table I show that sheet thickness has no 
appreciable effect if all ratios of dimensions are held 
constant. However, there is reason to believe that 
sheet thickness is an unimportant dimension, even 
among joint specimens that are not geometrically simi- 
lar. For, since a sheet is subjected to plane stress ¢, 
and since the ratio of this stress to the rivet bearing 
stress is independent of the sheet thickness, it appears 
that the sheet thickness has no influence upon the effi- 
ciency of a joint unless the rivets become critical in 
shear. This argument, of course, implies that the ulti- 
mate bearing stress for a sheet is independent of the 
thickness—an assumption that is not thoroughly sub- 
stantiated by test results (see Fig. 2). 

TABLE I 
YANDEM-RIVETED LAP JOINTS WITH TWO IDEWTIGAl ROWS OF RIVETS 
EFFECT OF SCALE FACTOR ON EFFICIENCY 


% cr 
SPEC, WIDTH HICK- BOLE DIA. PITCH sPac- w/p FL p/é a/t s/p = wr 
». i 1m 
oe 
67.8 
2.068 .0219 131.132 2375 31S 3/4 71300 3S 28/4 5/6 2165 
D 2,087 0212 .152 Sis2 (375 313 3/4 «71500 Ss 8/4 5/6 2140 BBD 
60 Web 
2a 40136 20404 0256 0256 4750 6625 3/4 70100 3S 25/4 5/6 8&2 
Z 4:13 0398 1256 .254 .780 .625 3/4 70100 3 25/4 5/6 6270 8 
900 (LF 
8,221 .0621 .800 .800 1,500 1.250 3/4 72000 3 25/4 5/6 3 
bo 8.236 20848 5300 500 1,800 1,250 3/4 72000 3 25/4 5/6 [5420 0.6 
ns 
fe 2.406 20214 «6132 132.438 «4313 «3/4 2S 72800 7/2 25/4 5/7 2625 
f> 2,406 0210 132,132,438 .313 3/4 71300 7/2 25/4 5/7 2640 Set 
13.8 
Sa 4,817 0596 284.254 .875 .625 3/6 70100 7/2 28/4 5/7 9850 s 
Sb 4,839 $0400 $254 254 875 .625 3/6 70100 7/2 28/4 8/7 9990 7S! 
72.4 
6a 9.596 .0815 500 ,500 1.750 1,250 3/4 72000 7/2 28/4 5/7 40800 
6> 9.629 50848 800.500 1.750 1.250 3/4 72000 7/2 25/4 5/7 40740 6%? 
7a 2.749 .02146 ,131.131 .500 .315 3/4 71300 4 28/4 5/8 3090 = 
7 2.756 0209 .131 131 .500 .313 3/4 71300 4 25/4 5/8 2990 Te 
Ba 5.463 .0898 254 .254 1,000 .625 3/4 70100 4 25/4 5/8 11200 be 
Bb 5.501 .0397 254 .254 1,000 .625 3/4 70100 4 25/4 8/8 11200 5 
9a 10,996 .0635 500 .500 2,000 1.250 3/4 72000 4 25/4 5/8 46900 af 
9 11,010 ,0620 .500 .800 2,000 1,250 3/4 72000 4 25/4 5/8 47100 


7300 9/2 25/4 5/9 2735 Theé 
71300 9/2 25/4 5/9 2550 673 


70100 9/2 25/4 5/9 10425 ret 
70100 9/2 25/4 5/9 10280 75+ 


72000 25/4 5/9 41660 69 
72000 * af 8/9 41400 «670 


10a 2.529 .0212 .131.151 .563 .313 3/4 
10 2.530 0210 .135 .135 .563 .313 3/4 


2254 .254 1.125 .625 3/4 
0254 4254 1.125 .625 3/4 


2500 .500 2,250 1,250 3/4 
2500 500 2,250 1,250 3/4 


lia 5.087 .0403 
ll 5.079 .0595 


12a 10,130 ,0618 
12> 10,137 0647 
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RATIO -OF RIVET DIAMETER TO SHEET THICKNESS 


Fic. 2. Chart showing maximum bearing stress for 17S-T pro- 
truding-head rivets in unclad 24S-T sheet. Crushing lines apply 
for joints with edge distance equal to 2d,. To obtain the ultimate 
bearing stress for rivets of the 7th row of a joint, select the smaller 
of the ordinates corresponding to the abscissas d,/t and d;/t. 





TABLE IT 
TANDEM-RIVETED LAP JOINTS WITH TwO IDENTICAL ROWS OF RIVETS. 
EFFECT OF SHEET THICKNESS ON EFFICIENCY 


5,501 .0397 .254 .254 1,000 .625 3/4 70100 11200 731 


sm, WIDTH THICK- MOLE DIA. PITCH SPAC- w/p F UL?. EFF‘CY APPROX. 
© mss | mG en 7 Y":) BEARING 

sy STRESS AT 

FAILURE 
le 4122 0210 .257 .257 .750 .625 3/4 71500 4555 70.2 80000 
T 4.123 .0209 4257 .257 .750 .6285 3/4 713500 4525 70,4 80000 
mm «4.136 0404 4256 4256 4750 4625 3/4 70100 8260 70.6 80000 
4.135 0598 .254 .254 .750 “.625 3/4 70100 8270 71.8 80000 
Se 4123 0812 9=257 2257 «4750 «= 6625S 3/4 72000 «16625 «= 6849 += 80000 
% 4122 ,0806 .258 .258 .750 .625 3/4 72000 16975 70.9 80000 
4 «4,807 .0210 4257 4257 4875 .625 3/4 71500 5100 70.9 97500 
@ 4,802 .0212 .257 .257 .875 .625 3/4 71300 S040 69.5 97500 
Se 4,817 0598 4254 4254 .875 .625 3/4 70100 9850 75.3 97500 
S 4.839 .0400 .254 .254 .875 .625 3/4 70100 9990 75.7 97500 
f 5495 .0212 .257 .257 1,000 .625 3/4 71300 5660 68,1 111000 
& 5.500 .0210 .256 .256 1,000 .625 3/4 71500 5940 72.2 111000 
fe «5465 40598 4254 2254 1,000 .625 3/4 70100 11200 73.4 111000 
nr 111000 


| 
| 
| 


In order to check the effect of sheet thickness upon 
joint efficiency, tests were performed on several sets of 
joints for which the sheet thickness was the only vari- 
able in a set. Data from these tests are given in Table 
II. 

The data of Table II indicate that the efficiency of a 
48-T joint does not depend appreciably upon the sheet 
thickness, if the rivets do not become critical in shear 
and if the bearing stresses do not exceed bounds that are 
determined by the maximum ratio d/t. For d/t = 12.5, 
itis apparent that this bound lies between 80,000 and 
100,000 Ibs. per sq.in. However, the ratio d/t = 12.5 
is exceptionally high, and it appears likely that when 
d/tis kept within practical limits (say d/t <6), sheet 
thickness will be an unimportant dimension when bear- 
ing stresses do not exceed ANC-5 allowables. 


33 Effect of Margin Distance 


Several specimens of Types A, B, and C (Fig. 3) were 
tested to determine the effect of margin distance on 
joint efficiency. Results of these tests are given in 
Table ITT. 





These data indicate that tandem-riveted joints bene- 
fit by large margin distance if the bearing stresses are 
low, but that the efficiency may decrease as margin 
distance increases if the bearing stresses are high. The 
influence of margin distance decreases as the overall 
width increases. The lesser effect of margin distance on 
wide specimens might be intuitively expected, and it is 
exhibited by Specimens Nes. 3, 6, 7, and 8 of Table III. 
However, margin distance nas a pronounced effect on 
joints with staggered rivet patterns. This is shown by 
Joints Nos. 9 and 10. Here the jeints with the larger 
margin distance (No. 10) gave considerably less effi- 
ciency. It was also noticed that Joints Nos. 9 and 10 
experienced dissimilar tears. Joints No.9 tore straight 
across the inner rows of rivets, while Joints No. 10 
sheared out the sheet at the end rivets. The latter 
“‘dog-eared”’ type of tear (Fig. 6) is characteristic of 
joints with a large margin distance and it no doub: re- 
sults when some critical value of the ratio m/e is ex- 
ceeded. It is also accompanied by a distinet drop m 
efficiency, as the data of Table III show. In fact, the 
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Fic. 3. Types of specimens used in pilot tests on margin dis- 
tance. 
TABLE III 


REFERENCE -- FIGURE S 
EFFECT OF MARGIN DISTAICE ON EFFICIENCY 


{ THICK- HOLE DIA, PITCH SPaC- wp F. ULT. EFF'CY APPROX. 


SPEC. TYPE WIDTH 
1» Ness “4 a ING ee) FEARING 
7 2 STRESS A? 
FAILURE 
le A 2.242 4.0390 4282 4253 2750 625 V2 71000 4170 67.1 86000 
lb A = 2e245 0894 *253 1253 750 .625 1/2 71000 4280 68.2 86000 
2a A 20433 40393 «6253 4253) «6750 =. 625 = 8/8 «71000 4810 «70,9 86000 
2b A - 24442 $0391 "253 4283 4750 .625 $/8 71000 4720 69.6 986000 
_ A 2634 0892 «25S 4253 4750 6625 43/4 971000 S210 71.1 86000 
3b A 2.618 ,039§ .253 .253 .750 .625 3/4 71000 5115 69,8 86000 
se A 24.985 90395 2535 .255 1,00 .625 1/2 71000 6345 75-6 117000 
a> A 2.976 20396 $253 $253 1.00 .625 1/2 71000 6230 74-4 117000 
74.5 127000 
Se A 3.255 .0391 .256 .256 1,00 .625 5/8 71000 6700 
sb A 3.248 ,0395 $253 253 1,00 .625 8/8 71000 6725 741 117000 
© 127000 
6a A 3,500 ,0590 .253 .253 1,00 .625 3/4 71000 6980 72. 
6d A 3.496 ,0396 .253 283 1,00 ,625 3/4 72000 7050 72.0 127000 
Te B 4.136 40404 4256 4256 .750 2625 3/4 70100 6260 70,6 80000 
% B 4135 ,0396 .254 .254 .750 .625 3/4 70100 8270 71.8 80000 
Be B 5.465 ,0396 .254 .254 1.00 .625 3/4 70100 11200 73,4 111000 
8 B 5.501 .0397 .254 .254 1,00 .62 3/4 70100 11200 75,2 111000 
Se C 3.761 04602 4257 .257 .750 .625 1/2 70800 7465 70.1 83000 
% C 3.741 0396 .257 .257 750 .685 1/2 70800 7480 71,2 85000 
l0e 6G) 44232 0368 4257 2257 2750 .625 3/4 70100 7390 66.6 83000 
10 C 4,140 .0387 .257 .257 .750 .625 3/6 70100 7400 66,6 83000 
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TALE IV 
REFERENCE — FIGURE 4 
EFFECT OF WIDTH ON EPFICIENCY.OF PERFORATED SHEETS 


SPEC. TYPE #WUMBER WIDTH THICK- PITCH SPAC- F ULT. EFF’cY 
©. OF HOLES NESS ING wu LOAD 
PER ROW 

le Tandem 2 1.237 _,0495 6/8 11/32 69000 2865 67.8 
1b ° 2 1.250 0493 5/8 11f32. 69000 «= 2885 67,8 
2a id 5 3.121 .0510 5/6 1/32 72040 7925 69.0 
2b ° 5 3,120 ,0512 5/6 1/32 72060 7925 68.9 
Se ° 10 6.248 0493 5/6 11/32 69000 15600 73.4 
3b ° 10 60246 = 0494 5/8 11/32. 69000 =: 15600 7323 
40 «= Tandem 2 1.235 0493 5/8 V2 69000 2870 68.4 
4b ° 2 1,243 0492 5/8 V2 69000 2890 68,5 
Se ° 5 3.130 623 5/8 1/2 71100 8340 71.5 
Sb ° 5 3.117 0523 5/8 V2 72100 8280 70.5 
6a ° 10 6.253 0485 5/6 V2 69000 14450 73.8 
6b e 10 6.258 0491 5/68 1/2 69000 15200 71,7 
7a «Tanten 2 1.999 0490 1.00 7/16 69000 5290 78.2 
a) ° 2 1.999 0490 1,00 7/16 69000 5165 1664 
Se ° 5 4.995 0501 1,00 7/16 71400 14100 78.9 
8b ° 5 5.016 0520 1.00 7/16 72030 14720 78.0 
9a e 10 9.970 ,0491 1.00 7/16 69000 26500 78.4 
% e 10 10,010 0695 1.00 7/16 69000 3263200 76.6 
10a Staggered 2 1.548 .0490 3/6 1/32 69000 3540 67.6 
1m * 2 1.543 0492 5/6 13/32 69000 3485 66.5 
lla ° 5 3.449 = -.9520 5/6 11/32 72500 9080 69.9 
ub * 5 3.467 0523 5/6 11/32 71280 8990 69.6 
l2a * 10 6.568  .0495 5/6 11/32 69000 16000 71.3 
> = 10 6.572 0485 8/8 11/32 69000 15300 69.6 
13a Stiggered 2 1.558 0492 5/8 V2 69000 3595 67.9 
is 0 * 2 1.553 40492 3/8 1/2 69000 3695 70.0 
l4e ° 5 3.430 0505 5/8 V2 72200 8800 69.9 
14d ° 5 30449 0522 5/8 1/2 72800 9100 69.8 
15a ° 10 6.573 0496 5/8 1/2 69000 16300 72,5 
sdb * 10 6.570 ,0188 5/8 1/2 69000 16000 72.3 
léa Staggered 2 2.487 0480 1.00 7/16 69000 6105 74.1 
16d « 2 2.500 0462 1.00 7/16 69000 6080 73.1 
178 ° 5 5.505 .,05]3 1,00 7/16 71800 15750 77.7 
17> ° 5 5.500 ,0522 1.00 7/16 71610 15540 75.5 
188 ° 10 10,500 0490 1.00 7/16 69000 26800 75.5 
18> ° 10 10,490  .0490 1.00 7/16 69000 27400 77.3 





following tests repeatedly show that any deviation from 
a single-row terrc on the front line of rivets is a mark of 
inefficiency. Tandem-riveted joints can also experience 
dog-eared tears but they are less susceptible to these 
failures than staggered joints. In laying out the 
specimens for the main tests of this program, 
proportions that would lead to dog-eared tears were 
avoided. 


3.4 Effect of Width 


The question arises, ‘‘How wide must a joint speci- 
men be, in order that it be representative of an infinitely 
wide sheet splice?’ This subject has been studied by 
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rivet holes. All holes were filled with AN-3 bolts 
Both tandem and staggered hole patterns were investi. 
gated. Corresponding specimens with two, five, and 
ten holes in a row were tested. Fig. 4 shows the five. 
hole patterns. 

Results of the tests of the perforated sheet specimens 
are shown in Table IV. These data exhibit a small in. 
crease of efficiency with the specimen width. This js 
believed to be attributable to the fact that margin 
strain concentrations, which cause premature failures, 
detract approximately a constant amount from the 
strengths of the specimens, and that the relative effects 
of these strain concentrations consequently disappear 
asymptotically as the width of a specimen is increased, 

With the exception of the first group of specimens of 
Table IV, there is little difference between the results of 
tests of specimens containing five holes per row and 
those containing ten holes per row. It may conse. 
quently be assumed that patterns containing five holes 
per row are suitable for test purposes, and that the 
efficiency of a five-hole specimen differs from the eff- 
ciency of an infinitely wide specimen by a percentage 
that is no greater than the test scatter in duplicate speci- 
mens. 
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SECTION A-A 


Fic. 5. Specimens used in pilot tests of comparative lap and 
butt joints. 














LAP JOINTS VERSUS BUTT JOINTS 


SPEC. TIPE %. RWS WIDTH THICK- RIVET F uLr, 4 xprct 
®. OF RIVETS Ness DIA. 2 86D 

la Lap 2 6.56 040 2256 70000 11760 «= GAs 
lb . 2 6.58 ,041 0257 68800 11675 629 
2a Butt 1 6.58 061 258 68900 11000 Si? 
2b . 1 6.57 041 2258 68800 11500 610 
3a Butt 2 6.57 .061 +256 68500 11655 Se 
sb . 2 6.58 041 2256 69500 11600 619 
4a Stiffened Butt 1 6.59 041 2285 69600 10990 565 
4d ° * 1 6.58 041 2260 69200 10%0 8657.5 
Se Stiffened Butt 2 6.59 041 2256 68200 11800 40 
1) ° . 2 6.58 041 «286 69900 «11450 G18 

- | ne 





means of tensile ruptures of 0.051-gage unclad 245.7 
sheets that were perforated with two transverse rows of 































35 I 


Lay 
aircra 
are St 
as bul 
inquit 
sentec 
fortur 
sentia 
laps b 
the sti 
curlin: 
to a le 

To. 
and bt 
stiffen 
embra 
with t 
in., an 
the fix 
Result 

It is 
efficien 
conclu 
test pt 
gives n 














4. 


This 
distinc 
and tri 
tested 


Eack 
test shi 
sheet t 
The de 
joints 1 
3.5, in 
equival 
equal g 
of the s 
for thes 
the con: 
type of 
heavier 
test sh 
heavier 
by the 1 
betweer 
Warrant 
occur w 
plate of 
is intro 

The 
Tivet s 
appreci 















d 245-7 
> TOWS of 
3 bolts, 
investi. 
ive, and 
the five. 


eCimens 
small in- 

This is 

margin 
failures, 
om the 
e effects 
sappear 
creased, 
mens of 
sults of 
ow and 
” conse. 
ve holes 
hat the 
the effi- 


centage 
e speci- 


lap and 


wr'ct 
0 
62.9 


oo 99 9398S 89S &¢o ~~ 
3 
te 


ff cf £2 € 


| 








3.5 Lap Joints Versus Butt Joints 


Lap splices and butt splices are both common in 
aircraft structures. Also, splice plates in butt joints 
are sometimes stiffened by adjacent structures, such 
as bulkheads or stringers. Therefore, it is pertinent to 
inquire if all of these types can be adequately repre- 
sented by lap-type specimens. It appears that this 
fortunate circumstance exists, since a butt joint is es- 
sentially the same as two lap joints in tandem—the 
laps being formed with a common splice plate. Also, 
the stiffening of a splice plate merely restrains it against 
curling, and this presumably cannot affect the efficiency 
to a large extent. 

To check these conjectures, several comparable lap 
and butt joints—the latter with both stiffened and un- 
stiffened splice plates—were tested. The specimens 
embraced single- and double-row tandem rivet patterns 
with ten '/,-in. rivets per row. The rivet pitch was °/s 
in., and the margin distance was '/3. in. Side views of 
the five joints that were tested are shown in Fig. 5. 
Results of the tests are given in Table V. 

It is seen that good agreement is obtained among the 
efficiencies of comparable specimens in Table V. It is 
concluded that lap and butt joints are equivalent for 
test purposes, and that the stiffening of a splice plate 
gives no appreciable benefit. 


4. Main TESTS—SPECIMENS AND PROCEDURES 


This section describes tests of approximately 110 
distinct joints representing seven basic single-, double-, 
and triple-row patterns. Most of the specimens were 
tested in duplicate. 

Each test specimen consisted of a lap joint between a 
test sheet, which was to experience rupture, and a jig 
sheet that was one gage thicker than the test sheet. 
The decision to restrict the program to a study of lap 
joints resulted from the pilot tests described in Sect. 
3.5, in which it is shown that lap and butt joints are 
equivalent for test purposes. Although joints with 
equal gage sheets were not explicitly tested, a number 
of the specimens had constant pitch and rivet size, and, 
for these specimens, the difference in thickness between 
the connected sheets is believed to be irrelevant since the 
type of failure would not have been altered had the 
heavier sheet been reduced to the same thickness as the 
test sheet. Also, the stiffening of the joint by the 
heavier sheet is an inconsequential effect, as is shown 
by the tests of Sect. 3.5. An extensive study of splices 
between equal-gage sheets was considered to be un- 
watranted, because important splices in aircraft usually 
“eur where sheet gages are changed or where a splice 
plate of somewhat heavier gage than the thinner sheet 
is introduced. 

The present tests only incidentally include cases of 
tivet shear. Therefore, the joint efficiencies were not 
appreciably dependent upon the sheet thickness (see 
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Sect. 3.2), and all of the test sheets were accordingly 
cut from 0.040-in. or 0.125-in. stock. 

All specimens were made of unclad 24S-T aluminum 
alloy. 17S-T aluminum-alloy flat-head rivets (pressed, 
not driven) were used in most cases. Some small stain- 
less-steel rivets were introduced in joints with mixed 
sizes of rivets when calculations indicated that the 
smallest rivets would shear prematurely if made of 
17S-T. All rivets were inserted with the manufactured 
heads against the test sheets. 

Photographs showing the six basic multiple-row pat- 
terns are presented in Fig. 6. In all cases, except the 
single-row patterns, the edge distance e on the test 
sheet was equal to 2d). 
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Fic. 6. Types of multiple-row patterns and representative 
failures. 


{ The widths of the specimens were conditioned to some 
extent by the rivet patterns, since it was considered 
desirable to have the rivets disposed symmetrically 
about the longitudinal centerline of a joint. It appears 
from the pilot tests of Sect. 3.4 that the specimens could 
advantageously have been built wider, but the efficien- 
cies that have been obtained lie within 1 or 2 per cent of 
those that could have been realized with extremely wide 
specimens. 

A tolerance of +0.015 in. was allowed for the rivet 
pitch, margin distance, edge distance, and row spacing. © 
The rivet holes were 0.003 to 0.007 in. larger than the 
nominal rivet diameters, with a maximum difference of 
0.003 in. in any one specimen. Diameters shown in the 
tables are average hole diameters prior to pressing of 
the rivets. 

Specimens were bolted to '/,-in. steel plates which 
were connected through pivots to bars that were gripped 
in the. jaws of a testing machine. Fig. 7 is a photo- 
graph of a typical test setup. The center pivots in the 
jig plates insured pure tension in the specimens; i.e., 
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Fic. 7. 


Typical test setup. 


transverse bending moments were eliminated. The 
specimens were tested in 200,000- and 400,000-Ib. 
Southwark-Emery machines that were accurate within 1 
per cent in the ranges used. Loads were increased 
slowly until the specimens failed. 

The ultimate tensile stresses of the sheet stock were 
determined from standard '/2-in. tensile coupons. 
Three coupons were cut from each piece of sheet stock. 
The average ultimate tensile stress obtained from these 
coupons was taken to be the F,, of all specimens cut 
from the parent sheet. 

5. Main Tests—RESULTS 

Joints have been arranged in series according to their 
types. The series are designated by a number and a 
letter; thus, Series 2A denotes a two-row joint of type 
A (shown in Fig. 6). Joints in a given series have been 
numbered consecutively, and the letters “‘a’’ and ‘‘b”’ 
have been appended to denote duplicate specimens. 

Tables VI to XII present the joint numbers, the meas- 
ured widths and thicknesses of the test specimens, the 
measured initial rivet hole diameters, the nominal values 
of the pitch and spacing, the numbers of rivets in the 
rear rows, the rivet factors C, of the front rows, the 
ultimate test loads, the F,,, values of the coupons, the 
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measured joint efficiencies, the predicted joint eff. 
ciencies, and references to the charts in which the data 
are plotted. 

The experimental charts of E versus C, which are 
shown in Figs. 8 to 20 have been plotted from the data 
of Tables VI to XII. Each chart covers a particular 
rivet pattern. Theoretical lines of sheet tear, rivet 
shear, and bearing failure for corresponding ideal joints 
(see Sect. 7) are superposed on these charts. 


6. ULTIMATE BEARING STRESSES 


Holback* has shown, by means of a large number of 
test data, that the load F per rivet, for rivet failures, is 
determined by an equation of the form 


F/d? = ¢(d/t) 


in which d/t is the ratio of the rivet hole diameter to the 
sheet thickness. This equation may be written: 


F/dt = ¥(d/2) 


Since the ratio F/dt has the significance of a bearing 
stress, it may be concluded that, for a given joint pat- 
tern, with edge distance and row spacings specified in 
terms of the rivet diameter, there is an ultimate bearing 
stress f which is a function of the ratio d/t. Theo- 
retically, this stress is the smaller of the two bearing 
stresses that cause crushing failures and shear failures; 
i.e., the theoretical ultimate bearing stress f; for the 
rivets of the ith row is the smaller of the values given 
by the following equations: 


fi = F,, and fi = Wd Fyy;/ At (1) 


The present tests have shown remarkably close agree- 
ment with these equations, although, as Holback has 
shown, the stress F,, decreases linear!y as the ratio d/t 
increases. Fig. 2 is a graph of f; versus d/t for multiple- 
row 248-T unclad lap joints with 17S-T flat-head rivets 
at an edge distance e 2d,. A graph, based upon 
Eq. (1) and ANC-5 allowables, is also given. The 
abscissa d,/t corresponds to the shear lines on these 
graphs. However, the abscissa d;/t corresponds to the 
bearing lines, irrespective of what the diameter of rivets 
in the ith row may be. This is because bearing failures 
are precipitated by the first (rear) row. 

The data from which Fig. 2 has heen obtained are 
given in Table VII. It is unnecessary to discuss these 
data at present; the justification of the graph is that it 


TABLE VI 
SERIES | 
SINGLE ROW JOINTS 
JOINT WIDTH THICKNESS RIVET HOLE PITCH NO OF. m/p e/d C, COUPON ULTIMATE OBSERVED TYPE OF CHART COMME NTS 
NUMBER DIAMETER RIVETS Ftu LOAD EFFICIENCY FAILURE REFERENCE 
1 6.75 .0ug - 256 - 750 9 -50 4Y =.658 466,900 11,950 54.0 TENSION FIG, 8 SHEETS OF EQUAL THICKNESS 
ib 6.76 -048 -257 . 750 9 -50 a -667 +900 11,950 56.4 TENSION FIG. 8 SHEETS OF EQUAL THICKNESS 
2a 9.06 -0u9 - 259 1.000 9 -50 yu =.743. «69,200 13,200 43.0 TENSION FIG. 8 SHEETS OF GM) UAL THICKNESS 
2b 9.01 -049 - 256 1.000 9 -50 4 «6.744 «69,500 14,600 47.5 TEINSION FIG. 8 SHEETS OF EQUAL THICKNESS 
3a 13.51 -048 - 257 1.500 9 -50 4 .829 69,800 17,600 39.8 TENSjON FIG. 8 SHEETS OF EQUAL THICKNESS 
3b 13.51 -048 -256 1.500 9 -50 4 = =.830 .69,800 15,000 33.2 TENSION FIG, 8 SHEETS OF EQUAL Lar 
Ma 6.75 -0u9 317 - 750 9 -50 3.2 .577 69,800 12,200 53.6 TENSION FIG. 8 SHEETS OF EQUAL va cee 
Uo 6.75 -0u9 319 - 750 9 -60 3.2 .575 ,»900 12,200 55.2 TENSION FIG. 8 SHEETS OF EQUAL May 1 SPLICE PUT 
5a 6.58 -Oul - 258 625 10 -75 2 .607 68,900 11,000 59.2 TENS1ON FIG. 8 BUTT SPLICE WITH .05 » SPLICE pull 
5b 6.57 -Oul .258 -625 10 -76 Ss. & 68,800 !1!,300 61.0 TENSION FIG. 8 BUTT SPLICE WITH .05! 
Se 
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TABLE Vii 
SERIES 2A 
gut. WIDTH THICKNESS RIVET HOLE PITCH SPACING NO. OF m/p C, COUPON ULTIMATE OBSERVED CALCULATED CHART COMMENTS 
MER DIAMETER RIVETS Ftu LOAD = EFFICIENCY EFFICIENCY REFERENCES 
a; do PER ROW 
ia «63-437—s«wOM0G—“(«é‘i«‘ CCRC 5 3/4 €30 70,106 6300 64.5 €5 0 FIG 9 
ip 3.455 «= .0N00—‘iéwzSOCLO5UOCCEZH «25 5 3/4 632 70,100 6350 65.5 FIG. 9 
tm 4.136 .ONON 1256 .256 .750 .625 5 3/4 .€90 70,100 8260 70.6 70 3 FIG. 9 
» 4.135 0398 «= «254 254 750) 625 5 3/4 .692 70;100 98270 71.8 FIG. 9 
% 4.817 0308 .254 .254 875 625 5 3/4 .736 70.100 9850 73.3 75.€ FIG. 9 
> «4.839 «=. 400 254.254 .875 .€25 5 3/4 .727 70,100 9990 73.7 FIG. 9 
ig 5.463 0398  .254 .254 1.000 625 5 2/4 .768 70,100 11,200 72.4 74.0 FIG. 9 
w 5.501 0397 z54 .254 1.000 .625 . 3/4 .770 70,100 11.200 73.1 FIG. 9 
& 5.087  .O40S 254.254 1.125  .€25 \ 3/4 600 76,100 10,425 72.€ 73 4 FIG. 9 
> 5.079 0395 8.254.254 «1.125. €25 4 3/4 800 70,100 10,280 73.1 FIG. 9 
fa 5.64! 0398 «=s_—«w 25 254 1.250 = 625 y 3/4 820 70,100 10,920 69.5 69.0 FIG. 9 BEARING FAILURE 
f> 5.650 0398 «=s-—«w 25s 251.250). 625 4 3/4 820 70,100 10,800 66.5 FIG. 9 BEARING FAILURE 
ta 6.200 Ou02 «= «« 25 —i«C BN L375) 625 \ 3/4 83€ 70,100 11.475 65.7 63.2 FIG. 9 BEARING FAILURE 
hb 6.209 0398 «6.256 «256 «(1.375 3.625 y 3/4 835 70,100 141,000 65.5 FIG. 9 BEARING FAILURE 
fa 6.740 0400 .25€ .256 1.500 ‘€25 u 3/4 BYE 70,100 11,100 5e.8 FIG. 9 BEARING FAILURE 
6.750 .O0N05 254 .254 1.500 .€25 u 3/4 850 70,5100 11,Nz5 59.6 57 8 FIG. 9 BEARING FAILURE 
9 3. uuy 399 255 .255 .625 1.125 5 3/4 620 70,100  62€0 65.0 65.0 FIG. 9 
0 416 0398 .255 .255 .760 1.125 5 3/4 690 70,100 8160 71.0 70.0 FIG. § 
ii =. 832 0396 -253 .252 875 «1.125 5 3/4 .738 70,100 9990 74 .€ 73.5 FIG. § 
2 §.483 0398 255 .255 1.000 1.125 5 3/K 7€7 70,10C 11,140 72.6 74.2 FIG. 9 
3 5.073 .OM05 .255 .255 1.125 1.125 y =/4 799 70,100 10,375 72.0 73.5 FIG. 9 
ih = 5.647, «ONO? Ss‘ BSC OBN 1.250 «1.125 \ 2/4 .820 70,100 11,025 69.2 €9.4 FIG. 9 BEARING FAILURE 
6 6.193 .OM05  .255 .255 1.375 1.125 ¥ 3/4 .63€ 70,100 11,470 65.2 €3.2 FIG. 9 BEARING FAILURE 
6 6.774% ONO! 254 .254 1.500 1.125 \ 3/4 .850 70,100 11.675 61.3 57.8 FIG. 9 BEARING FAILURE 
ia 5.060 .02I2 257 257 «(1.125 =. 625 y 3/4 .797 71,300 §110 €6.9 NOT PLOTTED BEARING FAILURE 
im 5.059  .0210 257.257 1.125 .625 4 3/4 .797 71,300 5000 66.0 NOT PLOTTED BEARING FAILURE 
ie = 3 -O410 = .257, «191 Ss 625s“ 25 5 1/2 .696 70,800 6450 71.4 71.4 FIG. 10 
i 3. -OUON =s-«. 257 1191 = «6625S 625 5 1/2 .696 70,800 €360 71.0 FIG. IC 
Boi ME Br cl ob CAS eee ene Re He 
; -O40 -257 .191 .750 = .€25— - : ; 4 ; 
te i: "Ou06 «= 257-1191 «875 1625 5 1/2 :782 70,800 9870 78.5 77.€ FIG. 10 
i |. ou0’u 257 191 .875 625 5 1/2 .782 70,800 9730 77.8 FIG. 10 
tla 8 040! .257 .191 1.000 .625 5 1/2 809 70,800 141,200 79.C FIG. 10 
2M. .0U02 0.257 191 «1.000 =. 625 5 1/2 .809 70,600 11,000 77.4 78.0 FIG. 10 
tm 3. 0405 257.191 .625 1.250 a 1/2 .694 70,800 645C 72.2 ; FIG. 10 
mb 3. .0407 257 .191 + .625 1.250 5 1/2 694 70,800 6800 72.4 71.2 FIG. 10. 
m3. .040Z = .257— 191 75C 1.250 5 1/2 .746 70,80C 8100 76.9 a FIG. 10 
a 3. oud! 297 191 780 1.260 8 112 748 70,600 79ec 7s. 75.% FIG. 10 
ls MLE -OU00 =. «. 257 191s 875s. \ 782 70,8 . : 
me y. Ou02  =.257. «191~=s -875—«1. 250 E 1/2 .782 70,800 9930 8C.0 77.7 FIG. 10 
eM. OuCO 8=—.-.257-—«=w 191,-—«s«1.000 «| 250 5 1/2 .806 ;600 11,430 61.1 79.0 FIG. 10 
> OM ou02 257.191 1.000 1.250 5 1/2 808 70,600 11,5° 61.3 FIG. 10 
120 0403 258 158 625 €25 5 1/2 747 71,500 6610 73.5 75.5 FIG. 11 
Hy 3120 Ou02  »=s-«.257 «. 163—Ss«w B25 —Ss« 25 5 1/2 740 76,860 6850 77.1 FIG. 11 
Me 3.747 0399 .256 163 .75@ .625 5 1/2 782 70,60C 8300 78.4 77.8 FIG. tI 
a 3.753 Ou00 )8=—:«257——«*EZ?-Ssis75O—s «B25 5 1/2 782 71,f00 8260 77.0 FIG. 11 
% 4.367 .OMOW .256 .162 .875 .€25 5 1/2 813 70,800 10,200 81.8 79.0 FIG. II 
28> © 4.377 OuO4 =. 25E «1E3)— 6875.25 5 1/2 .613 70,800 10,200 81.6 FIG. 11 
2 Biz \ 7025 72.2 FIG. 12 SHEARED RIVETS 
me Hee 6397 osu 130 628 00 374 B12 LS 7090 73.9 77.6 FIG. 12 SHEARED RIVETS 
4007 «© see, 183 180780 ooo 3 3u eur forloo «7D Os 65.5 FIG, 12 SHEARED RIVETS 
ile W023 0398 ‘25u 131.875 “B00 5 3/4 :864 70,100 7940 59.5 FIG. 12 SHEARED RIVETS 
3b 4.831 ouo2 253 130 875 500 5 3/4 866 70,100 7940 68.4 54.8 FIG. 12 SHEARED RIV 
% 5.523 0397 253 130 1.000 .£00 5 3/4 .882 70,100 7940 51.6 NOT PLOTTED SHEARED RIVETS 
% 5.500 0395 253 130 1.000 -50C 5 3/4 .882 70,10C 7790 51.1 48.8 NOT PLOTTED SHEARED RIVETS 
3 yy 78.4 79.4 FIG. 12 STEEL RIVETS IN FRONT ROW 
ae O308 oer st 608 £00 4 ay B10 70 108 Liver 9:7 FIG. 12 STEEL RIVETS IN FRONT ROW 
> wed oon $57 13! se £00 5 3/4 BuZ 70,100 8855 77.3 79.8 FIG. 12 STEEL RIVETS IN FRONT Row 
- * 2 g "25 é 2 70. FIG. 12 STEEL RIVETS IN FRONT ROW 
me Nise 6.050? a8 lst lens top § yu eo yorieo 9980 7s Fle. 12 STEEL RIVETS IN FRONT ROW 
Me OM. . . f : , 3 3 77.3 FIG. 12 STEEL RIVETS IN FRONT ROW 
3% 4.827 039! 257 .131 675 £00 5+ 3/4 864 70,100 10,355 La FIG. 12 STEEL RIVETS IN FRONT ROW 
% (5.523 .0392 257 1132 4.000 500 5 3/4 .881 70,100 10,850 71. . 12 S 18 peonT Roy 
% 5.503 0394 257.131 1.000 .500 5 3/4 881 70,100 10,655 70.0 71.2 FIG. 12 STEEL RIVETS IN FRO 
° 2. FIG. 13 
oe ee ee 8 oe ee ee: ee SS ee 
me 2.7N2 Se - ‘ 7340 73.6 74.0 FIG. 13 
% «3.017 0394 .191 .162 .562 ‘£00 5 3/4 .731 70,100 3. ho Sa 
3 70,100 7475 73.8 ; 
% 3.03 0398 191 .162 =. 562 £00 5 3/M ee 7210 73.9 76.5 FIG. 13 
% 3.uu2 .ONO! 191 163 .625 500 5 3/4 763 70,800 \ ‘ ; - 13 
%> = -3.439 o403 19 - 183 ‘-- £00 . Ay 7 a — oe 4. ‘3 SHEARED RIVETS 
x ° x ° ° e . 
oe 3787 Suge L191 “163 750 “£00 5 1/2 794 70,80C 8010 75.0 77.9 FIG. 13 SHEARED RIVETS 
TABLE Viti 
SERIES 28 
HINT WIDTH THICKNESS RIVET HOLE PITCH SPACING NO. OF m/p C, COUPON ULTIMATE OBSERVED CALCULATED CHART COMMENTS 
R DIAMETER 
VAMETE RIVETS IN Fu LOAD = EFFICIENCY EFFICIENCY REFERENCE 
oy dp REAR ROW 
la 3.745 0396 257.257 750 1.250 5 1/2 726 70,800 7510 7 2 
3 é ' ‘ ’ 1.5 72.7 FIG. 14 
. 3.750 .0395 -257 .257 +.750 ~=1.250 5 1/2 .726 70,800 7420 70.8 FIG. 14 
. 4,372 0397 -257 .257 .875 1.250 5 1/2 .765 70,800 9240 75.2 75.5 FIG. 14 
2 4.376 0298 -257  .257 .875 1.250 5 1/2 .765 70,800 9310 75.4 FIG. 1% 
& 4.996 Oude -257 .257 1.000 1.250 5 1/2 .794% 70,800 10,800 76.0 7£.5 FIG. 14 
% 4.995 0403 .257  .257 1.000 1.250 5 1/2 .794 70,800 10,825 76.1 FIG. 14 
4.50) -0402 -257 .257 1.125 14.250 4 1/2 .829 70,600 9840 76.6 76.5 FIG. 14 
%  Y¥.4gg 0400 .257 .257 1.125 1.250 4 1/2 .829 70,806 9690 77.6 FIG. 1% 
& 3.751 ouo2 -257 .257 750 .625 5 1/2 .726 70,800 7465 70.1 FIG. 14 
$3.74) 0396 .257 .257 .750 .625 5 1/2 .726 70,800 7460 71.2 72.7 FIG. 14 
& 4.376 0394 .257 .257 .875 .625 5 1/2 .765 70,800 9000 72.8 74.0 FIG. 14 
& 4.375 0397 .257 .257 .875 .625 5 1/2 .765 70,800 8850 72.1 FIG. 14 
hy 0 625 5 1/2 .79% 76,800 10,120 71.6 FIG. 14 DOG EARED TEAR 
-992 0400 .257 .257  «.00 62 2 
h oy 5 5 1/2 .794 70,800 10,150 71.9 74.2 FIG. 14 DOG EARED TEAR 
-998 .0399 -257 .257 1.000 «62 . 
& 4.50) 0402 -257 .257 1.125 .625 4 1/2 .829 70,800 9210 72.1 FIG. 14 DOG EARED Tea 
&® 4.500 040! -257 257 1.125 .625 4 1/2 .829 70,800 9120 71.4 74.8 FIG. 14 DOG EARED TEA 
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. 257 ‘ 2 5 ° 5 76.1 FIG. 15 
Ib 459 o408 257 257 -625 -625 5 3/4 .777' ~70,100 7375 74.5 FIG. 15 
2a 4, 130 -O408 267 .267 .750 625 5 3/4 .613 70,100 8950 75.6 76.8 FIG. 15 
2b 4.127 -o40e 257 «4.287 +750 -626 € 3/0 .813 70,100 9075 76.9 FIG. 15 
3a 4,820 -O405 267 .267 .€75 .62€ 5 3/4 .838 70,100 10,200 74.6 75.0 FIG. 15 
3b 4.810 -04C0 287 ct7 6756 625 5 3/4 .838 70,100 10,200 75.5 Fic. 15 
Ya §.620 O4C7 227 .267 1.000 625 5 3/4 .860 70,100 10,950 69.6 71.4 FIG. 15 BEARING FAILURE 
Ub §.£12 -O406 287 .287 1.000 -625 4 3/4 860 70,100 10,550 67.3 FIG. 15 BEARING FAILURE 
5a 3.745 -0403 -260 .260 - 750 - 750 € 1/2 .792 70,800 7875 73.8 76.8 FiG. 15 
5b 3.744 -O404 -260 .260 750 - 750 5 1/2 .792 70,800 7960 74.3 FIG, 15 
6a 4.372 -039¢ -260 .260 -875 -875 5 1/2 ..822 70,800 9410 77.3 76.8 FiG. 15 
6b 4,374 -0397 -260 .260 -875 -875 5 1/2 .822 70,800 9600 78.1 “FIG. 15 
7a 4,982 -O40! -260 .260 1.000 1.000 5 1/2 .844 70,800 10,890 76.9 75.4 FIG. 15 
7b 4.994 0399 257 .257 +4§.000 1.000 5 1/2 .845 70,800 10,850 77.0 FIG. 15 
fa 3.760 -0400 -$98 .191 -500 - 500 7 3/4 .796 70,100 8110 77.0 75.8 FIG. 15 
8b 3.756 o404 -191 191 - 500 - 500 7? 3/4 796 70,100 7900 74.3 FIG. 15 
9a 3.4uI -O404 -191 . 198 -625 - 500 5 3/4 .834 70,100 7400 75.9 76.5 FIG. 15 
So 3.458 o408 -191 . 191 -625 - 500 § 3/4 .834 =70,100 7390 74.7 FIG. 15 
10a 4, 130 0406 -19) = .1G1 - 750 - 750 £ c/4 .861 =70,100 8700 74.0 73.6 FIG. 15 
10b 4.112 0400 i9t 191 750 -750 £ 3/4 .861 70,100 8875 76.9 FiG. 15 
ila 3.115 0403 -257 =. 191 -625 - 7&0 1/2 816 70,800 6830 76.8 79.6 FIG. 16 
flo 3.119 0403 257 =. 191 625 -75 € 1/2 .816 70,800 6940 78.0 FIG. 16 
i2a 4.155 0408 257 —««. 191 -750 -750 5 3/4 .862 70,100 9350 79.1 80.1 FIG. 16 
12d 4.115 0398 -2&7 .193 -7£0 750 § 3/4 .869 70,100 8840 77.0 FIG. 16 
i3a 4.061 0402 257 198 813 750 17 1/2 .859 70,800 9180 79.2 80.0 FiG. 16 
13b 4.056 -O40! -257 +. 191 813 . 750 5 1/2 .859 70,800 9220 79.9 FIG. 16 
i4a 3.424 - 0400 -25€ .162 -625 ~7&0 5 1/2 .843 70,800 7120 80.5 81.0 FIG. 16 
uo 3.412 -0402 256 163 -62€ -750 € 1/2 .843 70,800 990 79.0 FiG. 16 
Asa 3.751 -O40! -257 .163 -75C -750 H ve ys Ls Go pore at 81.0 he 4 
ha) 3.736 o402 -256 .163 -750 750 1/2 .870 70,800 | » ° 
16a 4 133 0407 257.130 .750 .750 5 3/4 .906 70,100 9260 78.6 79.3 FIG. 16 STEEL RIVETS IN FRONT py 
160 = 1 C406 257.130 .750 .750 5 3/4 .906  7C, 100 8960 76.4 FIG. 16 ‘STEEL RIVETS IN: FRONT Rov 
17a 4. 136 . 0400 257.130 «=. 750. «575 E 3/4 .906 70,100 8550 73.7 FIG, 16 STEEL RIVETS IN FRONT Rov 
17> 6 0400 257 .130° «760 = «S78 £ 3/4 .906 70,100 84S0 73.1 FIG. 16 STEEL RIVETS IN FRONT Ry 
TABLE X 
SERIES 20 
JOINT WIDTH THICKNESS RIVET HOLE PITCH SPACING NO. OF m/p Ch COUPON ULTIMATE OBSERVED CALCULATED CHART 
NUMBER DIAMETER RIVETS IN Ftu LOAD EFFICIENCY EFFICIENCY REFERENCE 
dy dp REAR ROW 
la 5.990 -0397 -260 .260 1.000 1.000 6 1/2 .870 70,800 13,425 79.6 74.5 FIG, 17 
Ib 5.987 .0402 .260 .260 1.000 1.000 6 1/2 .870 70,800 1!3,335 78.3 FIG. 17 
2a 5.253 -0398 -256 .256 -875 -875 6 1/2 .853 70,800 11,850 79.9 78.6 FIG. 1? 
2b 5.251 -0399 -256 +;256 -875 -875 6 i/2 .853 70,800 11,870 79.9 FIG, 1? 
"Sa 4.501 -0405 -260 .260 .750 . 750 6 1/2 .827 70,800 10,290 79.6 77.9 FIG. 4 
3b 4.497 -0402 -260 .260 £.750 - 750 6 1/2 .827 70,800 9940 77.7 FIG. 
TABLE XI 
SERIES 3A 
JOINT WIDTH THICKNESS RIVET HOLE DIAMETER PITCH SPACING NO. OF m/p C, COUPON ULTIMATE OBSERVED CALCULATED CHART COMMENTS 
WUMBE! dy do a RIVETS Ih Ftu LOAD EFFICIENCY EFFICIENCY REFERENCE 
REAR ROW 
5 3/4 .628 71,000 6075 63.7 65.1 FIG. 18 
5 3/4 .628 71,000 6050 64.1 FiG. 18 
5 3/4 .680 71,000 8000 69.7 69.9 FIG. 18 
5 3/4 .690 71,000 7800 68.1 FiG. 18 
5 1/2 .726 71,000 9300 70.7 FIG. 18 
5 1/2 .727 71,000 9640 73.8 72.6 FIG. 18 
5 1/2 .772 71,000 12,150 78.5 75.7 FIG. 18 
5 1/2 .772 1,000 11,390 72.9 FIG. 18 
& 1/2 .804 71,000 11,500 79.3 77.6 FIG. 18 
y 1/2 .804 71,000 11,050 75.2 FIG, 18 
4 3/4 .849 70,800 13,680 75.1 FIG. 18 
4 3/4 .849 70,800 14,110 76.3 77.2 FiG. 18 
u 3/4 .669 70,800 16,890 74.6 75.5 FIG. 18 
o 3/4 .669 70,800 16,000 74M FIG. 18 
y 3/4 .886 70,800 16,350 66.8 65.3 FIG. 18 BEARING FAILURE 
u 3/4 .886 70,800 16,390 67.1 FiG. 18 BEARING FAILURE 
5 3/4 .767 70,800 8410 75.8 75.6 FIG. 18 
5 3/4 .768 70,800 8470 74.4 FIG. 18 
5 3/4 .823 70,800 12,110 gi.) FIG. 18 
5 3/4 .823 70,800 11,910 79.5 77.8 FIG. 18 
7) 3/4 .863 70,800 11,770 75.6 75.9 FIG. 18 
4 3/4 .863 70,800 11,680 75.0 FIG. 18 
yu 3/4 =.886 0.800 12,830 69.9 67.4 FIG. 18 BEARING FAILURE 
~ 3/4 .886 70,800 12,980 70.3 FIG. 18 BEARING FAILURE 
5 3/4 .676 71,000 3950 68.4 68.8 FIG. 18 
5 3/4 .676 71,000 3925 67.7 FIG. 18 
6 3/4 .760 70,800 6760 76.0 75.0 FIG. 18 
6 3/4 .760 70,800 6700 75.6 FIG. 18 
5 3/4 .811 70,800 7460 78.5 77.7 FIG. 18 
5 3/4 .811 70,800 7460 78.7 FIG. 18 
5 3/4 .605 71,800 26,930 61.1 63.0 FIG. 18 
5 3/4 .606 71,800 27,320 61.6 FIG. 18 
5 3/4 .654 71,800 33,850 67.2 67.2 FIG. 18 
5 3/4 .654 71,800 32,350 64.6 FIG. 16 
4 3/4 .729 71,800 36,950 71.8 72.8 FIG. 18 
y 3/4 .729 71,800 37,350 72.0 FIG. 18 
y 3/4 .774 «71,800 45,370 73.1 76.0 FIG. 18 
y 3/4 .774 71,800 46,270 744 FIG. 18 
y 1/2 .782 71,800 49,260 76.7 76.5 FIG. 18 
y 1/2 .782 71,800 48,650 75.6 FIG. 18 
3/4 .811 70,800 7390 78.2 78.0 FIG. 19 
H hy -810 70,800 7260 79.0 FIG. 19 
5 3/4 .865 70,800 10.5 3 80.3 4 ~ 
. 3/4 .865 70,800 10, . : RIVETS IN FRONT ROV 
° 70,800 12,060 79.6 77.5 FIG. 19 STEEL RI Rov 
Y —«-B/ 1907 70800 12.528 81.9 FIG. 19 STEEL RIVETS IN FRONT 
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SERIES 3B 

joint WIDTH THICKNESS RIVET HOLE DIAMETER PITCH SPACING NO. OF m/p C, COUPON ULTIMATE OBSERVED CALCULATED CHART 

WUMBER ' 2 3 RIVETS Ftu LOAD = EFFICIENCY EFFICIENCY REFERENCE COMMENTS 

REAR ROW 

la 4. 125 0395 257.257 «191 750 .750 5 3/4 .861 71,000 9110 78.7 FIG. 

ib M25 0396 257 + .257 ~. 191 750 .750 5 3/4 .861 71,000 9450 81.5 81.3 a 30 

2 «= 4.985 0390 257 +257. «.191 += 1.000 =~. 750 5 1/2 .885 71,000 10,60 77.1 81.3 FIG. 20 

2 4.99 0394 257 +.257. .191 ~=«§.000 =~. 750 5 1/2 .885 71,000 10,725 76.9 FIG. 20 

32 5-619 0394 257.257: «191 Ss1.125 Ss 875 5 1/2 .898 71,000 12,250 77.8 82.0 FIG. 20 

3 5.621 0390 256 .256 .191 1.125 .875 5 1/2 .898 71,000 11,860 76.1 FIG. 20 

Ya 6.247 0394 256 .256 .191 1.250 . .875 5 1/2 .90& 71,000 12,950 79.9 FIG. ZO 

yo 6.249 0396 256 .256 .191 1.250 .875 5 1/2 .908 71,000 14,100 86.2 81.0 FIG. 20 

5 6.867 0395 256 .256 .191 1.375 .875 5 1/2 .917 71,000 15,520 80.€ FIG. 20 

5 6.871 0390 256 .256 .192 1.375 .875 5 1/2 .917 71,000 15,310 80.4% 80.u FIG. 20 

a. %.1% 0399 257.257 ««. 163 750 .750 5 3/4 .882 71,000 9230 79.0 FIG. 20 

6 4.128 0394 256 .256 .163 750 .750 5 3/4 .882 71,000 r 81.7 82 0 FIG. 20 

ta 5.630 0396 257 + .257 162 «1.125 2—Ss«w 875 5 1/2 .91% 71,000 12,745 80.4 82.5 FIG. 20 

7 5.624 0393 257 +257 «.162)—Ci«=N«w125—Ss—«sw 875 § 1/2 .91% 71,000 12,540 80.0 FIG. 20 

g2 6. 244 0393 257.257 +163 «1.250 «=. 875 5 1/2 .922 71,000 13,950 80.0 FIG 20 

gb «6. 245 0396 256 .25€ .163 1.250 875 5 1/2 .922 71,000 14,250 81.2 62.0 FIG. 20 





yields results that agree well with the test data on rivet 
failures. In one respect, the graph disagrees with Hol- 
back’s results, which show a flat heel, rather than a 
sharp intersection, between the shear and crushing lines. 
However, Holback’s test points indicate that this heel 
is practically imperceptible for multiple-row joints with 
flat-head rivets. 


7. COMPARISON OF THE TEST RESULTS WITH THE 
ELEMENTARY THEORY OF JOINTS 


If it is assumed that there is no friction between the 
joined plates (Fig. 1), the load per in. taken out by rivets 
in the ith row is f¢(1 — C,), and this must be equal to 
the increment to,C; — to;_,C;_,; of the sheet load; i.e., 


oC; = o3-1C;,-1 + fi(1 — Cy) (2) 


This formula remains valid with 7 = 1 if C) is defined 
to be zero. Consequently, the solution of Eq. (2) is 


é 
oC; = Pa sea Cx) (3) 


These equations express the average sheet stresses on 
the various rows in terms of the rivet bearing stresses. 

It is assumed that in any set of consecutive rows, 
where yielding of the sheet is not occurring, the bearing 
stresses on all of the rivets are equal. This implies that 
the bearing stress on rivets of the front (mth) row is at 
least as high as the bearing stress on any other row. 
Consequently, for the usual type of joint, in which the 
tivets of the front row are no stronger than the rivets of 
any other row, a rivet failure should be expected when 
the bearing stress f, on rivets of the front row reaches 
the ultimate value f, that is given by Fig. 2. This 
hypothesis ignores the possibility of partial failure of a 
row of rivets; e.g., severe detrusion of the front row 
without complete shedding of the load of that row. 
The test data indicate that this phenomenon does not 
markedly affect the ultimate load of a joint that has 
rivets of a single material. 

The rivet rows are assumed to be so widely spaced 
that zigzag tears cannot occur. Then, in view of the 
preceding remarks, the most general type of failure 


TABLE XiJ 











consists of a tear on some row combined with a rivet 
failure on the nth row, or of a simple tear on the mth 
row. Hence, if stress concentration is neglected and if 
the sheet is assumed to yield indefinitely at the stress 
Fy, any mode of failure is represented by an equation 
of the type 


Cq = Fu, Sati = fags =... = fe =n 
in which @ is a number in the range 0, 1, 2, 3, ..., m. 
Then, fori = a, Eq. (3) becomes 
FyC, = Dfi( — Cy) (4) 


Also, since ¢,C, = o, Eq. (3) gives 
c= Shl-G)+fA-G) ©) 
k=l k=a+l1 


The first sum may be eliminated from Eqs. (4) and (5), 
and there results 


c= F.C. + fd eu. Cx) (6) 


A dimensionless factor ) is defined by \ = f,/Fy,. Then 
Eq. (6) may be expressed 


E, = Ca +r>X(1 — G) (7) 
k=a+l1 
The efficiency of the hypothetical joint is the smallest 
value that is obtained by successively setting a = 0, 
ey ae eB 
A more graphic form of Eq. (7) is obtained if the 
bearing areas of the various rows are represented by 
their ratios r; with the bearing area of the front (nth) 
row; i.€., 


1 — C, = r,(1 — G) (8) 
Note that, according to this definition, 7» = (1 — C,)— 


andr, = 1. Eqs. (7) and (8) yield 
E,=1- (ra kon AD) (1 ” Ch) (9) 
kea+l1 


If Eq. (9) is regarded as an equation of loci in the 
(C,, E,) plane, it represents two pencils of straight lines 
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Fic. 8. Experimental data for single-row lap joints. 
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Fic. 10. Experimental data for joints of Series 2A with d,/d, = 


80 .84 


—one with center at the point (1, 1), and the other with 
center at the point (1,0). Lines through the point (1, 1) 
are called ‘‘sheet-tear’’ lines because they represent 
failures that involve tearing of the sheet on some row of 
rivets. 

Lines through the point (1, 0) are called ‘‘rivet- 
failure’ lines because they represent failures that do 
not entail tearing of the sheet. These lines corre- 
spond to the value a = 0, and, therefore, their equa- 
tion is 


80 
76 
72 
E 
68 a 
d,-d, + '%4 
64 t + 0.040 
e= V2 
60 Mm/p = V4 
60 i7et SB" 12 16-30 BA 86 
C2 
Fic. 9. Experimental data for joints of Series 2A with d,/d, = 1. 
£ 
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Experimental data for joints of Series 2A with d,/d, = 
8/5. 
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Fic. 13. Experimental data for joints of Series 2A with d,/d: = 
6/5. 


4 d,= 6 d2= 5/32 
t = 0.040 ‘A 
e = 3/e 
s = Vz 


64 
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68 72 


E, = H(1 — CG), 


For any specified joint, the appropriate line of the 
pencil of rivet-failure lines is determined by the param- 


eter H (Eq. (10)). 


When the joint pattern is defined 
by preselected ratios r; (see Eq. (8)), there is a maximum 
value for H which is determined by the condition \ = 


Experimental data for joints of Series 2A with d,/d, = 2 
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Fy,/F mu (see Eq. (1)). The corresporiding rivet-failure 
line is called the ‘‘crushing” line, because crushing 
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lines of the pencil of rivet-failure lines are called ‘‘rivet- data are given in Table IV but most of the data on 
shear’’ lines. It should be observed that the rivet-shear perforated sheets have not been included in this report. 
line for a joint depends upon the sheet thickness, but The factor K is assumed to be a function of f,,/F,,. 
that the crushing line is theoretically independent of this Since the total bearing area in the joint is 
dimension. 
Theoretical lines of sheet-tear, rivet-shear, and wid( ~ Ca 
crushing failure, determined by Eqs. (9) and (10), are i=1 
plotted on the charts of the experimental data PP sete th e total load on the joint is 
8-20). Theoretical sheet-tear and crushing lines are : 
drawn with short dashes; rivet-shear lines are drawn Wte = WtEFy, = WtKRoC,F x 
with long dashes. These lines are based on the values of 
f, given by Fig. 2, and on the value F;, = 70,000 Ibs. the parameter f,/F,, is given by the equation 





per sq.in. Since back-row tears did not occur in the F.KRC 
tests, only the 45° sheet-tear line, E, = C,, (i.e., the fa/ For = — “ ord we 
line for which a = n), which represents sheet-tear on RT a= Co 


i=1 


the front row, has been plotted. 
It is seen that the theoretical lines are useful guides This muiy tut writaee f./¥,,. » Ka, in which 
for extrapolating the experimental curves. In general, 


the test data follow the theoretical rivet-shear and F RoC, 
crushing lines closely. The agreement of the data with emittte 5 ta (11) 
the sheet-tear lines is less favorable. For small values Fy, yi (1 — C;) 


of C,, the experimental efficiency exceeds the theoretical ; 
thus, the efficiency correction factor R exceeds unity Hence,the assumption that K is a function of the param- 
when C, is small. However, the ideal and the experi- eter f,/F,, is tantamount to the assumption that K 
mental curves usually cross in the interval 0.70 < C, < isa function of K¢. Explicitly, this means that K is a 
0.80. Furthermore, there is some interaction between function of ¢. 
the factor R and the bearing stress f. Consequently, The variable K has been plotted versus the abscissa 
there is not a sharp intersection between the line of tear by a process of curve fitting applied to the data of 
and the line of bearing failure. Rather, the intersection Tables VI to XII (Fig. 22). To substantiate Figs. 21 
formed by these two lines becomes well rounded. and 22, the efficiencies of all the test joints have been 
computed. Observed and computed efficiencies are 
tabulated side by side in Tables VI to XII. 
According to the preceding discussion, the method 
The theoretical efficiency of a joint is the smallest for computing the efficiency of a given 24S-T joint that 
value that is given by Eq. (7) when a@ takes the values is formed with protruding-head rivets is as follows: 


8. STRESS ANALYSIS OF JOINTS 


0,1,2,3,...,”. For a = 0, Eq. (7) becomes E = H X (1) Calculate the geometric rivet factors C; of all 
(1—C,). This is the equation of the rivet-failure lines rows. (Ratio of net area to gross area on a row.) 
shown on the experimental charts (Figs. 9-20). Since (2) Obtain the maximum bearing stress f, for rivets 


the test points follow these lines closely, it may be con- of the front row. Fig. 2 determines this stress for joints 
cluded that Eq. (7) is practically satisfactory for the that have 17S-T protruding-head rivets in the front 
case a = 0. Although the present tests did not include (mth) row and rivets at an edge distance 2d, in the first 
cases of back-row tears, it is likely that Eq. (7) remains row. For other cases, Eq. (1) may be used in lieu of an 
valid for these cases (a = 1, 2,3, ...,”—1), since such experimental chart, although allowance must be made 
failures are precipitated by failures of the leading for the variation of F,, with edge distance. In stress 
rows of rivets. However, it is evident that front-row analysis work, specified allowable stresses should be 
tears require special investigation, for here strain con- - used in conjunction with Eq. (1). (See reference 3 for a 


centrations exert significant effects. more complete discussion of these matters.) 
The efficiency of a joint that fails by a front-row (3) Calculate the joint parameter \ = f,,/ Fy. 
tear may be expressed in the form E = RC,. The fac- (4) Determine the factor RoC, from Fig. 21. If 


tor R depends partly upon the geometry of the rivet the front row is staggered with respect to the adjacent 
pattern and partly upon the bearing stress. Itisthere- row, an interpolation of this chart based upon the ratio 
fore convenient to split R into factors Ro and K, in S/p, is necessary. For patterns in which rivets are 
which Rp is the efficiency correction factor for the case omitted in the front row and S/p, > 0.50, the curve 
in which bearing stress is absent. for tandem patterns should be used, irrespective of 
In Fig. 21, the factor RoC, is plotted versus C, for whether the joint is staggered or tandem. 
tandem-riveted joints, and for staggered joints with (5) Calculate the factor ¢ by Eq. (11) and obtain K 
5/p ratios not less than 0.50. These curves have been from Fig. 22. 
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failures result when test points fall upon it. The other derived from tests of perforated sheets. Some of the 
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Fic. 14. Experimental data for joints of Series 2B with d,/d, = 1. 
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Fic. 15. Experimental data for joints of Series 2C with d/dz = 1. 
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Experimental data for joints of Series 2C with values 
of d,/d2 indicated above. 


Fic. 16. 


od,*d,-¥4 s-¥e mp- 
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C, 
Fic. 17. Experimental data for joints of Series 2D. 
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d,+*d2:d;-'4 $-5e mp-% t 
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Fic. 18. Experimental data for joints of Series 3A with d,/d; = 
z= 1. 
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Fic. 19. 


Experimental data for joints of Series 3A with d,/d; = 
2 and d2/ds = 3/2. ‘ 


a o d,=d,. V4 d;- 6 
a dd,» d= V4 d;3. S\6 
ad,+d,-% ds- ¥32 
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Fic. 20. Experimental data for joints of Series 3B. 
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Fic. 21. Efficiencies of perforated sheets. Curve A applies for: 
(1) all tandem patterns with p,—1 = pn; (2) staggered patterns 
with fr —1 = pnand S/p,n 2 0.80; (3) all patterns with p, = 
2p, _, and S/p, = 0.50. Curve B applies for staggered patterns 
with pn —1 = pnand S/ppn = 0.50. 
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(6) Determine the efficiency for front-row tear by 
the formula E = KR,C,. 

(7) Determine the efficiencies corresponding to 
other types of failure by Eq. (7), in which a takes the 
values 0, 1, 2, ..., m—1. The most important case is 
that in which failure results by shearing or crushing of 
all rivets (a = 0). For this case Eq. (7) reduces to 


B= (1- G) =H -G) 
=] 


(8) The actual joint efficiency is the smallest value 
that is obtained in steps (6) and (7). 


9. DESIGN OF EFFICIENT JOINTS 


The problem arises to determine the best proportions 
for a joint in which the ratios 7; (defined by Eq. (8)) 
are given constants. The optimum conditions for a 
joint of this type are 


fi=fp=... =fr=Sn 
Then Eq. (3) yields 


On = Fu 


Gy => (1 - CY 


i=l 
With Eq. (8), this gives 
C, = H(1 — C,) 
It follows that 
C, = H/(1 + H) (12) 


The theoretical efficiency is given by the equation 
E, = C,. This is the efficiency at the intersection of the 
sheet-tear line with the rivet-failure line (Figs. 9-20). 

Since Eq. (12) shows that C, is an increasing function 
of X, it is theoretically advantageous to choose \ as 
large as possible. This follows also from the charts 
(Figs. 9-20), since the steeper rivet-failure lines give 
higher peak efficiencies. However, the charts show 
that the abscissa of the maximum point on an experi- 
mental efficiency curve is somewhat smaller than the 
abscissa of the intersection of the theoretical sheet-tear 
and crushing lines. Actually, there is found to be no 
advantage in values of \ greater than 1.70 because of 
the strong interaction between bearing stress and rivet- 
hole strain concentration as the ultimate bearing stress 
is approached. Values of \ greater than 1.70 may be 
admitted without loss of efficiency, but the number 1.70 
is close to the upper limit that can be attained without 
exceeding Army-Navy allowable stresses. 

The procedure described below will generally furnish 
joints that lie on or near the maximum points of the 
experimental efficiency curves (Figs. 9-20). 

(1) Select the factors 7; of Eq. (8), noting that, for 
wide joints, 1 — C,; = d;/p;. Good joints can be ob- 
tained with all rows alike (r; = 1), and this condition 
should always be adhered to if the connected sheets 
have the same thickness. However, if a splice is made 
between unequal sheets, and if the utmost efficiency is 
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Fic. 22. Bearing correction factor for sheet efficiencies. 


required in order to show a positive stress-analysis 
margin or to permit a reduction of sheet thickness, the 
rivet factors C; should be varied from row to row. 
For maximum efficiency, the ratios of the factors 1 — C; 
in consecutive rows should be 4/3 or greater; i.e., 7, = 1, 
ta-1 = 4/3, tr-2 2 16/9, r,-3 2 64/27, etc. These 
conditions can be realized by varying the rivet sizes or 
the rivet pitches, or by a combination of both steps. 
However, the test data slightly favor tandem arrange- 
ments, in which the rivet diameters are varied. Values 
of (1 — C,-1)/(1 — C,) may be chosen arbitrarily above 
4/3 but it is practically inadvisable to make these 
ratios larger than necessary. Thus, the designer has a 
great deal of freedom in the selection of r-factors. 

(2) Select the rivet diameters and compute the 
joint parameter \ = fn/ Fu, being guided, as far as is 
feasible, by the condition that \ shall be at least as great 
as 1.70. Otherwise, efficiency will be sacrificed because 
of an early cutoff line (i.e., rivet-shear line). The stress 
F, is given by Fig. 2 for the case of 17S-T protruding- 
head rivets at 2d, edge distance on the first (rear) row. 
Eq. (1) approximately determines f, in all cases of pro- 
truding-head rivets, if the ultimate bearing stress F;, is 
properly determined (see also reference 3). If rivet 
diameters are varied from row to row, it may be desir- 
able to use small steel rivets in the front row in order 
to avoid getting too large rivets in the back rows. 

(3) Select the number of rows n. If the condition 
» 2 1.70 is satisfied, 80 per cent efficiency can be 
realized with a two-row joint and there is consequently 
no justification for using more rows. If \ is less than 
1.70, the number of identical rows of rivets required 
to obtain efficiency E is 


1.1E 


eee, BK OB X< 170 «(13 
in = (13) 


This equation is derived by introducing the empirical 
factor 1.1 into the formula obtained by solving Eq. (12) 
for n, with H = nd, and E = C,. (Observe that the 
relation H = nd applies when all rows are identical.) 
Eq. (13) is valid for values of E up to 80 per cent. 
Higher values of E are nearly always precluded by 
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strain concentrations. If Eq. (13) gives a noninteger 
value of n, the next whole number should be adopted. 
However, by varying the sizes or pitches of rivets in a 
joint it is occasionally possible to obtain the required 
efficiency with a smaller number of rows of rivets than 
that given by Eq. (13). 

(4) Calculate the rivet factor C, for the front row. 
This is accomplished by choosing the smaller of the two 
values given by the equations 


C, = H/(1 + A) 


and 


C, = 1.70 > 7r,/(1 + 1.70 Dor) 

i=l i= 
(See Eq. (12).) The factor C, determines the rivet 
pitch on the front row, and the pitches on the other 
rows are then determined by the general layout that is 
indicated by the preselected factors r; and the rivet 
diameters. 

(5) Determine the spacing between rivet rows. 
For tandem joints with all rivets of the same diameter 
and pitch, the spacing between the rows is irrelevant 
provided that the sheet does not shear through the rivet 
holes. Shear of the sheet through the rivet holes is 
ordinarily precluded if S/d > 2.5. 

For joints with unequal rivet diameters in the vari- 
ous rows, the row spacings should be fixed by the condi- 
tion 


S/p > 0.80 (14) 


For joints that omit alternate rivets in the front row, 
the condition, Eq. (14), may be modified to 


S/Pn > 0.50 (15) 


Eq. (14) should be observed for staggered joints with 
equal rivet pitches and diameters. When the foregoing 
criteria for row spacings are observed, it is irrelevant 
whether the patterns are staggered or tandem. Joints 
having different rivet pitches in the various rows are 
essentially staggered in any case, and for such joints 
the relative arrangements of the consecutive rows have 
no determinate effects. 

(6) Check the joint efficiency by the methods given 
im Sect. 8. 


10. DISCUSSION OF RESULTS 


(1) Joint Efficiencies —Joint efficiencies, calculated 
by the methods of Sect. 8, are shown side by side with 
the observed efficiencies in Tables VI to XII. The 
agreement between the observed and predicted effi- 
ciencies validates the proposed stress-analysis procedure 
for 24S-T joints with 17S-T rivets. 

Analytical stress-analysis methods for single-row 
joints have not been presented. Such methods are 
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unnecessary since Fig. 8 serves as a complete design 
chart for these joints when appropriate rivet-shear lines 
are introduced. It is observed that the maximum 
efficiency for a single-row joint is 60 per cent, and 
that this is obtained when the rivet factor C is 
0.60. 

(2) Design of Joints —The charts of Figs. 8-20 show 
that a curve of E versus C, approximately follows the 
line E = C, up to the peak efficiency. However, the 
curve bends down sharply and quickly joins a rivet- 
failure line when the abscissa for maximum efficiency 
is exceeded. The methods for joint design which are 
presented in Sect. 9 may be confirmed by the observa- 
tion that they invariably furnish joints that fall on or 
near the maximum points of Figs. 9-20. 


(3) Partition of Load.—The test data indicate that 
the distribution of load among the rivets of a tension 
joint at ultimate load is usually such that the bearing 
stresses on all rivets are equal. Differences in bearing 
stresses occur only if the sheet is yielding on some row 
of rivets. 

The frequent assumption that shearing stresses on all 
rivets are equal—which is true when all rivets have the 
same diameter—is incorrect when mixed rivet sizes 
occur in a joint. For example, Joint No. 32b, Table 
VII, failed by shearing the front-row rivets without 
other damage to the joint. By Fig. 2, the maximum 
bearing stress for this joint is 98,500 lbs. per sq.in.; 
i.e., this is the bearing stress to shear the smaller rivets. 
The total bearing area of this joint is 


5(d, + d2)t = 0.0756 sq.in. 


Hence, the calculated load to shear the front-row rivets 
is 


0.0756 X 98,500 = 7,410 lbs. 


Actually, the failure of the front-row rivets occurred at 
7,790 lbs. On the basis of equal shearing stresses in all 
rivets (which is often the assumption), this joint should 
have sustained 12,000 lbs. without rivet shear. 

The above example is only one in a large number of 
examples that verify the accuracy of the assumption 
that, in a set of consecutive rows where sheet yielding 
is not occurring, the rivets carry equal bearing stresses. 
This fact is illustrated by the close agreement between 
the test points and the theoretical rivet-failure lines on 
the charts (Figs. 9-20). 


(4) Staggered Joints Versus Tandem Joints —The 
data indicate that staggered joints with constant pitch 
will yield efficiencies equal to those of similar tandem 
joints if S/p is greater than 0.80. Strain concentrations 
cause deleterious effects in staggered joints with S/p 
less than 0.80, even though zigzag tears do not occur. 
Staggered joints in which some front-row rivets are 
omitted are equivalent to similar tandem joints if S/Ps 
is greater than 0.50. 





(5) 
data, 
are p 


(6) 
of rec 
to si 
slight 
ratio 
rows 
ever, 
alike. 

(7) 
stress 
the sh 
not gr 
eared 
prono 


Dear S 

I wor 
issue of 
body of 
ing cas 

I bel: 
the situ 


where 
There 


This is i 
For a 
compres 


This sg 


or simil; 








that 
sion 
ring 
ring 
row 


1 all 

the 
izes 
able 
out 
um 


ets. 


vets 


1 at 


uld 


r of 
ion 
ing 
es. 
een 


The 
tch 
em 
ons 
/p 
ur. 
are 


[Da 





24S-T RIVETED TENSION JOINTS 147 


(5) Spacing of Rivet Rows.—In the light of the test 
data, the following recommendations for row spacings 
are presented: 


(a) For tandem joints with constant p and d, the 
minimum spacing is 2.5d. 

(b) For other tandem joints, Eqs. (14) and (15) 
should be observed. 

(c) Staggered joiits with constant p and d should 
be designed according to Eq. (14). 

(d) Staggered joints in which some front-row 
rivets are omitted should be designed ac- 
cording to Eq. (15). 

(6) Effect of Variation of Rivet Sizes—The practice 
of reducing rivet sizes from large rivets in the rear row 
to small rivets in the front row generally results in 
slightly higher efficiency than variation of pitch. A 
ratio of rivet diameters of 4/3 in the consecutive 
rows is adequate to furnish peak efficiency. How- 
ever, good joints can be designed with all rows 
alike. 

(7) Effect of Margin Distance—When bearing 
stresses exceed the yield stress before tension failure of 
the sheet occurs, the highest efficiencies result with m/p 
not greater than 0.50. Excessive values of m cause dog- 
eared tears (Fig. 6). Such tears are accompanied by a 
pronounced drop in efficiency. 
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Letter to the Editor 


Dear Sir: 


I would like to discuss the comments given by W. R. Sears and Lester Lees in their Letters to the Editor in the December, 1946, 


issue of the JOURNAL. 


Dr. Sears agreed with my previous statement that to the first order the surface pressures for a very slender 


body of revolution are independent of the Mach Number, while Professor Lees disagreed and stated that this result was only the limit- 


ing case for a body of zero thickness. 


I believe that my following development, to be published in an early issue of the Quarterly of Applied Mathematics, will explain 


the situation more clearly. 





2x 
Hay —,* a oy\ Ill 
1 C\s = A ' 2x\ S 
en fig = [1+ met YE M~ tog2 ql (1-2) |r + 2(1- 2) 4 
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where S = rr? is the cross-sectional area of the body of revolution. 


I have shown that the surface pressures on a symmetrical body of revolution in subsonic flow are given by 


ae (r1)2 + O(r?) (1) 


Therefore, the first order approximation for any slender sharp nose body is given by Eq. (1) as 





C, = —(S"/m) log (r/L) + 0(r?) for M < 1 (2) 
This is independent of the Mach Number and is in agreement with Dr. Sears’ statement. 
For a symmetrical body the maximum pressure coefficient is obtained at x = L/2 where r! = 0 and at this point the second order 
compressibility correction is given by Eq. (1) as 
Cut) 7s log V1 — M? (3) 
Coo 1 + log (r/L) — [(L2/48) (S!Y/S¥) + ....] ; 


This shows that Professor Lees’ result (N.A.C.A. T.N. No. 1127) is a second order correction and is valid only for a prolate spheroid 


or similar body where S” = 0 for = 2. 


E. V. LAITONE 
Cornell Aeronautical Laboratory 








Room Temperature Tensile Properties of 
Aluminum-Alloy Sheet Following Brief 


Elevated Temperature Exposure 


J. T. LAPSLEY,* A. E. FLANIGAN,'* W. F. HARPER,? anv J. E. DORN™ 


University of California 


ABSTRACT 


The room temperature tensile properties of nine high-strength 
aluminum-alloy sheet materials have been determined after 
elevated temperature exposures covering a range of times and 
temperatures suitable for elevated temperature forming. Times 
ranged from 5 to 20 min. at temperatures up to 500°F. It is 
found that the naturally aged materials 24S-T, 24S-RT, and 
61S-W may be exposed (without deformation) for 20 min. at 
temperatures up to at least 500°F. with little or no loss in yield 
stress. Large increases may be expected at temperatures from 
450° to 500°F. The artificially aged materials 24S-T81, 24S- 
T84, 24S-T86, 61S-T, R301-T, and XB75S-T may be exposed as 
long as 20 min. (without deformation) at temperatures up to 
400°F. with little or no loss in yield stress but above 450°F. the 


loss becomes great. 

Additional tests made on material taken from deep-drawn 
cups formed at 300° and 450°F. indicate that the deformations 
involved in hot forming may induce an acceleration of the ele- 
vated temperature aging process. 


INTRODUCTION 


arora annt SHEET is often subjected to plastic 
forming. It has been found that some of the new 
high-strength materials exhibit poor forming character- 
istics at room temperature and that formability is im- 
proved at elevated temperatures. In general, the 
higher the temperature, the greater the improvement. 
As a result of such treatment, however, the precipita- 
tion hardenable materials undergo structural changes 
affecting their subsequent room temperature properties. 
The effects of elevated temperature aging have been 
studied extensively in connection with the development 
of the high-strength artificially aged materials! * 4 
but such investigations have not been concerned with 
the relatively short exposure periods used in hot forming 
and do not define limits of time and temperature suit- 
able for such operations. The investigation described 
in this paper was undertaken to evaluate the effects of 
short-time exposure and to serve as a guide in the selec-’ 
tion of permissible temperatures and times for hot 
forming operations. It is a portion of an extensive 
program directed by the War Metallurgy Committee 
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and sponsored by the Office of Production Research and 
Development, W.P.B., now the Production Research 
and Development Division, Department of Commerce, 
and is published with their permission.{f A related 
portion of the program deals with tensile proper- 
ties affecting formability at elevated tempera- 
tures.° 


REMARKS ON ELEVATED TEMPERATURE AGING 


As is well known, the high-strength aluminum alloys 
achieve their remarkable properties as a result of pre- 
cipitation hardening. Their production involves a 
solution heat-treatment followed by aging (precipita- 
tion hardening) either at room temperature or at suit- 
able elevated temperatures. A considerable increase 
in tensile and compressive yield stresses results from 
precipitation at room temperature (natural aging). 
The elevated temperature treatment (artificial aging) 
generally results in a much greater increase, the gain in 
yield stress being attended by a smaller increase in ten- 
sile stress and by a decrease in elongation. 

In the first stages of elevated temperature aging the 
yield stress may decrease slightly. Subsequently it 
increases steadily until it reaches a maximum value 
characteristic of the aging temperature employed. In 
this condition the material is said to be ‘‘fully aged.” 
Further elevated temperature exposure results in 
“overaging”’ attended by a decreasing yield stress. In 
the production of artificially aged aluminum alloys, ag- 
ing times and temperatures are selected to achieve full 
aging. The following important characteristics are 
associated with the process: 

(1) The higher the temperature, the more rapidly 
is the fully aged condition attained. 

(2) The higher the temperature, the lower the 
maximum strength attainable in general. 

(3) In the case of the alloy 24S, the yield stress 
attainable as a result of artificial aging may be in- 
creased by the application of cold work preceding the 
aging treatment. For this treatment, the greater the 
degree of deformation, the greater the yield stress 
attainable. 


tt The work described in this paper was carried out by Project 
NRC-548 at the University of California. 
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CLASSIFICATION OF MATERIALS INVESTIGATED 


The high-strength materials may be grouped as 
follows according to the anticipated effects of elevated 
temperature exposure upon subsequent room tempera- 
ture properties. 

Group A.—Precipitation hardenable alloys previously 
aged at room temperature. These materials are sub- 
ject to additional precipitation at elevated tempera- 
tures. Brief exposure at temperatures below approxi- 
mately 500°F. may be expected to increase the room 
temperature yield stress. With longer exposures and 
higher temperatures overaging may be anticipated. 

Group B.—Precipitation hardenable alloys previously 
aged at elevated temperatures. Additional exposure 
may be expected to induce little, if any, strengthening. 
Overaging may be expected at higher temperatures. 

The materials studied in the present investigation 
were : 








Group B 
Group A “ 
24S-T 24S-T81 61S-T 
24S-RT 24S-T84 R301-T 
61S-W 24S-T86 XB75S-T* 





* The former experimental material, XB75S-T, is now desig- 
nated 75S-T by the producer, the Aluminum Company of 
America. 


All materials were obtained from the producer’s 
regular stock. With the exception of R301-T, the 
material was bare (not clad). Nominal chemical com- 
positions and standard heat-treatments are listed in 
Tables 1 and 2. The 61S alloy exhibits strengths in- 
ferior to the remaining alloys and perhaps should not be 
classed as a high-strength material. It is included here, 
however, because it is subject to precipitation harden- 


ing. 
SCOPE OF THE INVESTIGATION 


This paper is concerned principally with tensile tests 
of specimens subjected previously to elevated tempera- 
ture exposure for short intervals without attendant def- 
ormation. Since simultaneous deformation (as in hot 
forming operations) may modify precipitation effects 
and thus cause changes in room temperature properties, 
it might be considered desirable to employ specimens 
subjected to controlled deformation at elevated tem- 
perature. The scope of such an investigation was con- 
sidered too broad for the present investigation and, 
also, the results obtained would apply only for the types 
and degrees of deformation employed. For these rea- 
sons deformation was not included in the preparation of 
the specimens. 

In addition to the tensile tests on undeformed speci- 
mens, a limited number of tests were made on material 
taken from the wall section of cups which had been 
deep drawn at elevated temperatures. These speci- 





TABLE 1 
Nominal Composition of Alloys Tested 


248 —4.5% Cu, 0.6% Mn, 1.5% Mg, balance Al 

61S —0.25% Cu, 0.6% Si, 1.0% Mg, 0.25% Cr, balance Al 

R301 —4.5% Cu, 1.0% Si, 0.8% Mn, 0.4% Mg, balance Al 

XB75S—The principal alloying elements are zinc and magnesium 
with smaller additions of copper and other elements. 
Further information is not available 














TABLE 2 
Standard Heat-Treatments Used in Production of Materials 
Tested 
24S-T 24S alloy solution heat-treated at 910-930°F. fol- 
lowed by quenching in cold water and aging at 
room temperature 
24S-RT 24S-T strain hardened by rolling to a reduction of 
about 5'/2 per cent 
61S-W 61S alloy solution heat-treated at 960-980°F. fol- 
lowed by quenching in cold water and aging at 
room temperature 
24S-T81  24S-T stretched 1 per cent and aged 10 hours at 375° 
F. or 12 hours at 365°F. 
24S-T84  24S-T stretched 3'/,-4 per cent and aged 6 hours at 
375°F. or 9 hours at 365°F. 
24S-T86 24S-RT aged 5!/2 hours at 375°F. or 8'/: hours at 
365°F. 
61S-T 61S-W aged 8 hours at 350°F. or 18 hours at 320°F. 
R301-T R301 alloy solution heat-treated at 930-950°F. 


followed by quenching in cold water and aging 6 
hours at 350°F. or 18 hours at 320°F. 

XB75S-T XB75S alloy solution heat-treated at 860-930°F. 
followed by quenching in cold water and aging 24 
hours at 250°F. 





mens, tubular in shape, were tested in compression since 
it was difficult to obtain tensile specimens from the 
cups. The results, which include the effects of de- 
formation during elevated temperature exposure, are 
compared with those obtained in the tensile tests. 


TEST PROCEDURE 


Tensile test specimens were prepared from 0.064-in. 
sheet and were stressed transverse to the direction of 
rolling. The reduced section was 1/2 in. wide, the 
specimen design conforming to A.S.T.M. Standard 
E8-44. All specimens of each material were obtained 
from a single sheet. 

Previous to testing, the specimens were exposed to 
temperatures of 300°, 400°, 450°, and 500°F. with 
exposure periods including 5, 10, and 20 min. Expo- 
sure was achieved by immersion in heated baths of 
Houghton No. 80 tempering oil. Steady-state tem- 
poral and spatial temperature variations within the 
baths were no more than +1°F. Specimens were 
treated in small groups to minimize the mass effect and 
thus permit rapid attainment of temperature. To 
determine the transient heating conditions, a survey 
was made with a check specimen containing thermo- 
couples peened into the gage section. Two minutes 
after insertion in the 500°F. bath this specimen was 
found to be within 3°F. of the desired temperature. 
After 3 min. the differential was reduced to 1°F. 

Tensile stress, tensile yield stress (0.2 per cent offset), 
and elongation in 2 in. were determined in the tensile 
tests. An averaging extensometer having a 2-in. gage 
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length and a multiplying factor of 2 was employed in the 
determination of yield stress. It was equipped with a 
“Last Word”’ dial gage, the smallest graduation corre- 
sponding to a strain of 0.00005 in. per in. After frac- 
ture elongation measurements were made to the nearest 
0.01 in. with calipers and scale. The tests were per- 
formed on a screw-driven testing machine. The cross- 
head speed was approximately 0.03 in. per min. up to 
attainment of the yield stress. Subsequently it was 
increased to 0.25 in. per min. and maintained at that 
rate until fracture. 

The test results are presented graphically in Figs. 1 
through 9, each point on the curves representing an 
average value based on two or more tests. 


DISCUSSION OF RESULTS 


Of the properties evaluated, the tensile yield stress is 
probably of greatest interest. Accordingly it will be 
used as the basis of the following discussion. 

Considering first the alloys of Group A, it may be 
seen that 24S-T and 24S-RT (Figs. 1 and 2) exhibit 
similar behaviors. In each case 20-min. exposure at 
temperatures up to 300°F. has no appreciable effect on 
the final properties, an indication that precipitation at 
these temperatures is sluggish. Temperatures from 
300° to 425°F. cause small decreases in the room tem- 
perature yield stress, an effect consistent with similar 
decreases known to occur in the preliminary stages of 
the artificial aging of 24S.' The use of temperatures 
above 450°F. results in large increases in yield stress. 
As might be expected, the maximum values attained by 
24S-T and 24S-RT (61,000 and 66,000 Ibs. per sq.in., 
respectively) are only slightly lower than the original 
values for 24S-TS81 and 24S-T86 (64,000 and 68,000 
lbs. per sq.in.). The small differences may be ex- 
plained on the basis of time and temperature. With 
exposure periods of from 5 to 20 min. the maximum 
values are reached at temperatures from 450° to 500°F., 
whereas in the production of 24S-T81 and 24S-T86 a 
temperature of 375°F. is employed with a much longer 
aging period (refer to Table 2). As noted previously, 
the maximum attainable yield stress decreases with 
increasing aging temperature. 

Neglecting the small decreases occurring at temper- 
atures from 300° to 425°F., Figs. 1 and 2 indicate that 
24S-T and 24S-RT may be held briefly at temperatures 
up to at least 500°F. with no loss in yield stress and that 
large increases may be expected when the temperature 
is greater than 450°F. In this connection, however, 
consideration should be given to corrosion resistance, a 
matter tHat will be discussed presently. 

The remaining member of Group A, 61S-W (Fig. 3), 
behaves in a manner similar to that characterizing 
24S-T and 24S-RT. The maximum yield stress 
attained (36,000 Ibs. per sq.in.) approaches the initial 
value for 61S-T (38,000 Ibs. per sq.in.). In this in- 
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stance, however, there is no evidence of decreasing yield 
stress at the lower temperatures. 

Since in their production the materials of Group B 
are aged to develop maximum yield stress, subsequent 
exposure to elevated temperature may be expected to 
result in overaging with an attendant decrease in 
strength. Such a trend is shown in Figs. 4 through 9. 
For all of the materials the loss in strength is negligible 
for temperatures up to about 400°F. Above 450°F. 
the loss in strength becomes large. In the case of 75S- 
T, the original room temperature yield stress is 68,000 
Ibs. per sq.in. After 20-min. exposure at 400°F. it is 
65,000 Ibs. per sq.in., and with the same exposure at 
500°F. it becomes reduced to 33,000 Ibs. per sq.in., a loss 
of more than 50 per cent. Of the materials considered, 
75S-T and R301-T suffer the most rapid loss of strength. 

In the case of R301-T (Fig. 9), an increase over the 
original yield stress is noted after brief exposure at 
300°F. This probably indicates that in its production 
the original material had not been completely aged. 


Corrosion Resistance 


It is known that in some alloys the first stages of 
artificial aging are attended by a pronounced increase 
in susceptibility to corrosion. With continued aging, 
however, corrosion resistance improves steadily until in 
the fully aged condition it approaches that of the 
original material. Overaging of moderate duration 
does not appear to be harmful. 

It may be assumed, therefore, that the Group B 
materials may be formed at elevated temperatures with 
little or no loss in resistance to corrosion. Of the Group 
A materials tested, the initial loss in corrosion resistance 
is known to be particularly severe in the case of 24S-T 
and 24S-RT.! It would seem advisable, therefore, to 
allow these materials to attain a fully aged condition if 
they are formed at elevated temperatures. This might 
be accomplished in either of two ways: 

(1) By suitable choice of temperature and exposure 
time for the forming process. 

(2) By subjecting the formed material to a subse- 
quent aging treatment. 


Effect of Deformation at Elevated Temperature 


For the more usual case of hot forming in which both 
deformation and exposure to elevated temperature are 
encountered, the tensile test results may require inter- 
pretation and possible modification since they are in- 
fluenced by elevated temperature exposure alone. It is 
known that deformation of 24S alloy prior to elevated 
temperature aging increases the strength attainable 
and also accelerates the precipitation process. It is 
possible, therefore, that deformation at elevated 
temperature may also influence the subsequent room 
temperature properties of similar materials. A pre- 
liminary investigation of this effect was made by deter- 
mining the compressive yield stress for materials previ- 
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TABLE 3 
Room Temperature Compressive Yield Stress Values 
Deformed Material from the Wall Section of Hot-Drawn Cups 
Time at temperature prior to drawing—approximately 5 min. 
(Each value represents an average of two or more tests unless otherwise noted) 








300°F. Temperature 





-—Compressive— -—— Tensile —— 
Yield Stress, Yield Stress, * 
Material Lbs. per Sq.In. Lbs. per Sq.In. 
24S-T 67,700** 45,700 
24S-T86 65,200 68,100 
R301-T 59,600** 58,000 
XB75S-T 64,000** 67,900 


450°F. Temperature 





——Compressive— 
Yield Stress, 


-——— Tensile --—~ 
Yield Stress, * 


Material Lbs. per Sq.In. Lbs. per Sq.In. 
24S-T 67,500 48,500 
61S-T 43,000 38,400 
24S-T86 59,000 66,000 
R301-T 51,000 53,700 
XB75S-T 59,000 64,600 





* Tensile values for 5-min. exposure without deformation, taken from Figs. 1 through 9. 


** Single test value. 


ously deformed at elevated temperature. The program 
followed was limited in scope and served merely to sur- 
vey the problem for possible future investigation. 

Material for test was taken from the wall section of 
drawn cups formed at temperatures of 300° and 450°F. 
The alloys considered were those of the tensile test 
program. The drawn cups, 3 in. in diameter, were 
formed in a single draw from circular blanks of 0.064-in. 
sheet. Prior to drawing, the blanks were preheated in 
the draw press with a total exposure interval at tem- 
perature of 5 min. The draw conditions were such as 
to produce maximum height of cup without ironing. 
Test specimens in the form of rings 1 in. high and 3 in. 
in diameter were cut from the cup wall section and 
tested in compression to determine the yield stress. 
The results of these tests are presented in Table 3. For 
convenience of comparison, tensile test results for the 
same temperature and exposure conditions are shown 
also for each alloy considered. 

A direct comparison of the tensile test results with 
those of the hot formed material is complicated by 
several factors: 

(1) Normal differences known to exist between the 
tensile and compressive properties of directionally 
worked material. 

(2) Incomplete knowledge of the history and degree 
of deformation in the hot formed material. 

(3) Specimen materials from different sheets of the 
alloys. 

It is felt, however, that such a comparison has some 
merit as an indication of the effect of deformation at 
elevated temperatures upon resulting room temperature 
properties and also as an aid in applying the tensile test 
results. 

In the case of 24S-T formed at 450°F., Table 3 indi- 
cates a compressive yield stress of 67,500 Ibs. per sq.in. 
which is much greater than the corresponding tensile 
value. According to the tensile test results (Fig. 1), the 
latter could be expected eventually to reach a maximum 
value of about 61,000 Ibs. per sq.in. Perhaps the ten- 
sile value attainable in 24S-RT at this temperature 
should also be considered in view of the deformation 
applied in its production. The maximum yield stress 


in this case—about 66,000 Ibs. per sq.in. (Fig. 2), is in 
better agreement with the compressive value for hot 
formed 248-T. 

In the case of the same material formed at 300°F., 
Table 3 shows a compressive value of 67,700 Ibs. per 
sq.in. and again the tensile value is much lower. If the 
preceding explanation is applied, it is clear that the 
accelerating effect of hot forming must be rather great 
since the change at 300°F. is extremely sluggish in the 
absence of deformation. 

The materials of Group B are seen in Table 3 to ex- 
hibit values of compressive yield stress generally lower 
than the corresponding tensile values. Excluding the 
two exceptions (R301-T at 300°F. and 61S-T at 450° 
F.), the differences are rather uniform, averaging about 
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alloy. [—As-received. ©C—5-min. exposure. @—10-min. ex- 
posure. A—-20-min. exposure. 
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6 per cent. An acceleration effect accounts for the 
generally lower values since the anticipated changes are 
in the direction of overaging. The two exceptions, of 
course, cannot be explained on this basis. 

It may be noted that the yield stress obtained by 
forming 24S-T at either 300° or 450°F. approxiinates 
the original tensile yield stress of 24S-T86. Thus the 
value obtained by hot forming is nearly the same as 
that obtained by cold prestraining and artificial aging. 
This suggests the possibility that hot rolling or some 
other form of elevated temperature deformation may be 
useful in the production of high-strength sheet. 


APPLICATION OF TEST DATA 


In many forming operations the blank material is not 
deformed uniformly and some material may receive no 
deformation. In applications where the properties of 
such nondeformed material are critical the results here 
obtained are directly applicable for determining limiting 
temperatures and exposure times for hot forming opera- 
tions. 

The results of the compression tests indicate that hot 
forming results in more rapid aging and more drastic 
changes in room temperature properties than occur in 
the absence of deformation. Fortunately the differ- 
ences are not great in the case of the B materials where 
hot forming may be expected to be most useful. As 
noted previously, the yield stress after hot forming is on 
the average about 6 per cent lower than the tensile 
values. Differences are greater in the case of the A 
materials but the tensile values are on the conservative 
side, being lower than those resulting from hot forming. 
It is thought, therefore, that the tensile test results are 
of value in indicating permissible conditions for hot 
forming. 


SUMMARY AND CONCLUSIONS 


The room temperature tensile properties of nine high- 
strength aluminum-alloy sheet materials’ were deter- 
mined after elevated temperature exposures covering a 
range of times and temperatures suitable for elevated 
temperature forming. Times ranged from 5 to 20 min. 
at temperatures up to 500°F. The following conclu- 
sions may be drawn from these tests in which specimens 
were not deformed during the exposure period. 

(1) The naturally aged materials—24S-T, 24S-RT, 
and 61S-W may be exposed (without deformation) for 
times as long as 20 min. at temperatures up to at least 
500°F. with little or no loss in yield stress. Possible 
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increases in susceptibility to corrison may require con- 
sideration, however. Large increases in yield stress 
may be expected when the exposure temperature lies 
between 450° and 500°F. 

(2) The artificially aged materials—24S-TS81, 24S- 
T84, 24S-T86, 61S-T, R301-T, and XB75S-T may be 
exposed without deformation for as long as 20 min. at 
temperatures up to 400°F. with little or no loss in yield 
stress. Above 450°F., the loss increases rapidly. 

A small number of additional tests were made to 
determine the compressive yield stress of material 
taken from deep-drawn cups formed at 300° and 450°F. 
and these values compared with the tensile yield stresses 
for the undeformed materials. Although the com- 
pressive tests were few in number and the comparison 
suffers from certain defects, the following tentative con- 
clusions may be drawn: 

(1) The deformation involved in hot forming may 
induce a considerable acceleration of the elevated 
temperature aging process. 

(2) Tests on 24S-T indicate that in the case of 
naturally aged materials, the yield stress values result- 
ing from the forming at temperatures up to at least 
450°F. are higher than those indicated in tests involving 
no deformation during the exposure period. 

(3) In the case of the artificially aged materials, the 
yield stress values resulting from hot forming may be 
lower (about 6 per cent) than values indicated in tests 
involving no deformation. 

(4) The compressive yield stress developed by hot 
forming 24S-T sheet at either 300° or 450°F. (with total 
exposure at temperature of about 5 min.) equals the 
tensile yield stress developed in 24S-T86 by cold pre- 
straining and artificial aging. This suggests that some 
form of elevated temperature deformation may be use- 
ful in the production of high-strength sheet. 
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Some Ballistic Contributions to Aerodynamics 


A. C. CHARTERS* 
Aberdeen Proving Ground 


ABSTRACT 


A systematic study of projectile aerodynamics has been made 
in the spark photography range of the Ballistic: Research Labora- 
tory, and the results on drag are now sufficiently complete to 
have general application to compressible flow problems. The 
performance of conical heads is analyzed by the first order theory 
and by the theory of Taylor and Maccoll, extended to additional 
cone angles. An empirical formula representing the variation 
of drag coefficient with Mach Number is presented. The quan- 
tity J 1 + KpM? is shown to vary linearly with M at supersonic 
speeds for a wide vatiety of projectile shapes. Results on the 
length of boattailing are given. For a fixed boattail angle, it 
is shown that the drag coefficient decreases linearly with base area, 
that the side and base pressures are constant and independent of 
boattail length, and that their difference is portional to the slope 
of the linear function. The effects of banding, scale (size), inter- 
action of head on base, and skin friction are discussed briefly. 

Application is then made of these ballistic results to the design 


_of a supersonic aircraft fuselage shape giving minimum drag, 


with a fixed body diameter and with the arbitrary choice of a 
conical head and conical boattail. A 50-in. diameter fuselage 
at 50,000 ft. and a 20-mm. model at sea level are considered. It 
is found that the design and performance of the two are similar, 
with the optimum head angle (semivertex) being near 4°. Ata 
Mach Number of 1.5, base pressure and skin friction are nearly 
equal and account for over 85 per cent of the drag. As the Mach 
Number is increased to 5.0, the skin friction alone accounts for 75 
per cent of the drag. The boattail is found to effect a saving of 
23 per cent or more in drag, even at a Mach Number of 5.0. 
The power required for the 50-in. fuselage flying at a Mach Num- 
ber of 1.5 is computed to 2.263 hp. at 50,000 ft., which increases 
to the impractical figure of 21,243 hp. at sea level. 








SYMBOLS 
a,b = empirical constants in Q function (Q = a + bM) 
A = total wetted area (of projectile or flat plate, as the 
case may be—does not include base area) 
B = diameter of base at rear of projectile 
B = boattail angle (semivertex angle of cone generating 
boattail) 
c = velocity of sound in undisturbed atmosphere (or in 
free stream) 
CG = total skin friction coefficient for a flat plate in a sub- 
sonic flow 
C total friction force on flat plate 
du 1/spU*A 
G = average effective skin friction coefficient for entire 
projectile 
= total friction force on projectile 
f 1/opU2A 
D = diameter of body of projectile (also known as the 
“caliber’’) 
Y = ratio of specific heats (y = 1.40 for air) 





Presented at the Aerodynamics Session, Fourteenth Annual 
Meeting, I.A.S., New York, January 29-31, 1946. Received 
October, 1946. 
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Kp = drag coefficient (usually designates drag coefficient of 
complete projectile), Kp = drag force/pU?D? 

Kpo = drag coefficient at zero yaw; unless specified: Kp = 
Kpo . 

Kpu = head drag coefficient (pressure drag only), Koy = 
(w/4)/(pu/pU?) 

Kpsr = skin friction drag coefficient, Kpsr = (C,/2)(A/D?) 

Kps_ = base drag coefficient (includes both base and boattail 
drag), Kop = (x/4)[(ps + ps)/pU?) 

AKp = reduction in Kpg due to boattailing (includes skin 


friction and pressure forces) 

AKopsr = increase in Kpg due to skin friction acting on 
boattail 

AKpp = reduction in Kpz due to pressure forces only resulting 
from boattailing 


l,m = empirical constants in boattail equation [Kp = / + 
m(B/D)?) 

M = Mach Number (M = U/c) 

Nr = Reynolds Number (Nr = R) 

0,0, = semivertex angle of cone 

fi = atmosphere pressure or free stream pressure, as the 
case may be (absolute value—note that pU? = 
ypiM?) 

Du = overpressure acting on head (py = Pu — pr) 

bs = under pressure acting on side of boattail (ps = Ps — 
pi) 

bp = under pressure acting on base at rear of projectile 
(ps = Pa — fr) 

Pu = absolute head pressure 

Ps = absolute boattail side pressure 

Pp = absolute base pressure 

(—) = quantity averaged over cross-sectional area in- 
volved 

Ap = pu for a conical head 

Q =Vi+ KpM 

R = Reynolds Number 





Z U X characteristic length 
R = Nr= 
v 
p = air density in undisturbed atmosphere or free stream, 


as the case may be 
cross-sectional area of head shape (at any point. along 


7) 
ll 





head) 
U = velocity of projectile (or of free stream) 
v = coefficient of kinematic viscosity 
( ordinary coefficient of visonty) 
v= ——— 
p 

s = distance along the axis of the head from the point 

INTRODUCTION 


arog HAVE STUDIED the aerodynamics of 
projectiles at supersonic speeds for 50 years or 
more. They have made some progress in solving the 
problems of compressible flow involved, but for the most 
part their approach has been from an engineering stand- 
point and their methods have been largely empirical.’ 
Their success has been limited by the lack of any general 
theory to guide their research. 
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Fic. 1. 


The picture of slow progress in ballistic aerodynamics 
is now rapidly changing. Developments in theory have 
progressed to a point of limited practical application. 
Experimental techniques have improved radically and, 
in fact, to such a degree that the experimental data can 
serve as a sound basis for empirical generalization. The 

‘advance in theory, on the one hand, and in experiment, 
on the other, is already sufficient to have produced re- 
sults of general interest to the broad field of supersonic 
aerodynamics. 

One of the modern instruments for experimental 
research in ballistic aerodynamics is the spark photog- 
raphy range at the Ballistic Research Laboratory.? A 
view of the first group of stations is shown in Fig. 1 
looking up-range toward the firing point. This range is 
being used for a carefully planned program of research 
to study the aerodynamics of projectiles in free flight. 
The projectile is considered to consist of the various 
components making up its overall shape, such as the 
head, body, band, boattail, and base. The effect of 
each component in turn on the complete projectile is 
studied by making a series of systematic variations in 
the shape, changing only one component at a time. The 
performance of the whole projectile can be synthe- 
sized out of the differential performance of its 
parts. 

The basic program is only in its initial stages, but 
enough research has been completed to permit the 
formulation of empirical laws for (1) the variation of 
drag coefficient with Mach Number (M > 1), (2) the 
effect of boattail length on the drag, and (3) the effect 
of banding on drag. These laws, combined with theo- 
retical calculations for head pressure and skin friction, 
are adequate to compute the drag of certain simple 
types of projectiles at supersonic velocities. It is the 
purpose of this paper to present these ballistic contribu- 
tions, together with the related theory as needed, and 
to apply the two to the design of an aircraft fuselage 
giving the minimum supersonic drag. 
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HEAD SHAPE 


The head shape has been dealt with more effectively 
by aerodynamic theory than any other of the projec- 
tile’s components. A general first order theory has been 
worked out by von Karman and Moore for long, slender 
head shapes.* The special case of an infinite cone in a 
supersonic flow has been solved completely by Taylor 
and Maccoll for a nonviscous, irrotational gas, com- 
pressed adiabatically. Using the first order theory, 
von Karman has worked out the optimum head shape 
for minimum drag with a fixed head length.* 

The developments in the theory to date leave a great 
variety of head shapes untouched. However, the exact 
solution of an infinite cone is a skeleton on which the 
characteristics of many of the unsolved shapes can be 
based. The head shape can be defined in terms of one 
parameter, the head length in calibers, and one func- 
tion, the nondimensional head shape (given as a fraction 
of the caliber head length). The family of functions de- 
scribing the head shape can usually be constructed so 
that the cone is one member of the family. All the re- 
maining shapes are then variants of the cone. 

Taylor and Maccoll calculated the properties of the 


flow over cones at a series of supersonic speeds for cones ° 


having semivertex angles of 10°, 20°, and 30°. Similar 
calculations were carried out at the Ballistic Research 
Laboratory for an additional series of cones having semi- 
vertex angles of 7.5°, 9.0°, 9.5°, 12.1°, 15°, and 20° (re- 
peated). The results of this method give a discreet 
value of the drag coefficient at a particular Mach 
Number for the cone in question. The values were first 
plotted on a graph with Kp as ordinate and M as ab- 
scissa, and smooth curves were faired through each set 
of values for each cone angle. The results for cone 
angles of 4.5°,* 10°, and 15° are shown in Fig. 2. The 
values on this graph were then cross-plotted to give a 
graph of Kp as a function of semivertex angle, with 
Mach Number as the independent parameter. This is 
the graph desired by the projectile designer. The 








* The curve for 4.5° was cross-plotted from Fig. 3. 
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curves for Mach Numbers of 1.5 and 5.0 are shown in 
Fig. 3. 

The first order theory gives results that agree satis- 
factorily with experiment and with the exact theory for 
the cone. The drag coefficient for a cone derived from 
the first order theory is given as an explicit function of 
the Mach Number and semivertex angle by the for- 
mula‘ 


where a = VM? — 1. 

The comparison with the exact solution is shown by 
the dotted curves in Figs. 2 and 3. 

The optimum head shape giving minimum drag for a 
fixed head length is probably different from a cone. von 
Karman derives the formula from the first order theory® 
in differential form as follows: 


dS/dx = const. vi — (2x-/1)? 


where the head length is / and the point is atx = —//2.* 

His result appears surprising at first since it gives a 
shape that has a blunt point. This may not be far 
from the truth, however. Some firing tests have been 
carried out with a series of ogival heads and with heads 
with different meplats,f and some elementary calcula- 
tions have been made to determine the optimum meplat 
diameter for a conical head. The firing tests and calcu- 
lations indicate that the optimum shape for minimum 
drag will not be far from the shape consisting (1) of a 
secant ogive with twice the radius of the tangent ogive 
and (2) of a rounded meplat with a diameter of about 
0.1 caliber. The results also indicate that for heads 
longer than 3 calibers the head drag of the optimum 
shape will hardly be less than 50 per cent of the drag 
of a conical head with the same length. 


Q FUNCTION 


Tests in wind tunnels and firing ranges usually con- 
sist of measurements made at a series of discreet values 
of Mach Number. On the other hand, the designer is 
interested in a formula that will represent all measure- 
ments over the whole range of Mach numbers concerned. 
It is customary to derive such a formula from the theory 
and to evaluate the empirical constants involved by ex- 
periment. 

Unfortunately, the theory has not progressed far 
enough to treat completely the case of projectiles at 
supersonic speeds. By considering the variation of drag 
coefficient with Mach Number, the theory can handle 
the head pressure drag and the skin friction, at least 


* The solution of this differential equation is 


r ei ves Se OTE 2x\2) */? 
aE god ~i Se oh. ae ox 
R vaqit cna .-(2y} 


where r = radius of head at x and R = radius of body. 


t The meplat is loosely defined as the flat or blunt section at the 
front of the head. For example, most P.D. Fuzes have meplats. 
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Fic. 3. Drag coefficient for an infinite cone vs. cone’s semivertex 


angle. 


approximately, but it leaves the problem of base pres- 
sure completely unsolved. The ballistician has been 
required to depend so far on empirical formulas. 

Many formulas for representing Kp as a function of 
M have been devised. They are customarily a series of 
one type or another and are constructed for convenience 
in computing trajectories. They have the following 
common defects, however: (1) They are not satisfac- 
tory for extrapolating much beyond the range of Mach 
Numbers covered by experiment; (2) they can be used 
for only one shape of projectile. It would not be reason- 
able to expect a single formula to hold over the entire 
range from low subsonic to high supersonic velocities. 
The physical nature of the flow changes in going from 
subsonic to supersonic, and the speed of sound is a 
singular point in the theory. The transonic region must 
always be treated as a special case. However, there are 
no singularities past the speed of sound, and it should be 
possible to formulate rather general laws for the drag 
at supersonic velocities. 

In late 1944, R. N. Thomas, at the Laboratory, de- 
vised a drag coefficient function that works remarkably 
well for a broad class of projectiles in the supersonic 
region. Some hydraulic research by Dr. Hunter Rouset 
stimulated Thomas’ discovery. Rouse was studying the 
characteristics of flow past an object in an open water 
channel and sent Thomas some of his measurements. 
Thomas found on analyzing them that the square root 
of the pressure acting on the front of the object changed 
linearly with a hydraulic parameter corresponding to 
the Mach Number in gas flow. The corresponding rela- 
tion for projectiles is translated into terms of the drag 


coefficient by the equation 
— 4 Tv Ps; AM? " /a 
V Py = y" por { (522 — C+ ae) + Kor = 


where a and §# are empirical constants. 





a+ BM 
The first two 


t Iowa Institute of Hydraulic Research, University of Iowa. 
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SCALE COUMLE 


Fic. 4. Sketch of model projectiles: E-1, E-12, E-13. 
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Fic. 5. Drag coefficient vs. Mach Number for E-1, E-12, and 
E-13. 
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Fic. 6. Q function vs. Mach Number for E-1, E-12, and E-13. 


terms on the left-hand side under the square root bracket 
involve the base pressure and the skin friction. Un- 
fortunately, not only are the magnitudes of these quanti- 
ties not known accurately, but also they probably vary 
with Mach Number. Faced with this uncertainty 
Thomas suggested that these first two terms be re- 


placed by the number, unity. The equation then be- 
comes 


V1+ KpM? =a+ 06M 


Thomas designated the quantity V1 + KpM? by the 
letter Q and tested the fit of the Q function to the drag 
coefficient functions of all the standard projectiles, 
types 1 through 8.° The fit of the Q function to the ex- 
perimental data is remarkable, to say the least. Three 
examples of the Q function are given in Figs. 4, 5, and 6. 
In Fig. 5 the drag data are represented by the conven- 
tional Kp vs. M graph. In Fig. 6 the drag data are 
transformed into a Q vs. M graph. It can be seen that 
the Q function fits the measurements within the experi- 
mental scatter. 

The advantages of a linear function with only two 
coefficients that need be determined by experiment are 
self-evident. However, it should be remembered that 
the Q function is no more than an empirical law that has 
been tested by a limited body of data. Its justification 
is pragmatic rather than rational, and its derivation is 
intuitive rather than rigorous. Furthermore, it has 

certain definite limitations. It holds only over the 
supersonic region. It represents best the drags of sharp 
pointed, square-based shell. It is not suitable for blunt- 
headed projectiles, as proof slugs or spheres, and it is 
slightly in error for long boattailed projectiles. Analyti- 
cally, it is a special case of the inverse quadratic, which 
can be seen by solving for Kp, as follows: 
— py 2, @- 3) 
a, 
It is interesting to note that the Q function predicts the 
asymptotic value, b?, for Kp at infinite Mach Number. 


EFFECT OF BOATTAIL 


Projectiles are the supersonic twins of dirigibles, over- 
looking, for the moment, differences in the conventional 
methods of stabilization. The rear portion of a dirigible 
body is closed by a gentle curve to avoid separation of 


Fic. 7. Spark photograph of E-12M3 (E-12 bandless, 20 mm.). 
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the airflow and the high drag that accompanies a large 
wake. A similar attempt is made to streamline the rear 
of a projectile by adding, behind the body, a section 
tapered to the rear. This section, which is usually the 
frustrum of a cone, is called the boattail. Ordnance 
engineers have held the opinion that the advantages of 
boattailing are outweighed at supersonic speeds by im- 
provements in the head shape, since the overall length 
is ordinarily limited by the stability available from the 
twist of rifling in standard guns. Consequently, most 
artillery shell are either square-based or have only a 
short section of boattail (1/2 to 1 caliber long). 

The matter of boattailing is a debatable point, how- 
ever, since the data available to date are incomplete. 
A systematic program has been started at the Labora- 
tory to study this subject. Conical boattails are being 
investigated first, for most service shell use this shape 
and there are some tests to show that a conical shape 
gives a lower drag than an ogival shape (of the same 
length). Two parameters describe a conical boattail: 
the length from the body to the base; and the semi- 
vertex angle of the cone to which the frustrum belongs, 
known as the boattail angle. The first part of the test 
has consisted of a study of the effect of length on drag. 
The boattail angle has been fixed at 7° 15’, believed to 
be near the optimum. Three shells having a square 
base, a '/, caliber boattail, and a 1!/2 caliber boattail 
were fired at Mach Numbers from 1.3 to 2.5. 

The drag measurements from this test are shown in 
Figs. 5 and 6 and are discussed in part in reference 7. 
It should be mentioned that all drag measurements 
given here are corrected to the drag at zero yaw from 
the actual drag during flight. @Q functions were then 
determined for each projectile (see Fig. 4), and these 
have served as the basis for further analysis. 

In analyzing the test results, Dr. T. Hailperin* felt 
that it would be interesting to investigate the variation 
of drag with the area of the base, since the component 
of drag due to the base is equal to the product of the 
average base under pressure with the base area. He 
cross-plotted Fig. 5 into a graph with Kp as ordinate 
against the square of diameter of the base (in square 
calibers) as abscissa, with Mach Number as the inde- 
pendent parameter. The results are shown in Fig. 9. 
The model, E-13, has a slightly greater wetted area than 
the other two, and, as he was interested in pressure drag 
apart from skin friction, he made a small correction to 
the E-13 drag using the difference in wetted area and 
subsonic flat plate formula for the skin friction coeffi- 
cient (see Fig. 8). 

The corrected drag coefficients for the three models 
are fitted remarkably well by a straight line at any par- 
ticular Mach Number from 1.3 to 6.0 (see Fig. 9). 
Three points are hardly sufficient to establish with 
certainty that Kp varies linearly with (B/D). The 
fit to a straight line was so good, however, that Hail- 


* Ballistic Research: Laboratory, Aberdeen Proving Ground, 
Md. 
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Fic. 8. Skin friction drag coefficients for smooth flat plates at 
subsonic velocities. 
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perin felt justified to propose the linear variation as a 
working hypothesis and proceeded to investigate its 
physical significance. His analysis can be summarized 
briefly as follows: 

The linear hypothesis is expressed by the equation 


Kp = 1+ m(B/D)? (1) 


where / and m are functions of Mach Number but not of 
boattail length (also, they may change with boattail 
angle, but in the present discussion we are working al- 
ways with the fixed angle period of 7° 15’). Kp can be 
resolved into the components due to (1) the average 
pressure on the head, (2) the skin friction, (3) the aver- 
age pressure on the side of the boattail, and (4) the 
average pressure on the base. All pressures are based 
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on atmospheric pressures and are really pressure differ- 
ences from atmospheric (or free stream static, as the 
case may be) rather than absolute pressures. 


Kp = Kon + K +211 -(F)']4 
aes DH DSF pU? D 


4 
ba «(B\? 
be. T (5) (2) 
pU? 4\D 
Combining Eq. (1) with Eq. (2) gives a relation for the 


difference between the average side and base pressures, 
as follows: 





(pp = Ps) = (4/3) pU*m (3) 


The equation above does not contain any terms in- 
volving (B/D), and, in other words, the difference be- 
tween the average side and base pressures is independent 
of the boattail length. There are two possible physical 
interpretations. Either (1) the average side pressure 
and the average base pressure are each independently 
invariant with boattail length, or (2) there is a unique 
relation between the two satisfying Eq. (3). It has 
not been possible to date to say which of the two inter- 
pretations is correct. The independence and constancy 
of side and base pressure seem plausible, however, and 
will serve as a working hypothesis until more definitive 
information is available. This hypothesis can be de- 
veloped further to show that, if the average pressure is 
not a function of boattail length, then the individual 
pressures everywhere are also constant with respect to 
boattail length. 

On the basis of the foregoing, the results of the boat- 
tail program so far can be stated succinctly as follows: 

(1) The pressure on the base is constant and inde- 
pendent of boattail length. (The base pressure at the 
end of the boattail is the same as the pressure on a 
square-based projectile.) 

(2) The pressure on the side of the boattail is con- 
stant and independent of boattail length. 

(3) The difference between the side and base pres- 
sure is given by the equation 


Pe — Ps = (4/m)pU*m 


where pz, and ps are really under pressure from either 
atmospheric or free stream static. 


(3.1) 


SKIN FRICTION 


The principal components of drag are the head drag, 
base drag (including the boattail), and skin friction. 
The skin friction is the greatest unknown of the three. 
A theory has been developed for laminar boundary lay- 
ers over flat plates in supersonic flows, but its application 
to projectiles is limited, since the Reynolds Numbers are 
so high even for small models that most of the boundary 
layer will be turbulent. The state of the theory for 
turbulent skin friction seems somewhat uncertain. 
A few qualitative measurements have been made from 
the spin deceleration of artillery shell. The results 





give plausible values for the skin friction, but the 
method of deriving the skin friction from the spin de- 
celeration is not altogether sound. 

The lack of knowledge concerning skin friction has 
not greatly concerned ballisticians up to the present 
time, since friction accounts for only 10 to 20 per cent 
of the total drag of the artillery shell in current use. 
The picture changes as the aerodynamicist becomes in- 
terested in streamlined shapes designed for minimum 
drag without regard to overall length. Here the friction 
drag assumes a dominant role and may account for 
three quarters or more of the entire drag. Also, the 
balance between skin friction and pressure drag 
controls the design of the overall dimension of the 
shape. . 

The diffictilties in calculating the skin friction were 
discouraging at the start, but there was some encourage- 
ment in the possibility that subsonic friction formulas 
might work. Some evidence existed to show that 
boundary layers did not change radically in character 
as the external flow velocity increased from subsonic 
through sonic to supersonic. It was decided to try ordi- 
nary flat plate skin friction formulas. It would have 
been possible to calculate approximately the actual 
condition of flow adjacent to the surface over the pro- 
jectile. However, for the initial trial, the conditions of 
the undisturbed air ahead of the projectile and the pro- 
jectile’s velocity (or that of the free stream) were 
used. 

The supersonic wind tunnel of the Ballistic Research 
Laboratory had tested a 2-in. diameter model-of the E- 
12 bandless projectile at the supersonic Mach Number 
of 1.72. The tunnel measured both the total drag force 
and the base pressure. Subtracting the two gave a 
measured value of the drag component due to head 
pressure and skin friction. The same component of 
drag could be computed using Taylor and Maccoll for 
the head pressure and flat plate formulas for the skin 
friction. The theory of Taylor and Maccoll has been so 
thoroughly verified that the head pressure calculation 
must be reliable. Consequently, a comparison of the 
measured and computed values is a critical test of the 
application of flat plate skin friction formulas to the 
supersonic flow over projectiles. 

The wind tunnel measurements gave, at M = 1.72: 


Koy (head) + Kopsr (skin friction) = 0.055 
and the theoretical computation gave 
Kp (head) = (2/4)(pz/pU?) = 0.039 
Kosr (skin friction) = '/2(A/D*)C, = 0.017 
where C; is computed from the flat plate formula’® 
C, = [0.455/(logio R)***] — (1,700/R) 


and shown as a graph in Fig. 8, where R = length of 
projectile X U/v(= Ne). Thus, the computed value 
is 
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Koy (head) + Kpsr (skin friction) = 0.056 


As can be seen, the measured and computed values agree 
very well. 

A similar test of another conical headed model (E-21) 
was made in the wind tunnel. This second model had 
a 19° 58’ included angle conical head and a short cylin- 
drical body (0.29 calibers). The measured and com- 
puted values were compared and again the agreement 
was excellent (measured Kpy + Kpsr = 0.063, com- 
puted Kpy + Kpsry = 0.062). 

This analysis of wind-tunnel tests shows that subsonic 
flat plate formulas are satisfactory as a first approxima- 
tion for supersonic skin friction. The formulas given 
above’ (see Fig. 8) were used for all friction calculations 
made in this paper. On the other hand, it should be 
emphasized that the treatment of friction given above is 
at best an expedient serving in the lack of proper scien- 
tific knowledge. The condition of flow in supersonic 
boundary layers is undoubtedly different from that in 
subsonic. Further theoretical development in this field 
is necessary before the skin friction on supersonic pro- 
jectiles can be computed with confidence. 


EFFECT OF BANDING 


The great majority of artillery shell are stabilized in 
flight by spin, in a manner similar to tops and gyro- 
scopes. A band of soft metal is set into the body of 
the shell, usually near the base. It is sufficiently thick 
so that it has a greater diameter than the caliber of the 
shell and extends beyond the cylindrical surface of the 
body adjacent, and it is engraved by the rifling of the 
gun tube as the shell is fired. In this way it imparts the 
twist of the rifling to the projectile. 

In free flight the sections of the band that were not 
engraved project out into the flow over the body and set 
up disturbances that affect the aerodynamic properties. 
The drag is increased, and the overturning moment is in- 
creased or decreased depending on whether the band is 
ahead of or behind the center of gravity. The condi- 
tions of flow surrounding the band are complex, and it is 
not possible as yet to compute the banding effect from 
theory. 

However, it has proved feasible to fire bandless shell— 
smooth bodied projectiles—and spin them even through 
they have no rotating band. This is done by placing a 
short banded shell or ‘‘sabot’”’ behind the bandless pro- 
jectile in the gun tube. The sabot drives the projectile 
down the gun tube and falls away in free flight. Friction 
alone is adequate to impart full spin to the projectile. 

Tests were carried out with one type of projectile at 
M = 1.5 and with another type of projectile at M = 
3.2, The first tests consisted of firings of the E-12 and 
of the E-12M3, which is a 20-mm. bandless E-12 (see 
Fig. 4). A spark photograph of the E-12M3 is shown 
in Fig. 8. Neglecting any effect of scale, the band in- 
creases the Kp from 0.135 to 0.141. 
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Some pressure measurements made during the war in 
a German supersonic tunnel were available, and it was 
possible from these to estimate the additional drag due 
to the pressures acting on the front and rear of the band. 
These measurements compared with firing results 
showed that the pressure drag of the band alone ac- 
counted for almost all of the additional drag due to 
banding. If true, this result would mean that the band 
did not change the boattail and base drags materially 
and, in fact, could be considered as a kind of parasite 
projectile with its drag adding to the drag of its band- 
less host. 

The next step in the argument is to suggest that 
the change in band drag with velocity is governed 
by a Q function, since the band is equivalent to a para- 
site projectile. This hypothesis is expressed in mathe- 
matical terms by saying that, if the Q function of the 
banded shell is 


(1) Q=a+bM = V1 + KyM? 
and of the bandless shell is 
(2) Q.=a+)0,M = V1 + KpM? 


then the drag coefficient of the band is given by another 
Q function 


(3) Oz = de +. boM +- Vi + Kp,M? 


where Kp, is the additional drag of the band, and the 
three equations are related by 


(4) Kp = Kp, + Kp, 


The values of a and 6 and of Kp, Kp,, and Kp, for M = 
1.5 had been determined by experiment. This left the 
four unknowns 4, };, d2, and 5; to be solved by the four 
simultaneous equations above. The solution showed 
that b. = 0 (closely) and 5; = 8, so that the Q function 
of the band is given by simply Q = constant, and that 
the slopes of the banded and bandless Q function are 
equal. 

Physically, this result means that the effect of the 
band on drag decreases as the velocity is increased. 
The band drag was calculated from the first test as dis- 
cussed above at the Mach Number of the second test, 
M = 3.2, and compared with measured value from the 
second test. The agreement is satisfactory, being 
within the experimental error. The check is open to 
some question, however, since the shapes of the projec- 
tiles in the two tests are very different. 

The analysis of banding given here appears plausible 
and consistent with the data available. On the other 
hand, it is the author’s opinion that the analysis should 
be accepted with some reservation. The experimental 
measurements are far from complete, and more firings 
must be conducted before any final conclusions can be 
drawn. 





OTHER EFFECTS 


Spin 

All of the aerodynamic properties of spin-stabilized 
shell may be functions of the spin, as well as the veloc- 
ity, shape, air density, etc. Certain of the forces are 
known to vary with the spin. However, ballisticians 
commonly assume that the drag force is not affected by 
the spin. There are some measurements made on spin- 
ning models in German wind tunnels (during the war) 
which indicate that the effect of spin on drag is negli- 
gible for the spin rates produced by most service 
weapons.* All drag measurements given in this paper 
will be assumed to be independent of spin. 


Scale 


The scale is another parameter entering the picture, 
usually as the Reynolds Number. For bodies at sub- 
sonic speeds it is commonly assumed that the skin fric- 
tion alone is a function of the Reynolds Number and 
that pressure forces are independent of R except insofar 
as they are influenced by boundary phenomena, such as 
separation. This assumption will be used as a working 
hypothesis for the supersonic bodies discussed in this 
paper. 

Some test data are available on the effect of scale in 
firings of full-scale and model 155-mm. artillery shell." 
Unfortunately, the drag of the full-scale projectile was 
about 10 per cent higher at supersonic velocities than 
the drag of the caliber 0.50 model, and this measurement 
is contrary to predictions based on the hypothesis 
stated above. The 155-mm. results may be modified 
somewhat by future firings made under more carefully 
controlled conditions, since the models were not exact 
reproductions of the service shell and the relative rough- 
nesses of the two were not measured. It is commonly 
assumed that the effect of scale is small at supersonic 
velocities and the 155-mm. results bear this assumption 
out. Scale is another of those factors whose effect is not 
well known but is believed to be small. 


Effect of Head Shape on Conditions at Base 


The head and base are treated as independent quanti- 
ties in this paper. This is certainly true for the head, 
since change in conditions at the base cannot affect the 
head if the flow over the shell is supersonic. On the 
other hand, the base is subject to changes in the flow 
over the head, and the treatment of the base indepen- 
dent from the head is not rigorous. However, it is prob- 
ably a good approximation if the body is longer than a 
caliber. Some tests bearing indirectly on head-base 
effect show that the drag is roughly independent of body 
length for bodies one to three calibers long.'' Some 
preliminary firing data indicate that changes in drag 


* The majority of guns are rifled with twists ranging in value 
from one turn in 25 calibers to one turn in 40 calibers; one turn 
in 30 is most common. 
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produced by small changes in cone angle are accounted 
for entirely by head pressure. In fact, from a qualita- 
tive standpoint, the flows over conical heads with 
slightly different cone angles should reach about the 
same condition over the body. As the head angle is de- 
creased, the pressure on the head is decreased, but at the 
same time the drop in pressure around the shoulder to 
the body is decreased, so the two should tend to cancel. 
The independence of base from head is again a working 
hypothesis that should be accurate to a first approxima- 
tion. 


DESIGN OF SUPERSONIC AIRCRAFT FUSELAGE 


Design Conditions 


To illustrate their practical application, the empirical 
formulas that have been given in the foregoing will now 
be used to design the fuselage of a supersonic aircraft. 
The purpose of the design will be to develop the shape 
having a minimum drag for a given body diameter under 
the flight conditions to be specified. No consideration 
will be given to methods of propulsion or to lifting and 
stabilizing surfaces. The fuselage designed here will be 
a simple body of revolution with both the front and rear 
ends closed surfaces. 

The head and boattail will be designed separately and 
will be considered independent insofar as pressure ef- 
fects are concerned. All pressure forces will be based on 
pressure differences from atmospheric so that head and 
base pressure forces will both act opposite to the diree- 
tion of motion. Furthermore, the pressure forces will be 
assumed independent of scale—that is, Reynolds Num- 
ber. The effect of scale will enter only in the skin fric- 
tion drag. The skin friction force will be computed 
from the total wetted area and from flat plate skin fric- 
tion formulas and a Reynolds Number based on the 
overall length. The total drag will then be computed 
by adding head, base, and skin friction drag com- 
ponents. 


Head 


The head shape is arbitrarily chosen to be a cone. 
This is probably not the best head shape for any given 
length, as discussed previously, but the improvement 
that could be realized by using a better shape would 
hardly be significant insofar as the total drag is con- 
cerned. The head drag (pressure component) is scarcely 
more than 10 per cent of the total, as will be shown, so 
that even a 50 per cent improvement would only de- 
crease the total drag by 5 percent. The choice of a cone 
has the great advantage that the pressure drag can be 
computed from theory, and, since the experimental 
data on other head shapes are too incomplete to develop 
any general empirical formula, use of the Taylor and 
Maccoll theory is essential at this point. 

As the head is made longer, the cone angle becomes 
smaller and the pressure drag is reduced. On the other 
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hand, increasing the head length increases the wetted 
area and so increases the skin friction drag. The opti- 
mum length is that length which balances head pressure 
against skin friction so that the added drag of the two 
is a minimum. The head pressures were computed 
from Taylor and Maccoll (Fig. 3) and the skin friction 
coefficients from the flat plate formula’ (Fig. 8). The 
head drag coefficient resulting from both pressure and 
skin friction is shown plotted against the cone’s semi- 
vertex angle in Fig. 10. The nature of the minimum is 
clearly evident. In order to determine the exact mini- 
mum point, a quadratic was fitted by least squares to 
the five computed values of Kp, nearest the minimum, 
and the value of 6, for the minimum was computed from 
the coefficients of the quadratic. 


Body 

The body is a circular cylinder 2.14 calibers long. 
The choice of length may appear somewhat arbitrary, 
but it is the body length of the 20-mm. E-12M3 model 
(tested at the Ballistic Research Laboratory) used in 
subsequent base pressure computations (see Figs. 4 and 
7), and it gives a reasonable appearance to the final 
design. Inclined flight is not considered, so the body 
has no pressure drag and only contributes wetted area 
to the total skin friction drag. Also, the body is for an 
aircraft and does not have a rotating band. This factor 
appears in the calculation of base pressure from ballistic 
measurements of projectiles. 


Boattail 


The design of boattail and base shapes is not so 
simple as the design of head shape. Again, a conical 
shape is chosen. No theory similar to Taylor and Mac- 
coll has been developed for conical boattails, but nearly 
all the experimental studies have been made on boat- 
tails with a conical shape (or frustrum of a cone, de- 
pending on the length). Here, however, the argument 
for a conical boattail as the best shape is somewhat bet- 
ter than in the head shape case, for there are a few 
measurements that indicate that conical boattails give 
a lower drag than ogival boattails. 

The conical shape is defined by two parameters, 
angle and length, and the optimum combination will be 
determined by a balance of side pressure drag, base 
pressure drag, and skin friction. _It is difficult to visu- 
alize the combined effect of two parameters changing si- 
multaneously, so it will serve initially to consider the 
length fixed at some arbitrary value and the angle 
varied. As the boattail angle is increased from 0°, the 
pressure on the side will decrease from its value at 0° 
(the value of the pressure on the body). The wetted 
area will also decrease and with it the skin friction. 
However, this effect will be small at small angles com- 
pared to the pressure effects and can be neglected in the 
present discussion. On the other hand, the pressure on 
the base (defined here as the flat at the rear of the body) 
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Fic. 10. Head drag coefficient for both skin friction and pressure 
drag vs. cone’s semivertex angle. 


will remain constant.* If one looks at the projected 
cross-section area of boattail and base, he sees a central 
circular area acted on by the base pressure and a sur- 
rounding ring-shaped area acted on by the side pressure. 
At small angles the side under pressure is less than the 
base under pressure, and the total base drag resulting 
from both boattail and base pressures will be less than 
the base drag of a square-based shell. As the angle is 
increased further, the projected boattail area increases 
but the side pressure decreases. If the angle is in- 
cfeased far enough, the side under pressure will increase 
to the value of the base under pressure, providing the 
flow still follows the boattail surface. At this angle the 
total base drag will equal the square base drag. It 
is clear that there is an angle between 0° and 
the angle for ps = pg giving a minimum total base 
drag. 

Unfortunately, for the purposes of this design, the | 
law giving the variation of side pressure with boattail 
angle is not known. All of the firings reduced to date 
from the boattail programs are for the single angle of 
7° 15’. In organizing the basic boattail program, this 
angle was selected for the initial study of length effect, 
because (1) the flow does not separate and (2) older 
firing tests showed that the boattail angle for minimum 
drag lay between 5° and 9°. Judging by the older 
firings, the angle of 7° 15’ would not be far from the 
optimum. Accordingly, the boattail angle was fixed at 
7° 15’ for the present fuselage design. 

* It has been shown that the base pressure does not change with 
length for an angle of 7° 15’, and it seems reasonable to believe 
that this relation will hold at other small angles. The base pres- 
sure is therefore equal to the pressure on a square-based body and 
is invariant with both length and angle. 
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Fic. 11. Change in drag coefficient due to boattail: increase due 
to skin friction; decrease due to pressure reduction. 


The length is the only parameter left, since the angle 
has been fixed at 7° 15’. As the length is increased, the 
boattail results show that the base drag is decreased ac- 
cording to the formula 


Kp = 1+ m(B/D)? 


Hence, the reduction in base drag coefficient, AKpp, 
from pressure drag only, due to the addition of a boat- 
tail having a base diameter, B, to a square-based shell, is 


AKpp = m{1 — (B/D)?] 


To balance this, the wetted area and skin friction drag 


increase as the boattail is lengthened. The increase in . 


the base drag coefficient, AK ps,, from skin friction alone, 


esse 


where £ is the boattail angle and C; is based to a close 
approximation on the Reynolds Number computed on 
the overall length of the whole body. 

So the resultant reduction in base drag coefficient 
from adding a boattail to a square-based shell is 


rf C. B\? 
ato=|m— 5s) ~ (5) 
It is clear from this formula that boattailing is an all-or- 
none proposition. If 
m > (/8)(C;/sin B) 
the boattail should be extended to a sharp point. If 


AKpsr - 


m < (x/8)(C;,/sin 8) 
no-boattail at all should be added. At M = 1.5, 










JOURNAL OF THE AERONAUTICAL SCIENCES—MARCH, 1947 





m > (/8)(C;,/sin 7° 15’) 
and a completely pointed boattail is called for. 


As the 
Mach Number increases, both m and C; decrease, but m 


decreases much more rapidly than C;. The values for 
AKpp and AKps, are shown plotted against the Mach 
Number in Fig. 11. The two curves approach one 
another, but out to a Mach Number of 6 the pressure 
contribution is still greater than the skin friction con- 
tribution. 

The highest Mach Number considered in the present 
design is 5. In view of the considerations discussed, the 
boattail was fixed at an angle of 7° 15’ and was ex- 
tended to a sharp point. 


Base Pressure 


The graph of AKpp in Fig. 11 gives only the reduction 
in drag realized by boattailing. It remains to calculate 
the base pressure. It is assumed that head and base 
conditions are independent, and, consequently, the base 
pressure for the fuselage is the same as the base pressure 
for the bandless E-12M3. The total drag for the E- 
12M3 is given by the Q function (see Fig. 4). 


Q = 0.9462 + 0.1304M (= V1 + KpM?) 


A spark photograph of the E-12M3 is shown in Fig. 7. 
The head pressure is given by Taylor and Maccoll (see 
Fig. 2) and skin friction by subsonic flat plate formulas 
(see Fig. 8). The base drag is then computed by sub- 
tracting the head pressure drag and skin friction drag 
from the total drag. 


Kpsp = Kp — Kon — Kose 
The nondimensional base pressure is related to the 
base drag coefficient by the equation 


Pe _ (4/x) yM*Kpp 


pi 

The quantity referred to so far as the base pressure is 
really the under pressure from the free stream pressure, 
pi. The absolute base pressure, Pz, is given by the 
equation 

Py = fi — bp 
and, in terms of Kpz, by 
(Pp/pi) = 1 — (A/w)yM*Kop 


P,/p1 is shown as a function of Mach Number in Fig. 12. 

The fuselage has no actual base area, of course, since 
the boattail is brought to a sharp point. The drag on 
the boattail is computed by combining the square-base 
drag, Fig. 12, with the drag differential from square base 
to complete boattail, Fig. 11. The formula is, simply: 


* es (boattail) = Kp, (square base) — AKpp (Fig. 11). 


The boattail drag will be subsequently referred to as 
the ‘‘base drag” in discussing the performance charac- 
teristics of the various fuselage designs, but it will be 
understood that the term applies to the pressure drag 
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Fic. 12. Base pressure vs. Mach Number for E-12M3. 
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SUPERSONIC FUSELAGE 
Fic. 13. Sketch of aircraft fuselage and its models. 


acting on the rear section of the fuselage, whatever the 
shape may be. 


Skin Friction 


The skin friction is computed from the total wetted 
area and a Reynolds Number based on total length and 
free stream conditions as discussed previously. The 
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total length depends, of course, on the optimum cone 
angle for the head. 


Resolution of Drag Components 


The total drag is now computed by adding the head 
pressure drag for the optimum cone angle, the base pres- 
sure drag for the boattail (or base, as specified), and the 
skin friction 


Kp = Kon + Kosr + Kop 


The optimum cone angle will depend on both Mach 
Number and size, so the fuselage head shape will be 
slightly different for different design conditions. All 
drag components are functions of Mach Number, and, 
furthermore, the head drag (as influenced by the cone 
angle) and skin friction will also depend on size. Con- 
sequently, the overall aerodynamic performance will 
differ for each particular size and each flight condition. 


Fuselage Design 


The fuselage design is determined by (1) the Mach 
Number, (2) the meteorological conditions, and (3) the 
body diameter. Discussions with aeronautical engi- 
neers led to the consideration of a 50-in. diameter fuse- 
lage flying at M = 1.5 at 50,000 ft. as a design that 
might be feasible for actual flight. The results of the 
design for minimum drag are shown in Figs. 13 and 14. 

The most striking feature of the design is the long 
pointed head required for the optimum balance of pres- 
sure and skin friction. The relative contributions of the 
head pressure, base pressure, and skin friction to the 

















SUPERSONIC FUSELAGE 

50" Fuselage 20mm Model 20mm Model 20mm Model 20mm Model 

Boat tail Boattail Boattail Square Base Square Base 
50,000 ft.* Sea Level Sea Level Sea Level Sea Level 
Ms 1.5 Ms 1.5 M= 5.0 Ms 5.0 Ms 5.0 
@, = 3.54° @ = 4.04° & = 4,22° @= 4.22° @= 4.22° 
Lng. — Lng. _ Lng. Ing. _ g, Ing. _ 
 llaiees B. = 33-2 =? * Aid Wa; = 12-9 
Kp = 0.066 Kp = 0.075 Kp = 0.038 Kn = 0.047 Ky = 0.061 
Head. scccccceellS Head...ceccseccl3h Headecseceesel5& Healle.scccscesleh Bodh cccccccccese Gh 
Base... sees 000 48S pe eT Bas@..ccsccceel$ a RE Bas@.cecees oe ee 30% 
Skin FrecccoectlZ Skin Fre ccccect6% Skin Phocccceeyey Skin Fr...ceee Skin WeeccceccceS 
“Power required > 2,263 hep. (at sea level = 21,243 h.p.). 





Fic. 14. Aerodynamic performance table. 
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drag are shown in Fig. 14, and the effect of the long 
slender head is evident. The head pressure assumes a 


minor role in comparison to the skin friction and base _ 


pressure. It is interesting to note that at this Mach 
Number, 1.5, the base pressure accounts for nearly one- 
half the drag, even though the projectile has a complete 
boattail. 

The power was computed for flight at altitude and 
amounts to 2,263 hp. The price exacted for supersonic 
flight can be seen in this tremendous power required 
to drive a body only 50 in. in diameter. The figure of 
2,263 hp. is trifling, however, compared to the power re- 
quired for sea-level flight of the 50-in. fuselage at M = 
1.5, that is 21,243 hp. Nearly ten times the power will 
be needed to maintain speed as the aircraft descends 
from altitude to fly near the ground. This figure bears 
out the opinion that supersonic flight at low altitudes is 
impractical with the fuels and power plants that are 
available today. 

The development of any new type of airplane is usu- 
ally accompanied by scale-model tests. In fact, the de- 
sign of the full-scale airplane is often based on the mod- 
el’s performance, and it is interesting to see the opti- 
mum design that would result from sea-level tests of a 
small (20-mm.) model and to compare the performance 
of the two. The results, given in Figs. 13 and 14, show 
a striking similarity between optimum model and opti- 
mum fuselage. The drag coefficient of the model is 
somewhat higher, and the cone angle is slightly greater. 
The relative proportion of head, base, and skin friction 
drags are about the same. This comparison shows that 
the practical designer would find small scale-model tests 
at sea level satisfactory for determining the performance 
of the full-scale airplane at altitude. 

The aircraft designer is interested in performance at 
higher and higher speeds, so, next, the optimum design 
was computed for a 20-mm. model at a Mach Number of 
5. The results show that the best head angle is slightly 
steeper and that the drag coefficient has decreased to 
about one-half the value at M = 1.5 (see Figs. 13 and 
14). The most striking difference is the change in the 
proportions of head, base, and skin friction drag. The 
skin friction now takes the dominant role, with head 
drag about the same and base drag almost negligible. 
The skin friction accounting for most of the drag re- 
sembles conditions at subsonic speeds, but, of course, 
the optimum shapes at M = 5.0 and subsonic bear little 
resemblance. The importance of skin friction in this 
design emphasizes the necessity for the development of 
boundary-layer theory for turbulent layers at high 
supersonic speeds. 

Ballisticians have asserted for many years that the 
advantages of boattailing decreased with velocity (in the 
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supersonic region). In order to investigate this point, 
the performance of two square-based models was com- 
puted at M = 5.0. The first model has the same head 
and body as the optimum shape but the boattail is 
omitted entirely by ending the body with a square base 
at the point where the boattail would start (see Fig. 13). 
The results show that the boattail is still effective at 
M = 5.0. Leaving the boattail off altogether increases 
the drag coefficient 24 per cent. 

The effect of skin friction can hardly be neglected 
when comparing different models at M = 5, so the per- 
formance was computed for a second square-based 
model having the same overall length as the optimum 
model (see Figs. 13 and 14). The additional skin fric- 
tion raises the drag increase from 24 to 61 per cent. 
The practical square-based design* would probably lie 
somewhere between the two cases considered here. The 
advantage of boattailing, insofar as the drag is con- 
cerned, is evident from the increase of 24 to 61 per 
cent in drag that will result if boattailing is not 
used. 
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Beam-Columns 


WM. R. OSGOOD* 
Navy Department 


ABSTRACT 


A method of designing beam-columns is presented which 
appears to be both safe and economical. The method consists of 
two steps. (1) The given transverse load is replaced by two 
equal and opposite couples equal to the maximum bending 
moment produced by the transverse load acting alone. The 
two couples act at such sections that they produce the same 
maximum deflection at the same section as the given transverse 
load acting alone. (2) The two couples and the given axial 
load are now replaced by eccentrically applied loads at the ends 
which produce the same maximum bending moment as the 
couples and the axial load acting alone. 


INTRODUCTION 


Mo ENGINEERS are probably not completely satis- 
fied with current methods of designing beam- 
columns—that is, members carrying both axial and 
transverse loads or, more precisely, heterostatic loads 
(see below). An approximate method is proposed here 
which is believed to be safe and economical, although 
only extensive tests can verify the belief. The pres- 
entation is confined to members simply supported at 
the ends and to cantilevers. It could probably be ex- 
tended to members overhanging one or both supports, 
but such members as beam-columns are perhaps not too 
important. Statically indeterminate members are 
important but are not considered in this paper. 

The author has felt for some time that the most hope- 
ful approach to the beam-column problem is to try to 
replace the beam-column for purposes of design by an 
equivalent eccentrically loaded column. The latter isa 
special beam-column that can be studied accurately 
and coimparatively easily in the laboratory. It is 
assumed that, as a result of laboratory tests or other- 
wise, a formula or a chart is available relating load, 
length, and eccentricity for columns of a given material. 
Then the equivalent eccentrically loaded column can be 
referred to the formula or the chart for design. A 
formula that has received substantiation by some tests 
conducted at the National Bureau of Standards for the 


‘ Bureau of Aeronautics, Navy Department, is 


A = Ao/(1 + 4p?) (1) 
where 


(1/m)(l/i) VS/E 
ec/1? 


P/(AS) 


o 
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E is Young’s modulus, S is the compressive secant yield 
strength (slope of secant was taken as ?/;E for materials 
of tests); / is the length of the column; 7 is the least 
radius of gyration of the cross-sectional area in the 
plane of bending; e¢ is the eccentricity of the load; c is 
the half-depth of the column in the plane of bending 
(symmetry with respect to an axial plane normal to the 
plane of bending is assumed); A is the cross-sec- 
tional area; VP is the load at failure; and Xp is the 
value of \ for a given o when p = 0 (centrally loaded 
column). 


DISCUSSION OF LOADS 


In discussing beam-columns it will be convenient to 
use the terms orthostatic, astatic, and heterostatic 
introduced by Westergaard.' In particular, ortho- 
static loads are loads under which the law of super- 
position applies. In beams shears, moments, slopes, 
and deflections are all linear functions of the (ortho- 
static) loads. Loads that, when increased in the same 
proportion up to critical values, produce pure buckling 
are astatic loads. The law of superposition does not 
apply. Ina centrally loaded column shears, moments, 
slopes, and deflections are not linear functions of the 
(astatic) load. Heterostatic loading is loading that 
can be resolved into two component sets of loads: one, 
orthostatic; the other, astatic. The law of super- 
position does not apply. LEccentrically applied loads 
on columns are an example of heterostatic loading. The 
loads can be resolved into couples and transverse reac- 
tions at the ends (orthostatic loads) and axial loads at 
the ends (astatic loads). 


PROCEDURE 


A crude, safe, but in general uneconomical eccentri- 
cally loaded column equivalent to a beam-column is 
one with an equivalent eccentricity, e,, of load 


QE = M;/P (2) 


where i, is the maximum orthostatic bending moment 
and P isthe axialload. In designing, one would assume 
a section, enter the formula or chart with the given 
length of the member and the axial working load multi- 
plied by the factor of safety, and compare the eccen- 
tricity resulting from the formula or chart with e,. If 
the former were greater than the latter, the assumed 
section would be safe but uneconomical, and another 
section should be tried. 
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It is easy to see that the preceding procedure would 
in general result in an uneconomical design because the 
given orthostatic loading is replaced by a much more 
severe one—namely, couples at the ends that subject 
all cross sections of the member to the maximum bend- 
ing moment, M,, and increase the deflections. In 
order to improve economy without doing violence to 
safety, the following method of determining an equiva- 
lent eccentricity is proposed for consideration. 


DETERMINATION OF AN EQUIVALENT ECCENTRICITY 


The procedure consists of two steps. First, the given 
orthostatic loads are replaced by orthostatic loads con- 
sisting of two equal and opposite couples equal to the 
maximum orthostatic bending moment, the couples 
being applied at such sections that they produce the 
same maximum deflection at the same section as the 
given orthostatic loads. The maximum orthostatic 
bending moment and deflection are therefore the same 
under the new loads consisting of two couples as under 
the given loads. Now, the new heterostatic loading, 
consisting of the axial loads and the two couples, is re- 
placed by eccentric loads at the ends that produce the 
same maximum heterostatic bending moment as the 
new heterostatic loading. The problem consists in 
locating the couples and in determining the equivalent 
eccentricity of the loads at the ends. 

We will determine first the sections at which the 
couples replacing the given orthostatic loads act. 
With the notation in Fig. 1 for a member simply sup- 
ported at the ends, by the ordinary methods of the 
strength of materials, when P = 0 the deflection is a 
maximum at the section, for which 


ky = */(1 + ki? — hp?) (3) 


and is equal to 





A = 1/.(Mgl*?/EI) (ky? — hi*) (4) 
Solution of these equations for k; and k, gives 
ki = Wks? — 2[EI/(Mzl*)] A (5) 





ke = V(1 — kg)? — 2[EI/(Mol*)] A (6) 


In the case of symmetrical loading for which k; = ky, = k 
Eqs. (5) and (6) reduce to 





k = V (3/4) — 2[EI/(Mgl*)] A (7) 
Eqs. (4) to (7) assume elastic action. Consideration 
will be given to plastic behavior later. When ky, A, 


and the maximum bending moment M, for the given 
structural member are substituted in the right-hand 
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sides of Eqs. (5) and (6), or (7), the resulting values of 
k; and kp, or k, determine the new orthostatic loading to 
produce the same maximum moment and deflection as 
the given orthostatic loading. Now, if the axial load P 
is applied, one would expect the beam-column with the 
substitute loading to behave in much the same way as 
with the original loading. 

The second step consists in evaluating an equivalent 
eccentricity for the beam-column with the substitute 
loading. The loads resulting from the substitution are 
shown acting on the member in Fig. 1. So long as 
ki S '/2and ke S 1/2, the deflection is a maximum at the 
section for which 


cos akyl — cos ak! cos al 





t kl = 
— cos akl sin al 


where a = V/ P/EI and the maximum moment, at 
this section, is 

Mwez. = Mpcsc al X 

V cos? akil — 2 cos aki cos aki cos al + cos? ak (8) 





This same moment can be produced by the forces at 
the ends alone applied with proper eccentricities, @,, 


Mer, = V 2Pe, esc al V1 — cos al (9) 


We take as the equivalent eccentricity we are looking for the value e, obtained by equating Mme:.. from Eqs. 


(8) and (9): 








_ Ms .jcos* aki — 2 cos aki cos ak, cos al + cos? ake) 
¢ = — (10) 
V/2P 


1 — cos al 
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When k; = ke = k, Eqs. (10) and (11) reduce to 
é. = (Mz/P) cos akl (12) 


Design now proceeds as explained in connection with 
Eq. (2). If ki > */2 or ke > '/2, the maximum hetero- 
static moment and the maximum heterostatic deflection 
would not in general, if ever, occur at the same section; 
and unless they do, the method loses much of its possible 
value, because the equivalent eccentrically loaded 
column would probably behave quite differently from 
the given beam-column. 

That the method proposed is more economical than 
that embodied in Eq. (2) is evident from the fact that 
the multipliers of M,/P in Eqs. (10) to (12) never ex- 
ceed 1. The question of safety cannot be answered 
simply, because most of the equations involved assume 
elastic action, whereas at failure the action is not elastic. 
Only tests can give the final answer, but there are some 
considerations that indicate safety. Let us look at 
them. 


CONSIDERATIONS OF SAFETY 


Eq. (3) is valid in the plastic range as well as in the 
elastic. In the plastic range, E in Eq. (4) would have 
to be replaced by a reduced modulus F,, which could be 
done accurately, for example, by the procedure of 
references 2 and 3. This E, would appear instead of E 
in Eqs. (5) to (7). The A in these equations, which is 
the maximum orthostatic deflection of the given mem- 
ber, could also be computed by the procedure of refer- 
ences 2 and 3. It would be greater than under the 
assumption of elasticity and would be inversely propor- 
tional to some other reduced modulus, £,, since in the 
elastic range it is inversely proportional to E. If FE, = 
E,, the second term in each of the radicals of Eqs. (5) to 
(7) is accurate as it stands when A is computed on the 
assumption of elasticity. The first term in the radical 
of Eq. (7) is correct in any case, but the first term in 
each of the radicals of Eqs. (5) and (6) may be expected 
to be slightly different in the plastic range, since it de- 
fines the distance to the section of maximum deflection 
in the given structural member. Looking at the overall 
picture, one has the feeling that the loading of the sub- 
stitute beam will not be much in error when the ’s, 
locating the sections at which the couples Mz act, are 
computed from Eqs. (5) to (7) as they stand. 

The k’s, just discussed, enter in the final Eqs. (10) to 
(12), which are likely to be insensitive to all but gross 
errors in k. The greatest error committed by these 
equations if applied in the plastic range is due to the E 
in the expression for a. This error can be reduced 


somewhat by substituting the tangent modulus Z’ or 


é. = (Mz/P) V cos aki cos ak + [(cos aki — cos aks) /(2 sin al/2)]? (11) 





the double-modulus E for E. Fortunately, the error is 
on the safe side in any case. The effective modulus 
that should actually take the place of E is smaller than 
E and would lead to a smaller value of e, than the value 
computed from the equations as they stand. 

There is, of course, an assumption involved in the 
development that will bear questioning—namely, the 
implication that failure occurs at the same maximum 
moment for the eccentrically loaded column of Eq. (9) 
as for the column of Fig. 1 and Eq. (8). The smaller 
the sum of the k’s, the better this assumption. The ice 
does not seem too thin here to the author, however, but 
let it be said again the proof of the pudding lies in tests 
before eating, or before swallowing whole at any rate. 


EXTENSION TO CANTILEVERS 


The extension to cantilevers is fairly obvious when 
the cantilever is considered as one-half of a symmetri- 
cally loaded member twice as long as the cantilever. 
Eqs. (7) and (12) become 


ke = V1 — 2[EI/(Mpl2)) A (13) 
é = (M3/P) cos ak, (14) 


where /, is the length of the cantilever, k,/, is the distance 
from the free end to the section at which M, acts in the 
substitute cantilever, and A is the deflection of the free 
end of the cantilever. 

Finally, before playing with the equations a little, 
perhaps it should be noted that exceptional cases of 
loading can no doubt arise to which the method is in- 
applicable, directly at least. It might be necessary to 
introduce more complicated orthostatic substitute loads 
to handle such cases, but this could probably be done. 


APPLICATIONS OF THE METHOD 


A few applications of the method may be of interest. 
The loads indicated are the orthostatic loads. The 
total transverse load is denoted by W. 

(1) Simple supports at ends—concentrated load at 
middle: 


Mz = ('/)Wl, A = (1/48)(WI/ED) 
k= V1/, — 1/, = 1/2 V3 = 0.2887 
& = (W1/4P) cos [I/(2 V3)]V P/(ED) 


(2) Simple supports at ends—uniformly distributed 
load: 


Mz = '/sWiI, A = (5/384) [WI*/(ED] 





k 


ll 


V (1/4) — (5/24) = 1/2V6 = 0.2041 
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& = (WI/(8P)] cos [1/(2 V6)] VP/(ED 


(3) Simple supports at ends—load increasing uni- 
formly from zero at the ends to a maximum at the mid- 
dle: 

Mz = 1/,WI, A = (1/60) (W/*/(ED] 
k= V1/,—1/, = 1/(2 V5) = 0.2236 


[W1/(6P)] cos [1/(2 W5)] VP/(ED 


e 


(4) Simple supports at ends—load increasing uni- 
formly from zero at one end to a maximum at the other 


(k,l is measured from the end at which the load is 
zero): 


ka = V1 —2-V2/15, Mz = 2W1/9 V3) 











2 ee 3 SR reer 
k= @#2(1-4@-— - 1-242 i Ay? = 0.1664 
: y ( \" 1 15) 5 N57 25 Vis s 








9 


owl ft P P 
& = 7=— @/cos 0.25761 @/ — cos 0.16641 
9V3P EI EI 


.. (5) Cantilever—concentrated load at free end: 


Mz = WI., A = '/;[W13/(ED] 


k, = V1 — 3/, = 1/V'3 = 0.5774 
@ = (WI,/P) cos (1,/W3) VP/(ED 


(6) Cantilever—uniformly distributed load: 
Mz = '/2Wl., A = '/;(W18/(ED)] 


ke = V1 —1/, = 1/V2 = 0.7071 
e, = [W1./(2P)] cos (I./V2) VP/(ED 


(7) Cantilever—load increasing uniformly from zero 
at free end to a maximum at the fixed end: 





. (= 0.25761 V P/(ED) — cos 0.16641 V Jy 


2sin (1/2) ~P/(EI) 





Mp =1/3sWl., A = (1/15)(W18/(ED)] 


k= V1 — 2/, = V8/, = 0.7746 


/ 


(W1./(3P)] cosl. WV 3/5 WV P/(ED 


Ce 
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An Application of IBM Machines to the 
Solution of the Flutter Determinant 






E. L. LEPPERT, JR.* 
Lockheed Aircraft Corporation 


ABSTRACT 


An electrical “punch card”’ procedure for the formation and for 
the solution of the complex flutter determinant has been de- 
veloped, using an adaptation of W. M. Bleakney’s trial-and-error 
method.2 The procedure is of particular value in obtaining 
answers for a complicated problem in a minimum of elapsed time 
and greatly reduces the possibility of human errors caused by 
fatigue and strain. The effects of variation of parameters are 
readily obtained. Amplitude ratios may be calculated with a 
small number of additional operations. 


INTRODUCTION 


DD cpaserss THE PAST SEVERAL YEARS the increase in size 
and speed of the modern airplane has resulted in 
the increased importance of a complete and rapid 
flutter analysis of the aircraft structure. The increase 
in speed makes necessary an investigation of almost 


‘ every aerodynamic surface in order to insure the safety 


of the airplane throughout its speed range. The in- 
crease in size generally entails the desirability of in- 
creased degrees of freedom to adequately define the 
possible flutter motions. All factors are reflected in a 
greatly increased volume of computations. 

Sufficient analysis must be carried along simultane- 
ously with the preliminary design to eliminate the 
possibility of costly changes after the general structural 
design has been determined. A method is desirable by 
which a change in any of the parameters under control 
of the designer can be readily evaluated. Computa- 
tions and results should be of such a nature as to be 
readily checked and quickly verified by the licensing or 
purchasing agencies. 

Following the design and construction of the proto- 
type airplane, changes and improvements are fre- 
quently made which necessitate an immediate evalua- 
tion of their effects on the flutter characteristics. Ade- 
quate answers must be obtained in a minimum of 
time. 

Although considerable information can be obtained 
from investigations using two or three degrees of free- 
dom to represent a more complex physical problem, 
wide experience and much ingenuity are required to 
assure adequate similarity. A computation method 
that allows the solution of a relatively large number of 
degrees of freedom quickly, economically, and accu- 
rately is desirable. 
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A trial-and-error method of analysis, incorporating 
most of the desirable features mentioned above, is 
described in a previously published paper.?2 The IBM 
machine method herein described was designed to 
follow as closely as possible along the line of attack 
described in that paper. 

Tabulating machines have previously been used to 

calculate the three-dimensional mechanical and aero- 
dynamic terms that make up the flutter determinant. 
Iteration methods* ° have been successfully applied to 
its solution but do not satisfy all the desired features 
previously mentioned. The following application of 
IBM equipment retains most of the advantages of the 
“starring’”’ procedure of the trial-and-error method. It 
gives results to the same order of accuracy as the original 
data without the use of an increased number of signifi- 
cant figures in the intervening processes and without the 
necessity of converting to a number times ten to a power 
for multiplying and reconverting to decimals for addi- 
tion. It allows solution for any desired root without 
previous solution for, and elimination of, any other 
roots. Direct application of the machines to the 
method is not practicable primarily because of the 
limitations of the machines in dividing. By the simple 
expedient of saving until the final operation each divi- 
sion involved in the successive reductions in order of the 
determinant, the machines are usable for the major por- 
tions of the computations. Preselection of the order of 
starring allows an IBM machine procedure that is a 
simple pattern of repetitive operations. 


STATEMENT OF PROBLEM 


The general flutter equations as expressed by Bleak- 
ney” ? are as follows: 


(An + Mu)ai + (Arz + Miz)a2.... + (Ain + Min)an 


= 0 
(An + Moja: + (Ase + Mo)d2.... + (Aon + Mon)dn 
= 0 
(An + My) + (A n2 + Myz)d2.... + (Aus + Man), 
= 0 
where 
A, = three-dimensional aerodynamic coefficient 
My = —my[1 — 2yX(1 + igy)] 
= three-dimensional mechanical coefficient 
a; = complex amplitude 
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For a steady-state solution, the determinant of the 
coefficients of the amplitudes must equal zero. 

The material furnished the tabulating department by 
the engineering flutter analyst includes the basic data 
pertaining to the problem and additional information 
that determines the solution desired. 

The basic data include: 

(1) A matrix of the complex aerodynamic co- 
efficients for each 1/k to be investigated. 

(2) One or more completely real inertia matrices. 

(3) One or more completely real stiffness matrices to 
correspond to the inertia matrices of (2). 

It is usually convenient to furnish two sets of stiffness 
and inertia matrices, one set for symmetric modes and a 
second set for antisymmetric modes. The additional 
information includes: 

(1) Designation of the specific aerodynamic matrix 
for which a solution is to be obtained. 

(2) Designation of the particular stiffness and in- 
ertia matrices. 

(3) A factor P.= po( ‘V1 — M?)/p which determines 
the altitude and Mach Number at which the solution is 
obtained. 

(4) Designation of the “solution spot,” generally 
one of the main diagonals. 

(5) An estimated value of X and appropriate values 
for the g;,’s. 


SOLUTION 


The first step in solution is that of forming the com- 
plex flutter determinant. This involves multiplication, 
in the IBM machines, of the basic inertia and stiffness 
matrices by the requisite P, PX, and Pg factors, adding 
the resulting stiffness and inertia matrices to the aero- 
dynamic matrix and shifting the rows and columns so 
that the solution spot appears in the determinant (11) 
spot. Following these operations there are mn? cards 
punched with complex numbers and one additional card 
(the purpose of which will be explained later) punched 
with a completely real number. 

Each of these cards is identified with its position in 
the determinant by the subscript numbers and contains 
the numerical values for use in the next step of the solu- 
tion. The form of the determinantal equation is 
now 


Vu+mQnZ Viv 


Va V2 


Vis. oe Vin 
7, 
a ae 


Vu Vio Vins. ‘ Vian | 
where V is a complex number, mQo, is completely real, 
and Z is the complex unknown X(1 + ig). 

The second step is that of reducing the order of the 
determinant. This is called a ‘‘starring’”’ process and is 
accomplished by forming the product Vij Van (where Viy 
refers to all terms not included in the mth row) and sub- 
tracting the products V;,,V,;. The resulting vector is 


Vis’ = VisVan — VinVas (2) 
and the resulting determinantal equation 


Vii + VinmMQZ Vis! Vis' . 2. Vie! 0 


Vor" Veo! Ves'...Vop1' 0 
Ves Var as. — Veu.t” 0 = 0 (3) 
Van Via Vas tae Vue-1 Vinw 


where the superscript J refers to the value resulting 
from the first starring and defines a determinant of one 
less order, namely, the minor of Vy. 

Similar starring operations are then made on each of 
the resulting reduced determinants, until, after 1 — 1 
starring operation, only terms in the (11) spot remain. 
The required condition that the determinant equal zero 
is that the sum of these terms be equal to zero; the re- 
sulting equation is 


ad ” + [( Van) ( Fine 1) ( Vaca wal eee ( Va Si aj x 
[mM%X(1 + ig)] =O (4) 


from which the value of Z = X(1 + ig) can be obtained. 
In order to have effected a solution, this calculated 
value must agree with the trial value that had been sub- 
stituted in the stiffness terms. 

The operations required to reduce the determinant 
by one order are accomplished in the IBM machines 
through the use of a “master deck”’ of [4(m — 1)? + 1] 
cards properly punched? to allow for the multiplication 
and collection of complex products as required by the 
starring process outlined above. Four detail cards are 
provided for each term of a determinant of one less 
order. The first card provides for the multiplication of 
real times real and real times imaginary; the second, for 
imaginary times real and imaginary times imaginary. 
The third and fourth cards provide for similar multipli- 
cations but with signs reversed. 

A ‘‘working deck’’ copy of the master deck is first 


. gang-punched. Numerical values from the cards of 


step one are then punched into the correct fields of the 
working deck cards. The multiplications are carried 
out and checked, and the values of the four detail cards 
are summarized to a single card for each spot in the 
reduced determinant. The single additional card previ- 
ously mentioned, designated as the (01) card, is proc- 
essed according to the same pattern. It is multiplied 
in each reduction by the starring term. 

One cycle of reduction is thus completed, and a new 
set of cards punched with numerical values has been 
obtained defining a determinant of one less order. This 
determinant is the minor of the V,,, term of Eq. (3) and 
must be equal to zero since V,, ¥ 0. 

For each order determinant a master deck is provided 
which determines the operations required to reduce the 
order by one. Cycles of operations are continued until 
the determinant is reduced to first order, the deck at 
that time consisting of two cards. The (01) card is the 





produc 
times ; 
(mo) 
the (01 
unknov 

A fi 
(11) ca 
results 
denom: 
solutio’ 
readily 


A s 
values 
practic 
three t 
value ¢ 
analyst 
Bleakn 
and in 
lated v 

In n 
value s 
trial is 
particu 
In som 
modes 
separat 
in the ] 
methoc 

It sh 
tion all 

Addi 
readily 
exampl 
tion of 
able w: 
some a 
ods are 
determ 
equipn 
large. 

One 
vestiga 
X vary 
determ 
of the « 
zero ar 
plane ; 
and da 
require 
ampliti 
above 
oscillat 
axis ca 
torque. 


(2) 


(3) 


ing 
ne 





SOLUTION OF THE FLUTTER DETERMINANT 173 


product of the progressively obtained starring terms 
times a real number which is the stiffness coefficient 
(m%,) for the solution spot. The (11) card divided by 
the (01) card determines the negative of the value of the 
unknown Z. : 

A final operation of multiplication of the (01) and 
(11) cards by the negative conjugate of the (01) card 
results in a complex numerator and a completely real 
denominator and completes the IBM operations for 
solution of the determinant. The value X(1 + ig) is 
readily obtained from the single manual division. 


DISCUSSION 


A solution is obtained only when the calculated 
values of X and g agree with the substituted values. In 
practice this requires an average of between two and 
three trials. The method of converging to the correct 
value depends on the preference and experience of the 
analyst in charge. A preferred method, developed by 
Bleakney, assumes linearity in the region of the solution 
and involves substitution of two assumed and two calcu- 
lated values as described in the Appendix. 

In many cases the first guess results in a calculated 
value so close to the correct solution that the second 
trial is merely a proof of the correctness. This occurs 
particularly when the modes are not strongly coupled. 
In some cases, when solution frequencies of two different 
modes are almost identical, it is extremely difficult to 
separate the roots. However, this difficulty is inherent 
in the problem and is probably experienced in any other 
method of solution. 

It should be noted that the proposed method of solu- 
tion allows for other than equal g solutions. 

Additional information of considerable value is 
readily obtainable with the IBM method, one particular 
example being the investigation of the effects of varia- 
tion of parameters. Amplitude ratios are also obtain- 
able with the same machine pattern of operations and 
some additional calculations. These specialized meth- 
ods are possible because of the ease with which the basic 
determinants can be formed and solved with IBM 
equipment when the number of solutions required is 
large. 

One example of a use of this nature occurs in the in- 
vestigation of control-surface-driven modes. Values of 
X varying from 0 to infinity are substituted in the basic 
determinant and the values of M = —m + mQ(X + ig) 
of the control surface for which the determinant equals 
zero are obtained. Values are plotted on the complex 
plane and represent graphically the stiffness, inertia, 
and damping properties of the control surface which are 
required to maintain harmonic oscillations at constant 
amplitude. (See Fig. 1.) Where the curve extends 
above the real axis, damping is required to prevent 
oscillations of increasing amplitude and the imaginary 
axis can be evaluated in in.Ibs. of required damping 
torque. Points below the axis determine the control 
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Curve of mQ vs. mQg. All other g’s = 0; 1/k = 0.5. 
Representative values of X noted on curve. 


Fie. 1. 


surface properties for harmonic oscillations of decreas- 
ing amplitude. It should be noted that each curve is 
for a particular set of structural damping terms for the 
other modes. Different curves at the same 1/k value 
can be obtained by varying the altitude factor P, or the 
inertia, stiffness, and/or structural damping parameters 
of the other modes and so evaluating the effect of these 
changes. 

Each set of IBM calculations is directly usable in this 
type of investigation. There is no trial-and-error solu- 
tion or iteration involved. 

Calculation of amplitude ratios is readily possible as a 
continuation of the process from which a root was deter- 
mined. The procedure of the starring process is 
identical with that used in the solution of simultaneous 
linear equations. That is, the array 


ile 0 ae 
hie Ww 0 ee 0 
Va" -3 Ves" -3 Va" -3 0 | = 0 (5) 


— , 0 

Vin Vas Ven ive Pa 

the values for which have already been obtained in the 
solution for X is equivalent to the set of linear equations 


Ves" ~*a, oa Ve" ~*a, -o Ods o> 0d, = 0 
Var a + Via" "s, + Var “te. ; 0d, = 0 
ist rey; cises ve 0a, =0 (6) 
Vin + Vinge + Vngls oe Vann = 0 


In these equations all amplitudes except a; and a, 
have been eliminated from the second row of the original 
equations by a systematic procedure easily recognizable 
as the usual method of solving simultaneous linear 
equations. It remains only to eliminate the unwanted 
amplitudes in the other rows. This is accomplished by 
a continuation of the starring process after dividing 
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through by the amplitude a;. The resulting equations 


are then 
tle oo alae) es =0 
ay 
Ty Y-3 49 4+ ee. 0 md \@ 
Me 0 =0 
N-ly.4Q0 +0 YTV = 0 


where Vx" ~? is the first starring term, and’ V3;,¥~? = 
Vo -3 Veg -2 — Vy"~2 Vyq' -3, 

From the resulting Eqs. (7) the amplitude ratios are 
obtained. The IBM procedures and operations re- 
quired to produce the above result are identical in 
principle with that of reducing the determinant. All 
that is necessary is an additional set of master decks 
defining the operations required to continue the process 
to its final form. 

The entire success of the method hinges on a card 
processing technique that allows computation and com- 
bination of the required numerical values with a min- 
imum number of machine operations and card manipu- 
lations. Development of the technique is due to 
William D. Bell of the Lockheed Tabulating Depart- 
ment. Two Multiplier Boards are permanently set up, 
one for multiplication, and one for proving. One addi- 
tional Reduction Summary Board is required. The 
machine operators need no special training. Modifica- 
tions have been found to reduce operations still further, 
but the resulting complications to the machine oper- 
ator’s procedure make them impractical. Additional 
improvements in technique may be expected with im- 
proved IBM equipment. 
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APPENDIX 
T; = assumed complex root substituted in first 
trial 
C, = calculated complex root resulting from sub- 


stitution of T; 


T, = assumed complex root substituted in second 
trial 
C, = calculated complex root resulting from sub- 
stitution of T; 
AC=Q-G 
AT = T2 — Ti 
Zo = probable root assuming linearity in region of 
solution 
AC 
Zo = C; ricer Ife me Zs 
0 2+ AT = ac $ 2 2) 
AT 
= T: momteompoonees! (Cy xe iT; 
2+ AT — rc | 2 2) 
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Research Studies Directed Toward the 
Development of Rational 
Vertical-Tail-Load Criteria 


LAWRENCE A. CLOUSING* 


National Advisory Committee for Aeronautics 


ABSTRACT 


During recent years the need for revision of methods of com- 
puting the maneuvering vertical-tail-design loading conditions 
has been demonstrated. The advance in airplane performance 
and the increase in the powerfulness of the controls have both 
caused an increased need for more rational vertical-tail-load 
criteria. 

The thought is expressed that vertical-tail-load criteria should 
be more closely related to the stability and control characteristics 
of the airplane. It is pointed out that, for the purpose of flight 
verification of methods of computing the loads, measurement of 
sideslip angle and rudder angle rather than vertical-tail load may 
provide sufficient information for a check of the essential parts 
of the methods of computing load. 

Three maneuvers are discussed as being possible to use as a 
basis for specifications for the computation of critical vertical- 
tail loads. These maneuvers are the rolling pull-out maneuver, 
the steady sideslip, and the fishtail maneuver. For each ma- 
neuver, information on methods of computation of load is pre- 
sented and a comparison is made with the results of flight tests. 
The possible relation of each maneuver to the critical load distri- 
bution is discussed. 


INTRODUCTION 


geet FAILURE of vertical-tail surfaces on air- 
planes in maneuvers during recent years has dem- 
onstrated the need for revision of methods of computing 
maneuvering limit-load values. There is, at least to 
the author’s knowledge, no recent history of failures of 
vertical-tail surfaces due to gust loads. It appears, 
therefore, that the gust loading conditions that are now 
considered in quite a rational manner are satisfactory, 
but that the rapid advance in airplane performance 
and the greater emphasis being placed on effective 
control with reasonable control forces even at high 
speed, particularly in military attack and fighting types 
of aircraft, have made the need for more rational ma- 
neuvering vertical-tail-load criteria necessary. Further 
evidence of the need for revision of vertical-tail-load 
criteria is the considerable difference between the 
Army, Navy, and Civil Aeronautics Administration 
specifications for maneuvering limit loads. 

It is the purpose of this paper to present information 
leading toward improved maneuvering vertical-tail- 
load criteria and not to present finished specifications. 

Presented at the Annual Summer Meeting, I.A.S., Los Angeles, 
July 18-19, 1946. 

* Head of Flight Research, Ames Aeronautical Laboratory. 





The nature of the load distribution, which is also a 
large factor in the design of the vertical tail, will be 
touched upon in this paper only to the extent to which 
the subject comes up in relation to the selection of criti- 
cal loading conditions on various parts of the vertical 
tail. 

The N.A.C.A. research work, which will be discussed, 
is related to the manner of specifying limit values of 
maneuvering vertical-tail load for structural design 
purposes. The research on this subject has been di- 
rected (1) toward selecting the types of maneuvers 
which may occur in flight and which, if considered, will 
completely establish the critical maneuvering design 
loads on the tail; (2) toward setting up methods based 
on the stability and control parameters of the airplane 
for the computation of vertical-tail loads; and (3) to- 
ward carrying out for each maneuver an experimental 
verification of the predicted loads by measurements in 
flight of load or important load parameters on various 
aircraft that differ as widely as possible with regard to 
pertinent stability parameters. If this program is car- 
ried to its logical conclusion, it will ensure adequate 
strength and, by reducing the factor of ignorance, tend 
to eliminate unnecessary weight. 


COEFFICIENTS AND SYMBOLS 


b = wing span, ft. 

S = wing area, sq.ft. 

S: = vertical-tail area, sq.ft. 

S,; = rudder area, sq.ft. 

m = mass of airplane, slugs 

W = weight of airplane, Ibs. 

I, = tail length, ft. 

Iz = moment of inertia about the Z-axis, slug-ft.? 
p = air density, slugs per cu.ft. 


V = velocity of airplane along flight path, ft. per sec. 
= free-stream dynamic pressure (1/29V?), Ibs. per sq.ft. 


qd 

qt = dynamic pressure at tail, lbs. per sq.ft. 

be = 2m/pSb 

n = load factor 

t = time, sec. 

B = angle of sideslip (positive when right wing is forward) 
p = rate of roll, rad. per sec. 

r = rate of yaw, rad. per sec. 

L, = vertical-tail load, lbs. 

L = moment about X-axis, ft.lbs. 

N:; = normal force on vertical tail, Ibs. 

Cy, = vertical-tail normal-force coefficient (N:/g¢S:) 
Ci = lift coefficient (nW/g5S) 

C; = rolling-moment coefficient (L/gq5Sb) 
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load using measured values of 8 and 6,. (From reference 3.) 


Cp average chord of rudder, ft. 


Ch rudder hinge-moment coefficient (H/qS,c,) 

AC; = rolling-moment coefficient due to lateral-control de- 
flection 

C, = yawing-moment ccefficient (N/gq5Sb) 

AC, = yawing-moment coefficient due to lateral-control de- 
flection 

Cg = OC,/0B 

Crip = OC,/28 

GQ, = 0C,/d(pb/2V) 

6, = rudder deflection, degrees 

btab = rudder-tab deflection, degrees 

K_ = ratio of rudder-pedal force to rudder hinge moment, 
lbs. per ft.Ibs. 

HH = rudder hinge-moment, ft.lbs. 

F, = rudder-pedal force, Ibs. 


COMPUTATION OF VERTICAL-TAIL LOADS 


Fundamental to the computation of vertical-tail 
loads is the equation shown in Fig. 1: 


OCy; oCy; 
n= 50[0(G)+0(%)] 
It is thought that, at present, literature on flapped 
airfoil characteristics is extensive enough to permit the 
constants OCy,/O08 and OCy,/06, to be estimated 
fairly accurately when the geometry of the tail and the 
end-plate effects of the horizontal tail in increasing the 
effective aspect ratio are taken into consideration. 
Although the validity of Eq. (1) has been demon- 
strated in numerous works, Fig. 1 gives further confir- 
mation to its validity. Values of vertical-tail load, as 
obtained in rolling pull-out maneuvers from pressure- 
distribution measurements in flight on a pursuit-type 
airplane, are compared with those computed by Eq. 
(1). It is emphasized that, in this computation, meas- 
ured values of sideslip angle and rudder deflection were 
used. The values of 0Cy,/08 and OCy,/06, were com- 
puted from the geometry of the vertical tail, consider- 
ing the end-plate effects of the horizontal tail in in- 
creasing the effective aspect ratio of the vertical tail. 
In general, in the application of Eq. (1) to the calcu- 
lation of limit loads, the speeds that will be significant 
with regard to critical loads will be so high that the 
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values of g, and g will be very nearly the same and that 
the use of g rather than g, is permissible. It is evident 
that the determination of critical loads on the vertical 
tail by use of the foregoing equation is largely one of 
detetmining critical combinations of sideslip angle 8, 
rudder angle 6,, and dynamic pressure g, which may be 
obtained in the maneuver. This is so not only because 
of the use of these factors in Eq. (1) but also because, 
in considering conditions for critical vertical-tail loads, 
it is necessary to consider the type of chordwise loading 
as affected by rudder position. For example, critical 
rudder load conditions may not necessarily occur during 
conditions that produce maximum total vertical-tail 
load or during conditions that produce maximum fin 
load. Therefore, this paper will be largely devoted to 
a consideration of methods of computing sideslip angle 
8, rudder angle 6,, and dynamic pressure g, in ma- 
neuvers that may lead to critical combinations of these 
parameters on either rudder, fin, or the entire vertical 
tail. 


THE ROLLING PuLL-OuT MANEUVER 


Because of the large yawing moments produced by 
abrupt deflection of the aileron in accelerated flight, 
especially on ‘the more maneuverable-type airplanes, 
the idea was put forth in a recent N.A.C.A. Bulletin 
by Gilruth' that a turn or a dive pull-out, in which 
ailerons were abruptly deflected, would be a maneuver 
in which critical vertical-tail loads might be encoun- 
tered. Accordingly, Gilruth made calculations of the 
approximate sideslip angles and vertical-tail loads 
that might be produced in this maneuver. For an air- 
plane with a rolling ability measured in terms of pb/2V 
for 50 Ibs. stick force at various indicated air speeds, 
as shown in Fig. 2, computations for a rudder-fixed 
condition resulted in sideslip angles of the order shown 








in Fig. 3. In Fig. 3 is also shown the equation Gilruth 
used for the computation of sideslip angle 8, 
B = (C,/8)(pb/2V)[1/(0C,/08)] (2) 
"ae { 
wy L ! | 





aes 
08 ad aaa 50 POUND STICK FORCE 




















AILERON -06;-———— we ae 
EFFECT- Pgh | 
IVENESS 94|__ | 
pb /2V | 
Bs | cee } 
° | 
100 200 300 400 500 








'NDICATED AIRSPEED, mph 


Fic. 2. Aileron effectiveness assumed in computation of 
vertical-tail load possible on fighter airplane in rolling pull-out. 
(From reference 1.) : 


The portion of the equation (C,/8)(pb/2V) is equiva- 
lent to the value of the yawing-moment coefficient 
C, because of the combined effects of ailercn deflection 
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Fic. 3. Angles of sideslip computed for fighter in rolling pull-out 
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Fic. 4. Vertical-tail loads computed for fighter in rolling pull-out 
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As may be seen from information presented in 
reference 2, it can be assumed, in arriving at this value, 
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that for elliptical span loading the yawing-moment 

coefficient due to aileron deflection and the yawing- 
moment coefficient due to rolling velocity are each 
equal to a value of (C,/16)(pb/2V) or one-half the 
total yawing-moment coefficient. 

It may be seen from Fig. 3 that at some intermediate 
speed a sideslip angle sufficient to stall the vertical tail 
of the airplane may be obtained in this maneuver. Fig. 
4 shows the results of computing the vertical-tail load 
by means of Eq. (1), using values of sideslip angle 8 
such as are presented in Fig. 3. It is apparent from 
Fig. 4 that vertical-tail loads may be obtained in abrupt 
aileron rolls at high accelerations, which can account 
for vertical-tail failures experienced on some airplanes 
in abrupt rolling maneuvers. It is apparent also from 
this figure that skillful use of the rudder on the part of 
the pilot in rolling maneuvers to minimize sideslip can 
greatly reduce the vertical-tail load. 

In a more recent N.A.C.A. Technical Note by White, 
Lomax, and Turner,’ the analysis of Gilruth has been 
compared with flight tests, and refinement of the meth- 
ods of computing sideslip angle in the rolling pull-out 
has been made. In the course of preparation of this 
Note, data obtained from flight tests of a single-engined 
pursuit-type airplane were analyzed as noted in Fig. 5. 
Sideslip angles are used in this figure rather than tail 
loads as the sideslip angles afford the same accuracy of 
comparison in the rolling pull-out maneuver. 

It may be seen that these tests showed computed 
values, using Eq. (2), which were in the neighborhood 
of 50 per cent of those measured. An explanation of 
this was found in examining the time histories of the 
rolling pull-out maneuver as recorded in flight on several 
airplanes. The examination revealed yawing accelera- 
tions of considerable magnitude at the time of maximum 
sideslip. By including the effects of yawing accelera- 
tion in the equations for computation of sideslip angle, 
the second equation of Fig. 6 was derived. 


_ (CG,/8)(pb/2V) + (Iz/qSb) (dr/dt) 
“: dC,/oB 


In this equation the second term of the numerator is 
the term that allows for the effect of yawing accelera- 
tion, and the first term corresponds to the allowance for 
adverse yawing moment previously presented in Eq. 
(2). For the pursuit airplane on which tests were 
made, this equation gave better agreement with results 
measured than did Eq. (2), as may be seen from Fig. 
6. 

Fig. 7 shows a comparison of computed values of 
sideslip with those measured in flight on a pursuit-type 
airplane. In this case the comparison is made for two 
configurations of the airplane. In the case where the 
airplane has a large vertical tail, better agreement is 
shown for Eq. (3) in which yawing acceleration is con- 
sidered than for Eq. (2), which underestimates the side- 
slip angle actually obtained. In the case of the air- 
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tween computed and measured values of sideslip angle obtained 
in rolling pull-outs using approximate equations. 
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plane with a small vertical tail, Eq. (2) gives better 
agreement with flight results than Eq. (3), although 
Eq. (2) still somewhat underestimates the sideslip 
angle, whereas Eq. (3) greatly overestimates it. 

The inability to obtain agreement between computed 
and measured values of sideslip in the rolling pull-out 
to a more accurate degree is considered in detail in refer- 
ence 3. In this Note, further refinements of the calcu- 
lation of sideslip angle in the rolling pull-out maneuver 
are made by setting up and solving the rather well- 
known complete equations of lateral motion of an air- 
plane, such as are described, for example, in reference 4. 

Fig. 8 shows one of six charts in reference 3 which 
present the results of a numerical solution of the com- 
plete equations of lateral motion over a considerable 
range of variables for the purpose of providing more ac- 
curate means for the prediction of sideslip angle in the 
rolling pull-out maneuver on airplanes of differing sta- 
bility parameters. Separate charts are drawn in refer- 
ence 3 for various values of the yawing-moment-coef- 
ficient change with change in angle of sideslip Cng and 
for two values of change in rolling-moment coefficient 
with sideslip Cig. As shown in Fig. 8, any one chart 
permits interpolation for the effects of yawing moment 
due to lateral control deflection AC,, the effects of the 
density factor we, and the effects of the relationship 
(AC,/Ci)(C,/Crg), a term that may be rewritten with 
little change in accuracy as (pb/2V)(C,/Cng). In the 
analysis that led to these charts it was shown that, 
within the limits assumed to occur in practice (of the 
variables aspect ratio, taper ratio, and moment of iner- 
tia about the various axes), the effect of these variables 
on the computation of the sideslip angle 8 was of second- 
ary importance. Consequently, the effect of variations 
of these quantities was not considered in the prepara- 
tion of the charts. 

In considering the numerical analysis over the range 
of conditions investigated, it may be seen in Fig. 9, 
as well as somewhat in Fig. 8, that the sideslip angle is 
approximately expressed as one-fourth the relationship 
(AC,/Ci,)(C,/Cnrg) plotted as the abscissa of the charts. 
In other words, 


B= 0.25( So) =) ~- CG po sian I a (4) 
Cw /\Cug 4 2V (0C,/08) 

Eq. (4) gives a value of sideslip angle in the rolling 
pull-out maneuver twice that given by Eq. (2) and 
very nearly that given in Eq. (3). 

Fig. 9 permits an evaluation of the amount of error 
which might result from use of the approximate solution. 
Lines are shown on this figure for the highest and lowest 
values of sideslip angle which would be obtained in the 
rolling pull-out maneuver for practical limits of vari 
ables. It is seen that, for airplanes of high values of 
loading mye, high change of yawing moment with side 
slip Cng, and large rolling-moment change with side 
slip Cig, the sideslip angle 8 is smaller than given by the 
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approximate solutions; whereas, when these values 
are small, 6 is larger than given by the approximate 
solution. 

It should be noted that the data of Fig. 9, as well as 
of Fig. 8, were computed considering Cxg constant over 
large angles of sideslip. Inasmuch as this may not al- 
ways be the case, the analysis was further extended, as 
shown in Fig. 10, to show the effect of nonlinear varia- 
tions of C, with 8 on the maximum sideslip angle. It is 
seen that a considerable variation in the computation of 
8 is possible because of nonlinearity. This chart has 
been used in reference 3 to develop methods for utilizing 
the linear charts previously discussed to compute £6 for 
nonlinear cases. 

It is realized that the evaluation of the variation of 
C, with 6 is difficult to establish accurately. It is be- 
lieved, however, that, if the values of this relationship 
were determined for a number of airplanes from flight 
tests and if these measurements were more completely 
compared with values as computed from the geometry 
of the airplane, the difficulty of accurately establishing 
this relationship by computation would be greatly less- 
ened. It might be mentioned here in this regard that 
a recent N.A.C.A. Technical Note by Bishop and 
Lomax® should be helpful in assisting the gathering of 
accurate flight data on the parameter Cng for various 
airplanes. 

In evaluating the place of the rolling pull-out ma- 
neuver in vertical-tail limit-load specifications, it should 


be noted that the maneuver gives loads only for a high- 


angle-of-attack, small-rudder-angle condition in which 
a gust-load-type distribution occurs which does not 
necessarily cause critical loads for either fin or rudder 
considered separately, although it might cause total 
load on the vertical tail of critical magnitude and there- 
fore critical fuselage torsional loads. 

In summarizing’research work done toward using the 
rolling pull-out maneuver for determining vertical- 
tail loads, it may be said that the computation of 8 by 
Eq. (4), in general, gives sufficiently accurate results, 
that the rudder angle may be regarded as zero in this 
maneuver, and that the tail load as computed using the 
proper value of 8 will agree well with values that actu- 
ally occur. It may also be said that the loads that may 
be computed thus may be larger and in closer agree- 
ment with actual maximum flight loads than those that 
would be computed by present specifications, particu- 
larly for airplanes with powerful lateral control and 
relatively low directional stability. 


THE STEADY SIDESLIP 


The steady sideslip maneuver frequently occurs in 
flight. It has, within the past several years, been 
brought into use as a basis for specifying limit-load condi- 
tions on the vertical tail in one set of vertical-design- 
loading specifications. In this set of specifications it is 


required that the vertical-tail strength shall be sufficient 
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to withstand 5° of yaw at limit diving speed with zero 
rudder deflection and with rudder deflected to maintain 
5° yaw. 

Objections have been raised to the specification on 
the basis that it leads to excessive tail strength. In 
general, however, measurements that were made have 
indicated that this requirement is not unduly conserva- 
tive, as may be seen from measurements of sideslip 
angle on a typical pursuit airplane in Fig. 11. The 
specification of limit load on the basis of a set angle of 
yaw at limit diving speed is objectionable, however, 
in that it encourages the use of small vertical tails 
which leads to directional-stability difficulties and, con- 
sequently, to the probability that sideslip angles greater 
than the amount considered in the computation of load 
may be encountered in actual flight. In other words, 
the specification is not so rational as it might be. 
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If the steady sideslip is to be used as a basis for a 
limit-load condition, it would seem to be more rational 
to define the limit load and the rudder angle at this 
limit load in terms of an increment of rudder force 
which a pilot might reasonably be expected to apply 
from a condition of trim at high speed. An equation 
for the computation of vertical-tail load under such 
conditions is readily derived from the following basic 
equations: 


C, = F,/KqS,c, 


a os) (=) (2) 
Cr a af dB + 6, 26, + Stab “eR 


and Eq. (1) previously presented in Fig. 1. 
tion for tail load is 


RG ie. Si oc] 7 
oil Ze Nactan| Sabu @) 


or, as shown in Fig. 11, ? 


AF, oe 
pas, fe foe bod 


(5) 


and 
(6) 


The equa- 








(8) 
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It would seem that computation of AL, could be 
fairly accurately accomplished from this formula, inas- 
much as 0C,/0Cy, may be quite accurately estimated 
from the large amount of literature available on the 
aerodynamic characteristics of flapped surfaces of dif- 
ferent types. On Fig. 11, for a pursuit-type airplane a 
computed increment of tail load due to 180 Ibs. rudder 
force is compared with the increment measured. The 
value of 0C,/OCy; used was estimated from published 
reports. In this comparison, values for substitution in 
the equation were K = 2.0, S, = 11.07, c, = 2, S; = 
19.01, and (0C,/OCy,) = —0.116. 

Reasonably good agreement is shown here for power- 
off conditions. In the few cases checked, agreement for 
power-on conditions has been good. Fig. 11 also shows 
that, in keeping with what is indicated by the preceding 
equation, the increment of vertical-tail load due to an 
increment of rudder force is about a constant value 
regardless of speed. This would be true, however, only 
at speeds above those at which slipstream effects are 
large and only at those angles of sideslip below stalling 
of the vertical surfaces or at speeds above those at 
which the rudder would be against the stop. 

Although the effect of offset of the fin on the vertical- 
tail load is not considered in this equation, it is believed 
that the introduction of this factor is an unnecessary 
complication, inasmuch as with proper directional trim 
at the various speeds, which is the usual condition main- 
tained by a pilot, the load on the vertical tail may be 
expected to be reasonably close to zero for zero rudder 
force. 

With regard to the value of rudder force which should 
be used as a basis for a specification of this type, ex- 
perience has shown that 180 Ibs. force is about as much 
as a pilot is likely to impose on a rudder and that this 
might be a proper value to use. It is also possible that 
a value to use would be the limit rudder-pedal force 
for which rudder control systems are designed. By 
C.A.A. rulings the value of F, would then be 200 Ibs. 

For the purpose of computing pressure distribution, 
expressions for the evaluation of sideslip angle and rud- 
der angle due to an increment of rudder force are neces- 
sary. An expression for the increment of sideslip angle 


1S 


AB = ALI,/Sbg(OC,/08) (9) 


which upon substituting for AL, becomes 


cia E33 Beare lance oe 


An expression of the increment of rudder angle as de- 
rived from Eqs. (5), (6), (8), and (10) is 
~ KqS,¢,(0C,/26;) Sb(0C,/OCy,) 


In evaluating the place of the steady sideslip ma- 
meuver in vertical-tail limit-load specifications, it should 








Aé, 


be noted that the maneuver gives loads with a balancing 
load type of distribution in which the fin is at large 
angles of attack with large opposite rudder, in which 
case critical fin loads might result with but moderate 
rudder loads. 

In summarizing study relative to the sideslip ma- 
neuver, it is seen that fairly simple formulas may be de- 
veloped to specify the tail-loading conditions in this 
maneuver as a function of rudder force and that the 
maneuver may be one that would cause critical fin loads, 
However, a fundamental criticism that may be di- 
rected against the steady sideslip as a basis for vertical- 
tail limit-load determination is that loads may be ex- 
ceeded for equal values of rudder force in other ma- 
neuvers. For example, it is apparent that higher loads 
would occur if the rudder were abruptly applied against 
the sideslip. It is also apparent that higher loads would 
occur in dynamic maneuvers resulting from abrupt dis- 
placement of the rudder or in a fishtail maneuver. If 
such dynamic conditions were used for determination of 
limit load, consideration of the steady sideslip condition 
would appear to be unnecessary unless means were de- 
vised or restrictions were set up preventing excessive 
motion in lateral maneuvers of this type. The loads 
that may be imposed during dynamic conditions 
such as exist in the fishtail maneuver and in vari- 
ous modifications of this maneuver will next be con- 
sidered. 


THE FISHTAIL MANEUVER 


During the course of measurements of vertical-tail 
loads of a pursuit airplane in various maneuvers the 
largest vertical-tail load actually measured, a load that 
exceeded the design load, was obtained in a fishtail 
maneuver in which the rudder :novement was of the 
same period as the natural period of the airplane in 
lateral and directional oscillations. However, in this 
maneuver neither large rudder forces nor rudder angles 
were present. Objections have been raised to using a 
fishtail maneuver as a basis for setting up a design spec- 
ification for the vertical tail because it is a maneuver 
that produces excessive values of sideslip, while it is 
not a tactical maneuver or a maneuver that would be 
performed in any ordinary type of flying or stunting. 
This maneuver or some modification of it, or both, ap- 
pear, however, to be the most rational means of setting 
up a specification for the vertical-tail load that may be 
obtained in maneuvers. 

A recent N.A.C.A. Technical Note by Boshar and 
Davis® has presented the results of analytical and ex- 
perimental studies of the fishtail maneuver and of one 
modification of it, the instantaneous deflection of the 
rudder. It is shown in this Note that these maneuvers 
may be sufficient to define critical load conditions for 
the whole tail, as well as for the fin and the rudder. 
Equations are presented in the Note for prediction of 
the loads. 
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Fig. 12 presents results of computations made in ref- 
ference 6 on the nature of the loads that might be ob- 
tained on the vertical tail of a large flying boat at vari- 
ous times following an abrupt rudder deflection at zero 
yaw. It is indicated that, at two points in an abrupt 
rudder-deflection maneuver, critical vertical-tail load- 
ing conditions may exist. At point (a) critical rudder 
loads may exist, and at point (0) critical fin loads may 
exist. 

It is apparent, as is shown in Fig. 12, that the motion 
following an instantaneous deflection of the rudder 
leads to a greater angle of sideslip and, consequently, a 
greater magnitude of load than would be obtained in 
the steady sideslip for the same value of rudder force or 
rudder angle. Unless consideration of the relative 
probability of occurrence of the two maneuvers leads 
to a different conclusion, it would appear that the 
steady sideslip should be ruled out as a maneuver for the 
consideration of critical vertical-tail loads. 

Fig. 13 presents the results of the study in reference 6 
dealing with the loads that might be obtained in the fish- 
tail maneuver. It is shown here that an oscillating 
rudder deflection of small amplitude causes large loads 
that are reached in a short time. When the period of 
rudder motion is the same as the natural period of the 
airplane, the rudder motion or the rudder forces neces- 
sary to maximize the load are moderate. In the ex- 
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ample shown, a rudder amplitude of but 9° would cause 
failure of the tail with this type of deflection. 

The initial cycle of a fishtail maneuver may be con- 
sidered as a rudder reversal, and its critical nature, if 
executed at the natural period of the airplane, is shown 
by the fact that 86 per cent of the load corresponding to 
final resonance may be attained with only one cycle 
of rudder motion. 

The load distribution at maximum angle of sideslip in 
this maneuver is for high angle of attack and small 
rudder deflection corresponding to the gust-load condi- 
tion or the rolling pull-out condition. 

In the employment of the maneuvers just presented 
as a basis for limit-load specifications, a need is seen to 
exist for specifications of what is required of the airplane 
with regard to yawing maneuverability, particularly at 
high speeds where the maneuvers discussed produce 
critical loads. As mentioned in the discussion of the 
steady sideslip, it is believed that a rational method of 
specifying vertical-tail limit loads in the fishtail ma- 
neuver or in the abrupt rudder deflection would be on the 
basis of an increment of rudder-pedal force, a value of 
180 Ibs. being suggested as reasonable for initial con- 
sideration. However, further study is needed before 
quantitative values of this type can be agreed upon. 

For example, with the maneuvering capabilities of 
typical airplanes of the present day, the application of 
180 Ibs. of rudder force in a fishtail maneuver at high 
speed would no doubt result in failure of the vertical 
tail. This is graphically illustrated in Fig. 14, which 





3000 DESIGN LIMIT LOAD 


LL 


© FROM FAIRED DATA 100 80 





2000} $e @*——_sa0 — 
— | RIGHT 


1000; 






TOTAL o PEDAL 
VERTICAL- ___ FORCE 


TAIL LOAD. ° a 


LB ; 
$1000) 7 
“80 Lert 


= ——., 
——_—_ = 


— 











-2000 i 198- 
AIRPLANE LIMIT 
act MACH NUMBER: M=0.8 ~ 








100 200 300 400 500 
INDICATED AIRSPEED, mph 


Fic. 14. Maximum vertical-tail load obtainable in fishtailing 
maneuvers as a function of air speed for several values of rudder- 
pedal force. 


presents the results of some measurements of vertical- 
tail load on a typical pursuit airplane. It is apparent 
from the extrapolated curves that, at high speed, 180 
Ibs. rudder force would cause the vertical tail to fail. 
The figure also indicates that greater loads for a given 
rudder force may be expected at higher speeds. It is 
believed that the trend is important to consider in 
setting up a design limit-load specification in view of the 
trend toward higher airplane speeds and in view of the 
greater difficulties with regard to dynamic stability at 
higher speeds. 
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In summarizing study on the fishtail maneuver and 
the abrupt rudder deflection, it may be stated that it 
appears possible that these two maneuvers may be 
sufficient to completely define the critical loading con- 
dition for the vertical tail and its component parts and 
that equations for computing the motion and vertical- 
tail load in the abrupt rudder-deflection maneuver have 
been shown to check experimental results reasonably 
well on the one airplane for which comparison has been 
made. No check of the equations for computing ver- 
tical-tail load in the fishtail maneuver has yet been 
published to the author’s knowledge. 
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CONCLUSIONS 


Figs. 15, 16, and 17 are presented to show the nature 
of the chordwise load distribution that is most apt to 
produce critical loads on the rudder, the fin, and the 
whole vertical tail. These figures list, in addition, the 
maneuvers that are most apt to lead to critical loads on 
the rudder, the fin, or the total tail. Mention is also 
made of the maneuvers that conceivably might cause 
critical loads. It may be seen from Fig. 15 that it is 
considered that critical rudder loads probably would 
result from the abrupt rudder-deflection maneuver at 
zero yaw. Fig. 16 indicates that critical fin load prob- 
ably would result from the latter part of the abrupt 
rudder-deflection maneuver, and Fig. 17 indicates that 
critical total load probably would result from a fishtail 
maneuver. 

The research work touched upon in this paper, leading 
to improved methods of determining limit vertical-tail 
loads, is continuing with the aim in mind of furnishing 
more experimental information ‘for aiding in increasing 
the accuracy of computation of loads, and for estab- 
lishing yawing maneuvering limitations, or for estab- 
lishing practical limits of rudder deflection or rudder 
force to consider in the specification of limit-load condi- 
tions. 
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Aerodynamic Center and Center of Pressure 
of an Airfoil at Supersonic Speeds 


H. W. SIBERT* 
Unwersity of Colorado 


SUMMARY 


Simple formalas (based upon the Ackeret theory for supersonic 
flow) are presented for obtaining the aerodynamic center and 
center of pressure of both symmetrical and cambered airfoils. 


INTRODUCTION 


i IS WELL KNOWN that according to the Ackeret (or 
linearized) theory for supersonic flow both the aero- 
dynamic center and the center of pressure of a sym- 
metrical airfoil are always at the midpoint of the chord. 
It appears, however, that no general formulas have as 
yet been given for determining the aerodynamic center 


and center of pressure of a cambered airfoil. Such 
formulas will be developed herein. 
SYMBOLS 
x axis = in direction of relative wind 
y axis = perpendicular to x axis 
z axis = along the chord 
t axis = perpendicular to z axis 
y’ = dy/dx 
y" = dt/dz 
a = angle of attack in radians of chord line with respect 
to the x axis 
c = chord length 
M, = free-stream Mach Number 
qs = free-stream dynamic pressure 
Ds = free-stream static pressure 
p = local static pressure 
= pressure coefficient = (p — p,)/ds 
M =VM?-1 
M./2 = moment per unit span about midpoint of chord (posi- 
tive when stalling) 
N = chord force per unit span 
Subscripts 
u = upper surface 
L = lower surface 


AERODYNAMIC CENTER AND CENTER OF PRESSURE 
It is well known that according to the Ackeret theory 
P = 2y'/# (1) 


for both surfaces, provided the positive y axis for the 
upper surface points up and the positive y axis for the 
lower surface points down. Since the Ackeret theory is 
valid only when P is small, y’ must not be large. Thus 
the airfoil must be thin and at a small angle of attack. 
Moreover, the leading and trailing edges must be sharp 
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because y’ becomes infinite at a rounded leading or 
trailing edge. 

Let the chord pass through the leading and trailing 
edges, and let the positive ¢ axis for the upper surface 
point up and the positive ¢ axis for the lower surface 
point down. Since both y’ and ?¢’ are tangents of small 
angles, they may be taken equal to the corresponding 
angles in radians. Thus ¢,’, t,’, y,’, and y,’ may be 
taken as the angles shown in Fig. 1. By virtue of the 
positive directions assumed for the two y and the two ¢ 


axes, these four angles are all positive in Fig. 1. Since 
a is also positive in this figure, it is evident that 
w =h' —a, wy =’ +a (2) 


Substituting these values of y,’ and y,’ in Eq. (1), 
Py, = 2(t,’ — a)/M, Py = At,’ + a)/# (3) 
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Fig. 2 shows the pressures acting on two small seg- 
ments.of the upper and lower surface whose two ends 
are at the same distances from the ¢ axes. Since dz is 
positive for these two surfaces, it follows from Fig. 2 
and Eq. (3) that 


dN = p,ds, cos 6, — p,ds, cos 4, prdz — p,dz 
dN = [(pr — Ps) — (bu — bs)\dz = QAP, — Py)dz 
dN = q,(2/44) (2a + t,’ — ty')dz 
Since the airfoil is thin, the moment of the chord 
forces can be neglected. Hence, 
Me = — Slj2 dN = qs(2/¥) S72 (tu — tr’) dz 
(4) 


ll 
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because the integral of 2az dz between the limits z = 
—c/2 and z = c/2 is zero. 
Note that for either surface 


"tg dz = Sl, (dt/dz) dz = Sia iq zt 
Since t is zero at both z = —c/2 and z = c/2, it follows 
from an integration by parts that 
Sjot'2 da = [tL — 
Note that the last integral is the area S between the 
chord and the surface in question. Hence, 


Sipe ty't de = —Sy Sip tr'2 dz = —S, (5) 
Hence, from Eqs. (4) and (5) 
Me sa (29,/A4) (Sz — S,) (6) 


lin tdz = — pas B tdz 


Since the right-hand side of Eq. (6) is independent of 
the angle of attack, the aerodynamic center of any air- 
foil satisfying Ackeret’s assumptions is at the midpoint 
of the chord, and the moment about its aerodynamic 
center is the value of M,,, given in Eq. (6). 
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From Eq. (6) the moment about the midpoint of the 
chord is zero when S, = S,;. Since the center of pres- 
sure is located at the point about which the moment is 
zero, it follows that the center of pressure of an airfoil 
is at the midpoint of the chord when the area bounded 
by the chord and the upper surface equals the area 
bounded by the chord and the lower surface. Note 
that a symmetrical airfoil is a special case of the more 
general condition that S, = S,. If S, does not equal 
Sz, the distance of the center of pressure from the mid- 
point of the chord can be found from Eq. (6) and the 
well-known formula for obtaining the moment coef- 
ficient at a given distance from a point at which the 
moment coefficient is known. 


EXAMPLE 


Since there may be some doubt in the reader’s mind 
as to the validity of Eq. (5), the integral of ¢’z dz will 
be found for a surface composed of the three straight 
lines: (1) a line with ¢ increasing from zero at z = —c/2 
tohatz = —a’, (2) aline witht = 4 from z = —ato 
z = 6b; (3) a line with ¢ decreasing from h at z = b 
to zero at z = c/2 (see Fig. 3). The area between the 
chord and this surface is 


S=(¢c+a+td)h/2 


while 


_ h fo 
, ie ee dz — 
Pee = OE I 


h c/2 
aso zdz = —(c+a+b)h/2 = -S 
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Simple Analytical Equations for the Velocity 
of an Airplane in Unaccelerated Level, 


Climbing, and Diving Flight 


HUGH B. FREEMAN* 
Consolidated Vultee Aircraft Corporation 


SUMMARY 


It is shown that the velocity of an airplane in unaccelerated 
flight may be expressed as a simple function of the thrust and two 
constants of the airplane; the parasite area (f) and the thrust 
corresponding to the configuration for maximum (L/D). 

For level flight the equation for the velocity is: 


y = 14.02 To\? |” 
use \: E y'- i - (7) | (m.p.h.) (9) 
where 


T = effective thrust 


Ty = 1.128W/b +/f/e = thrust for configuration for maximum 
(L/D) = constant 
f = equivalent parasite area 


A simple method of solving Eq. (9) is presented. The con- 
ventional thrust horsepower required curves are obtained follow- 
ing the solution of this equation by the simple process of multi- 
plication. 

Similar equations are derived for the conditions of unacceler- 
ated climb and diving flight. 


INTRODUCTION 


Pp sinmpen METHODS of analysis are always tedious, 
their accuracy is limited, and they often obscure 
the essential physical relationships involved in a prob- 
lem. These difficulties very early became apparent to 
aerodynamicists working with performance problems, 
and many attempts were accordingly made to develop 
analytical methods to replace the graphical ones avail- 
able.” F 

The difficulty of expressing the performance of an 
airplane analytically stems partly from the complexity 
of the equations resulting when power is introduced into 
the analysis prematurely. If the thrust is used instead 
of power, extremely simple analytical equations result. 
From these equations the power relationships follow 
quite simply. The equations in the following pages 
were developed in an attempt to place the performance 
analysis of the airplane on a simpler and more analytical 
basis. While the author has not succeeded in this 
objective, it is believed that these equations are of 
sufficient interest and use to warrant their publication. 
This is especially true for performance computations 
using power plants, where the propulsive force is inde- 
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pendent of speed, as in rocket propulsion and in some 
forms of jet units. 

It is hoped that someone will investigate the possi- 
bility of extending this analysis with the object of ob- 
taining a more simplified method for the overall per- 
formance analysis. 


DERIVATION OF EQUATIONS 


Equation for Level Flight 


Assume with Millikan and Oswald that the polar of 
an airplane may be represented by an equation of the 
form 


Cy = Cop + Cr?/weAR (1) 
where 


Cp, = effective parasite drag coefficient (Millikan) 
e = airplane efficiency factor (Oswald) 
AR = aspect ratio = b*/S 


For unaccelerated level flight the thrust is equal to the 
drag: 


ll 


T = D = (Cp, + C,?/meAR)(p/2) V?S 
or 


(p/2) V? 
W/S 





T/W = (Cop + C,?/reA R) 
Writingt W/S = (p/2)V,", where (V;) is the velocity 


corresponding to C; = 1.0, and writing C, = 1/(V/Vi)’, 
the thrust becomes: 


T/W = Cp,(V/Vi)? + Way 


(V/V) a 
Transforming and multiplying through by (V/V)? 
Cp, (V/Vi1)* — (T/W)(V/Vi)? + 1/reAR = (3) 


This is a quadratic equation, the solution of which is; 


(T/W) + V (T/W)? — (4Cp,/weAR) 


2Co, 








(V/V;)? = (4) 


t In the original derivation the author made the substitution 
W/S = Cimaz.(e/2)Vs? and Cr = Cimar./(V/Vs)? in place of 
those shown above. The simplification of using W/S = (p/2) Vi?, 
was pointed out to the author by C. G. Wolcott of the Con- 
solidated Vultee Engineering Department. 
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This is a nondimensional equation for the velocity of 
an airplane in steady level flight, independent of alti- 
tude. Eq. (4) may be further simplified, however, 
after a consideration of the speed for level flight at the 
maximum (L/D) configuration. 

Returning to Eq. (2), differentiating with respect to 
the velocity and equating to zero: 


d T 1 2(V/V, 
: any (; ) = 2Cp,(V/Vi) — RM aS = 0 
d(V/Vi) \W mweAR (V/V;)4 
or 
(V/V;1)ot = 1/Cp,reAR (5) 
This is the velocity corresponding to the minimum 
thrust and to the configuration for maximum (L/D). 


Substituting this value of (V/Vi)o into Eq. (2) yields 
the minimum thrust: 


(T/W)o = 2Cp,(V/Vi)o? (6) 


Remembering that 


Vi? = (842/c)(W/S) [where 842 = (2/po)] 
Co, = f/S 
AR = }?/S 


Substituting these values into Eqs. (5) and (6) and 

transforming, we obtain 

Vo = (14.88/0')(W/b)'"(1/fe)'“* (m.p.h.) (7)* 

Ty = 1.128(W/b)(f/e)'” (8) 

T> is the thrust or drag of the airplane in the configura- 

tion for maximum (L/D) and is a constant for a given 
airplane. 

Returning to the general equation for level flight, 
Eq. (4), and introducing equivalent values for (V1), 
(Cp,), and (AR): 

842 (W\ TT, HES oe 
——{—]|7/W «= iW} = —— 
y2 = "e AS dL Asan joshi 7 
2f/S 


or 


y= hip VP 127K /UTe) 
"OR. 2 ra 2 1 979(W/h\2 F/ 
aa 272(W/6)*(f/e)| 
from Eq. (8) 
1.272(W/b)*( f/e) = To? 


substituting 
21 eeneera 
7 Gnckane (T ot VT? — T>?) 
of 


or, rearranging and changing to engineering units, 


, 14.02 ,/T _ (T\? |" , 
a V; : « y is (| (mph) (9) 








* This equation may be derived directly from Oswald’s equa- 
N.A.C.A. Report No. 408. 
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This is the equation for the velocity of an airplane in 
steady level flight in terms of the thrust (7) and all of 
the significant airplane parameters that are contained 
in the two airplane constants: f = CS and 7) = 
1.128(W/b) ~/f/e. 

Eq. (9) should be extremely useful in the performance 
calculations of jet airplanes where the thrust is moré 
convenient to use than the thrust horsepower. Com- 
pressibility corrections may also be easily incorporated 
so long as the initial assumption in regard to Eq. (1) 
remains valid. 

A later paragraph will present a simple and rapid 
method for the solution of Eq. (9). 


Equations for Climb 


"7 6 «= angle of climb 


The equations for climbing flight are derived in a 
manner similar to those for level flight. For un- 
accelerated climbing flight 


mo . we 2) ") 72 
T= D+ Wsind = (6, + te (5 V2S + 
W sin 6 


dividing through by W, substituting as before for W/S 
and eliminating C, by substituting its equivalent, which 
in this case is 


C, = cos 6/{(V/V;)"] 


we obtain 


+ sin6@ (10) 





cos” 6 
T/W = Cp,(V/V1)? + —— 
/ Dp / 1) + (V/ V;)°*reAR 


Transforrhing and multiplying by (V/V1)? we obtain 
the quadratic 


P : ‘ ae 

Cp, ( V/ Vi) - (T/W — sin 0)(V/ J Fi “4 an =' () (11) 
Solving for the velocity 
(V/Vi)? = q : po vee 
(T/W in 6) r/ 4Cp, cos? 6 

/W — 38 + STU => MRS a ee 
ancestry nat sigs Tse eo gla 

2Cp, 
(12) 


This is a nondimensional equation for climbing flight 
corresponding to Eq. (4) for level flight. 

Returning to Eq. (10) and differentiating with re- 
spect to the velocity, for a constant angle of climb: 


d 2(V/V;) cos? 6 _ 


eee es oh Fre = 
d(V/V;) (=) Co,(V/Vi) (V/V;)4xeAR 
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Equations for Unaccelerated Descending Flight—Terminal 


























- solving for the velocity Velocity 
7 4 = 2 
1 of (V/Vi)e* = cos? 0/Cp, reAR (13) The equations for unaccelerated descending flight or 
ned The subscript (C) denotes conditions corresponding to the terminal velocity at any angle of dive (6) may be 
= the maximum (L/D) configuration in climb analogous obtained simply by changing the sign preceding (sin @) 
to the subscript (0) for level flight. This velocity in all the equations for climb. The general equation 
nce corresponds to the minimum thrust. Substituting for descending flight corresponding to Eq. (19) is 
oré back into Eq. (10) the minimum thrust for climb at the 14.02 T+ Wsinod 
ym- angle (6) is oii, SRA ae 
e 4 , 
ted (T/W)e¢ = 2(Cp,/meAR)'* cos8+ sin (14) : 
(1) ? eile Ta ue W sin 0\2 (20) 
Transforming and substituting equivalent values for \ T+ Wsine 
pid (Cp,) and (AR) into Eqs. (13) and (14), the minimum a ; 
, , where 7, = minimum thrust or drag corresponding to 
thrust and the corresponding velocity are d P ; Sq? , 
/ the maximum (L/D) configuration for diving flight. 
To = 1.128(W/b)(f/e)'* cos8@+ Wsin6@ (15) Remembering that 7; = 7) cos @ — W sin @, this equa- 
1 tion reduces to 
14. - hk 
Ve = (W/b)'” (=) (cos 6)" —* (16) — —— 
fe . 2 fT W sin 0 
The relationships between the thrust and velocity for a f ak SMG oJtR GRP 
level flight and climb at maximum (L/D) are (compar- T,cos@ \?|” 
ing these equations with Eqs. (7) and (8)): l=qi- (2 W sin 6 (m.p-h.) 
| a To = Te cos 6 +)/W sin 6 (17) The terminal velocity in a vertical dive then, where 
acs = Vo Vcos 6 (18) 6 = 90°, sin @ = 1.0, and cos @ = 0, reduces to 
Returning to Eq. (12), substituting for (Vi), (Cp,), Vi = (14.02/0°7*) V(T + W)/f V2 
and (AR) their equivalent values, the velocity in climb ai 
becomes ioe emadanrare’ sll: 
i yo = (19.80/0"") V(T + W)/f (21) 
/S 421 : — *4.979 2(W/b b)2 cos? t 8 or, when the thrust is zero the terminal velocity 
ch — (T — Wsin @) l--—— —— vol i 
f (7. — W ein 0) V, = 19.80 VW/of (m.p-h.) (22 
From Eq. (15) it is seen that 
4 272(W/b)? cos? 0 = (To — W sin 6)? Solution of Equation for Level Flight 
aN l/s 
9 \2 
; then i hee y) E oe \t _ (7) (m.p.h.) (9) 
y — 14.02 T — Wsiné ‘ f 
: gt f A simple and rapid solution to this equation is 
= a gig a eee obtained when it is noted that the term within the 
hs + y:- = ( ae —wa at | (m.p.h.) (19) bracket may be solved once and for all. 
— oy Let T/T = A and (1 + V1 — A*)? = B. Note 
1) This is the equation for steady climbing flight corre- also that the ratio (7)/T) has a maximum value of 
sponding to that of Eq. (9) for level flight. When unity. By assuming a series of values for (A) we may 
(9 = 0) this equation reduces to Eq. (9). solve for the corresponding values of (B). These 
a TABLE 1 ins cot 
aa _ Tabulation Chart for Rapid Solution of Eq. (9). 
= = Vi e/T A | = [1+ Vi - G/T" 
T/T = . = 0.6 0.7 0.$ 10 0.9 08 O07 O06 05 04 0.3 
») ‘ = 0.4475 0.535 5 aa 0.751 1.0 1.199 1.265 1.310 1.342 1.366 1.385 1. - 398 
(1) Zz = To ‘A = | 
it (2) T/f = 
(3) V = 14.02/o'/ V/T/f-B = 
p- Speeds below that pan es _ Speeds above that 
! for maximum L/D | for maximum L/D 
(4) T.-hp. = TV/375 = 
Note: Jo = 1.128W/b Vf/e = (a constant for a given airplane). { = CppS = (a constant for a given airplane). 
) Eq. (9) is valid only over that portion of the airplane polar for which the initial assumption in regard to Eq. (1) is valid. 
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Illustrating solution to Eq. (9) and the significance of the parameter (7) in the cases of thrust required curves and 


the thrust horsepower required curves. 


values may then be placed permanently on a working 
chart, as illustrated in Table 1. The subsequent steps 
of the computation are outlined on the chart and are 
self-explanatory. 

Eq. (9) is plotted in Fig. 1 for two altitudes. All of 
the thrust required curves for different altitudes are 
tangent to the line (J) = constant). 

If the conventional thrust horsepower curves are 
desired, only one additional step in the calculations is 
required as indicated by step (4) in Table 1. In this 
case, all of the power required curves are tangent to the 
straight line whose slope is equal to (Zo) (Fig. 2). 
This line represents the overall maximum speed an 
airplane can attain in level flight for a given thrust 
horsepower (below speeds at which compressibility 
occurs). In other words, the maximum speed at which 
an airplane can fly, for given thrust horsepower, is that 
corresponding to the maximum (L/D) configuration or 
(Vo) as given by Eq. (7). In practice this speed can 
only be realized when the cabin is pressurized or other 
provisions allow the pilot to fly at the extreme altitude 
corresponding to this value of (Vp). 


Solution of Equation for Climb 


y = 1402 [7 — Wind 
~ ae f 


7 ° '/s 
To — Wsin oy 
1 = 4/1 — (————_—_ m.p.h.) (19 
| y (await apn.) co 
This equation may be solved in exactly the same 
manner as Eq. (9), letting 


Te - W sin 6 _ 
T—Wsind 


: '/3 
To — Wsin 2) 
lw alps (LL en 
| \ ( T — Wsin@ | - 
where (A) and (B) have the same numerical values as 
in Table 1, 





and 





However, once the level flight characteristics are 
known from the solution of Eq. (9), a much simpler 
solution for the climb equation is available. 





Let 
14. = 
Vetimb = - \" asm B (23) 
a f 
and 
Vievet = (14.02/0') »/T/f B (24) 


Assuming that (B) has the same value in both equations 
and dividing 





Vatimb we T’ — Wsin@ (25) 
Vievet - 


Let (7’) represent the thrust in climb, temporarily, to 
avoid confusing it with (7), the thrust for level flight. 
Now 








To — W sin 6 at To 
T’ — W sin 0 T 
hence 
Vesimay _ yi" — Wsin 0)/A 
Vievel T/A 
and since 


Te = Ty cos8+ Wsin@ (from Eq. (17)) 
then 
Vievel V cos — 

In other words, for corresponding values of (A) and 
(B), the velocity in climb is equal, for all practical 
purposes, to the velocity in level flight since the angle 
(6) is always small so that the square root of (cos @) will 
approximate unity. 

Hence, once the level flight characteristics are known 
it is only necessary, in order to obtain the climb charac- 
teristics, to compute the thrust (required) for climb for 
corresponding values of (A) from 


(26) 


Vetimb = 
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T’ = (Ty) cos 6/A) + Wsin @ (27) 


and plot this thrust against velocity, for corresponding 
values of (A) and (B), from Table 1. Or, if the thrust 
horsepower required for climb is desired, multiply the 
values of thrust from Eq. (27) by the corresponding 
values of the velocity as follows: 


T.hp. - T’ Vieves/375 





(28) 


Curves of (T.hp.) required may then be plotted against 
velocity in a manner similar to that for level flight. 


CoNCLUSION* 


It will be noted from the foregoing that the equations 
presented herein may be used for general performance 
computations only in cases where the thrust is inde- 
pendent of speed. For conventional propulsive means 
the equations are useful in determining the thrust or 
power required relationships, in the comparison of the 
relative merits of various designs, and allow the effects 
of changing the individual airplane parameters to be 
easily and rapidly determined in a more rational manner 
than by conventional, graphical methods. 


* After submitting this paper for publication, the author’s 
attention was called to a report by Lippisch, ‘Performance 
Theory of Airplanes with Jet Propulsion,” Translation Report 
No, F-TS-685-RE, October, 1946, Headquarters Air Matériel 
Command, Wright Field, Dayton, Ohio. Some of the equations 
developed by Lippisch appear to be almost identical with those 
presented here. His treatment is much more comprehensive in 
scope, however, than that of the present paper. 
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APPENDIX I 


A performance equation, more general in form than 
Eq. (3) for the level flight condition, may be derived by 
introducing power and propeller efficiency into the 
initial equations given under the heading ‘Derivation of 





Equations.” The resulting equation is 
5507 
Cop(V/Vi)! = Far V/Vs) + 1/neAR = 0 = F 


where 7 = propeller efficiency and P = brake horse- 
power. 

The general solution to this equation depends upon 
evaluating the propeller efficiency (7) in terms of the 
velocity. If the efficiency is considered constant, the 
equation may be solved simply by graphical methods. 
A similar equation may be set up for the climbing con- 
figuration. 


APPENDIX II 


Velocity in a Steady Horizontal Turn 


The equation for the velocity in a steady horizontal 
turn may be derived in a manner similar to that for 
level flight. The equation is 


a 1/3 
re ile= OTT 


where (¢) is the angle of bank. 
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Wind-Tunnel Turbulence Effects 


DIMITRI P. RIABOUCHINSKY 
French National Center of Scientific Research 


Epiror1AL NotE—The following historical account of early work on wind-tunnel turbulence effects is published at 


the request of the author. 


force on an airfoil are found under the conditions of highest turbulence. 


It should be noted, however, that according to modern experience the highest lift and normal 


Unfortunately, in the pioneering work of the 


author, no method was available either for quantitative measurement of turbulence or even for determining whether such 
a change as removing an elbow from a Gottingen return-flow wind tunnel increases or decreases the turbulence. 


INTRODUCTION 


| ie HIS INTERESTING Wright Brothers Lecture (Decem- 
ber 17, 1938), Dryden! mentioned that the discovery 
of wind-tunnel turbulence effects began with spheres 
and that a similar influence of turbulence on maximum 
lift coefficients was found subsequently. The author, 
who, 34 years ago, took an active part in the elucida- 
tion of this problem, would like to recall certain re- 
searches that seem to be completely forgotten, in spite 
of the influence they exerted and the historical interest 
they present. 


RESEARCHES 


Early in May, 1905, Rateau? in France and Mallock? 
in England established that the forces and moments 
that act upon a rectangular plate, measured in an arti- 
ficial wind, may have two different values when the 
angle of incidence is between 29° and 36° (Rateau) 
and about 40° (Mallock). 

Repeating this experiment in his wind tunnel at the 
Aerodynamic Institute of Koutchino, near Moscow, 
the author did not find‘ the discontinuity observed by 
Rateau and Mallock. Rateau® ascribes this discrep- 
ancy to the damper used in those experiments, which 
he considered masked the phenomenon. However, 
even without the damper, the author® never observed 
in the wind tunnel of the Koutchino Institute two dif- 
ferent positions where the plate remains in equilibrium 
during ‘‘several seconds.’ The position appeared more 
or less clearly, but it was always single. 

In 1910, Eiffel? and Prandtl* noticed a considerable 
increase of the resistance of a square plate in the vicin- 
ity of the critical angle of incidence corresponding to the 
change of the eddies’ configurations, but, while Eiffel 
presents the results obtained by a continual curve, 
Prandtl observed between 38° and 40°, the indetermi- 
nate zone of Rateau and Mallock. 


EXPERIMENTS 


Not having succeeded in finding in the wind tunnel of 
Koutchino either such an important increase of pres- 
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sure or such discontinuities, the author decided, in 
order to elucidate the causes of such discrepancies, to 
study carefully the influence of the conditions under 
which the experiments were carried out, which were 
considered a priori as secondary. In view of the said 
researches, he constructed, first, a wind tunnel with free 
air stream; second, a small Eiffel tunnel; and, third, 
a Prandtl-type return-flow tunnel, which, at that time, 
was not of open working section. 

These researches enabled the author to throw light 
upon the influence exerted upon the results of observa- 
tions by: (1) certain peculiarities of different types of 
wind tunnels; (2) the value of the ratio d:D, where d 
is a linear dimension of the body placed in the stream 
and D that of the stream section; (3) the means of sup- 
porting the body in the stream; and (4) the degree of 
turbulence of the latter. 

The chief results of comparative researches carried on 
by the author in the tunnel of Koutchino and in that 
of the Géttingen type are the following:® (1) By utiliz- 
ing the three methods for supporting a square plate in 
a stream, shown in Fig. 1 (a, 6, and c), the author ob- 
tained in the Koutchino tunnel the three curves (Fig. 2) 
of the orthogonal resistance coefficients of a square plate 
as function of the angle of incidence. Those three 
curves are continuous, and their maxima in the vicinity 
of 35° are comparatively weak. 

(2) In the Géttingen-type tunnel the author ob- 
tained, in applying the supports a and c of Fig. 1, the 
corresponding curves of Fig. 3 with strong maxima and 
discontinuous change of eddy configuration at an angle 
of incidence of about 40°. 
those of Géttingen, but the author also established 
that with the support 0} (Fig. 1) the strong maximum 
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Fic. 2. Influence of turbulence and of supports on orthogonal 
force coefficients of a square plate in the Koutchino wind tunnel. 
The curves a, b, and c were obtained, respectively, with the 
supports a, b, and c of Fig. 1. 
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and the discontinuity disappear (curve b of Fig. 3). 
This was a new and important result. This experi- 
ment put forward clearly the fact that a small pertur- 
bation of the flow in vicinity of the body can cause a 
considerable decrease of the force. 

(3) Noticing that the oscillations of the obstacle 
during measurement were somewhat stronger in the 
Koutchino tunnel than in the Géttingen return-flow 
tunnel, the author came to the conclusion that the 
stream fluctuations were somewhat smaller in the 
Géttingen-type tunnel, and he ascribed the discrepancy 
between the a and b curves in Figs. 2 and 3 to this dif- 
ference of the degree of turbulence in these two tunnels. 
Thus, the stalling phenomenon appeared as sensitive 
to turbulence. 

(4) It follows from remarks 2 and 3 that a local 
perturbation of the stream on the surface of the ob- 
stacle acts, in certain cases, similar to an increase of 
the stream turbulence. 

(5) If one intentionally increases the turbulence of 
the stream, for instance in removing an elbow of the 
Gottingen return-flow tunnel, so that the air is drawn 
directly into the tunnel, the strong maximum and the 
discontinuity of the resistance curve disappear. The 
curve d of Fig. 3 has been obtained under these condi- 
tions by mounting the square plate on the support a 
of Fig. 1. 

The author concluded his work by a query sign, thus 
inviting his colleagues to express their opinion upon 
his explanation of the causes which led to observations 
peculiarly contradictory in different aerodynamic 
laboratories. 
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Fic. 3. Influence of turbulence and of supports on orthogonal 
force coefficients of a square plate in a Géttingen-type tunnel. 
The curves a, b, and c were obtained, respectively, with the sup- 
ports a, b, and c of Fig. 1. The curve d, obtained with the support 
a of Fig. 1, corresponds to the case when an elbow of a Gottingen 
return-flow tunnel was removed, in view of increasing the turbu- 
lence of the stream. 


LATER RESEARCH 


None of the papers which the author has just quoted 
and which appeared in the course of the four years 
1909-1912, were mentioned in the historical sketch on 
the wind-tunnel turbulence effects given by Dryden.' 
He mentions only the following three publications issued 
in that period of time: (1) the value 0.18 of the drag 
coefficient of a sphere published in 1911 by Eiffel;'® 
(2) the statement made in 1912 by Foéppel,"' of the 
Aerodynamic Institute of Géttingen, that Eiffel’s 
published value 0.18 was obviously an error, probably 
a misprint, and that the true value was 0.44, or nearly 
three times as great; (3) a Note of Eiffel'® presented to 
the French Academy of Science, December 30, 1912. 
In this Note three curves connecting drag coefficients 
of spheres and speed are shown, and Eiffel adds that 
he was not mistaken in his previous paper,'® either in 
the measurement or in the calculation of the coeffi- 
cients in giving the value 0.18 and that, if the Gottingen 
Laboratory had used sufficiently great speed, it would 
have equally found this value, the coefficient 0.44 cor- 
responding to comparatively small velocities. No al- 
lusion to the influence of the stream irregularities is 
made in these three papers. 

In 1913 in the article ‘‘Fliissigkeitsbewegung,’’ in the 
first edition of the Handwérterbuch der Naturwissen- 
schaften, which appeared also as a separate booklet 
under the title Abriss der Lehre von der Flissigkeits- 
und Gasbewegung, Prandtl,'* speaking of the strong 
maximum and the discontinuity of the resistance curve 
of a thin plate, writes: ‘In a turbulent flow the phe- 
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nomenon disappears (Riabouchinsky),” thus confirm- 
ing the sound basis of the author’s explanation. It was 
natural to extend the author’s explanation to other 
contradictory results obtained in wind tunnels, and 
such an extension to the problem of the resistance of 
spheres was published the next year (in 1914) by 
Prandtl"* and Wieselsberger™ in Germany, and by the 
author himself'* in Russia. In the fifth Bulletin of the 
Aerodynamic Institute of Koutchino, published in 1914, 
the author'® assembled the researches he carried on 
during the years 1912-1913 on frictional resistance of 
rotating discs of which the surface was smooth or cov- 
ered by linen, on skin friction and momentum, on 
turbulent flow around a rotating cylinder, and on re- 
sistance of spheres and ellipsoids in the Koutchino wind 
tunnel. 

He has also given, in the same Bulletin, an 
illustration in which he represented all the torque 
coefficients of spheres, in air and water, with which he 
was acquainted at that time, beginning with those of 
Mariotte and Newton to end with those of Costanzi,!” 
Hesselberg and Birkeland,® Eiffel,!2 Maurain, and 
his own,’ as function of the Reynold’s Number. The 
curves published by Wieselsberger of the Aerodynamic 
Institute of Géttingen are not on this illustration, for 
the fifth Bulletin of the Aerodynamic Institute of Kout- 
chino was ready for print when the author received the 
Zeitschrift fir Flugtechnik und Motorluftschiffahrt of 
May 16, 1914, including Wieselsberger’s paper, but the 
author summarized the leading points of this paper in a 
footnote. 

Prandtl’s experiments on the change of the flow 
régime around a sphere caused by the application of a 
thin metal ring on the surface of the sphere and that of 
the intentional increase of the flow turbulence by plac- 
ing a metal net before the sphere, published in 1914, 
are certainly convincing. But those of the author with 
square plates, published in 1912 and recalled here above 
(Figs. 1, 2, 3), are equally so, and one ought to take into 
consideration that, when Prandtl carried on his experi- 
ments, he was already acquainted’? with the results 
obtained by the author. Thus, it would be only fair 
to modify Dryden’s assertion! that ‘‘the discovery of 
wind-tunnel turbulence effects began with spheres and 
a similar effect of turbulence on maximum lift coeffi- 
cients was found subsequently,” as follows: “The 
discovery of wind-tunnel turbulence effects began 
with square plates and a similar influence of turbulence 
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on drag coefficients of spheres was found subse 
quently.” 
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